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On Korn’s Inequality 


L. E. PAYNE & H. F. WEINBERGER 


Communicated by R. A. TouPIN 


1. Introduction 


Let %(x,, 2, x3) be a continuously differentiable vector field on a region R 
with boundary B. If an isotropic homogeneous elastic medium in R is subjected 
to an infinitesimal displacement proportional to #, the classical strain energy is 
proportional to 


(1.1) E(t) =S (#) + [[f laivaeax 
where : 
(1.2) ee (4,5; + 4,3) (wu ij + M7, i) 4% 


Here u,, U2, us are the components of # in the rectangular coordinates (x,, %5, %3), 
ax is the element of volume, the constant o is Poisson’s ratio, and the symbol , ; 
indicates differentiation with respect to x, We use the summation convention 
throughout. 


We introduce the Dirichlet norm 


(1.3) Di) ah J eta naa 

Korn’s inequality asserts that if # satisfies the integral conditions 
(1.4) dBA idee = ON cee tongs 

if 

(1.5) C=, 


and if R is sufficiently regular, there exists a constant K, depending only upon 
o and RF such that 


(1.6) D(é) SK, E(u). 
It is easily seen that 
(4:7) OS fff |divé|?dx<3S(d). 
R 


Thus if (1.6) holds for o=0, 2.e., if 
(1.8) D(é)< K,S(u%), 
then (1.6) clearly holds for all o> —1, with Ky)=K, for o20 and Ky=(1 +0) K,/ 


({—20) for —1<o0<0. (For negative values of o we have made use of (1.7).) 
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The inequality (1.6) is useful in proving the existence and uniqueness of 
the solution of the second boundary-value problem of classical elasticity. An 
explicit knowledge of the constant Ky also leads to explicit arbitrarily close 
bounds for the solution. This follows from results of [7] if B is star-shaped 
with respect some point in R, and from [8] if R hes exterior to B. 


We now let 
(1.9) R (uz) =Z JJ (U4; ; = Us; i) (4; 5 Se U; i) dx. 
It is easily shown that the inequalities (1.8), 
(1.10) R(é) < (4 — Ky’) D(a) 
and 
(Ae) R(#) < (K,—“) S(u) 


are completely equivalent. Korn [4, 5], FRIEpDRICHS [3], and BERNSTEIN & 
TouPIn [1] define the constant (K,—1) of (1.11) to be the Korn constant K. 
It is clear that K,=1. 

Korn’s inequality was first proved by A. Korn [4, 5] by means of integral 
equations on the boundary. A more direct proof involving integration over the 
interior of R was given by K. O. FriEpRIcHs [3]. In a recent paper in this 
journal, B. BerNsTEIN & R. Toupin [Z] gave a simplified version of the FRIED- 
RICHS proof which permitted them to prove the inequality (1.8) with the explicit 
value 
(1.12) K,=21 
for the case when R is a sphere*. 


In this paper we reverse the procedure of BERNSTEIN & TouPiIN. We first 
prove the inequality (1.8) for the sphere by explicitly determining the spectrum 
of the quadratic form D(«%) with respect to S (#), thus showing that it is bounded. 
In this way we find the best possible constant 


(1.13) Ky= 3 


in (1.8) for the case of the sphere. We also find the best value of Ky for any 
o> —1 in the inequality (1.6) for the sphere. 

In Section 3 we extend Korn’s inequality to any region R that can be 
mapped onto the sphere by a twice differentiable mapping. 

In Section 4 we show that if Korn’s inequality holds on several domains, 
it also holds in their union. In this way we obtain Korn’s inequality on a rather 
general class of regions R. 


In Section 5 we obtain for the two-dimensional circular region the best 
possible constant K 
= 4. 


2. The sphere 


Let & be the unit sphere centered at the origin. Let % be a continuously 
differentiable vector field satisfying the normalization conditions (423)s “teas 


* It is easily seen that (1.8) is invariant under a uniform dilatation. Thus, the 
constant fx, for a sphere is independent of its radius. 
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easily seen that under these conditions both S(#) and D(u%) are positive definite, 
_ providing we identify vector fields differing only by a constant. We seek the 
spectrum of D(w) with respect to S(#). The inequality (1.8) asserts that this 
_ spectrum is bounded. It follows immediately from integration by parts that if 


2.1) : div@=0 inR 
and if 

f2.2) =O) fava g 

then 

(2.3) DV OS (ay. 


Thus, 2 is an eigenvalue of infinite multiplicity. 
Also, if y is any twice differentiable function, then 


(2.4) D (grad y) = S (grad y). 


Thus, 1 is also an eigenvalue of infinite multiplicity. 


To obtain the remainder of the spectrum, we derive the Euler equation in 
the usual manner. We obtain the differential equation 


(2.5) (2—K)u,,,—Ku,;,,=0 inR; (aie As) 

and the boundary conditions 

(2.6) [((2— K)u,;— Ku, ,)*,;=0 on B; Sebo 

Here K is the eigenvalue parameter. We have used the fact that the radius 


vector ” with components 4; is a normal vector on B. 
Applying the operator 7- curl to (2.5), we find for K +2 


(2.7) A(?-curl%) =0. 
Since the operator 7 - curl involves only tangential differentiation on the sphere, 
we can apply it to (2.6) to obtain 


(2.8) (2—K)r 2 curl) + K? curl =o on B. 


It follows from (2.7) that the left-hand side of (2.8) is harmonic in R. Hence 
it vanishes identically. Thus, 7- curl # is a regular harmonic function which is 
homogeneous of degree — A(2— 4K)? in R. Either it vanishes identically, or 
this degree must be a positive integer. (Since 7- curl w= }div[(r?—1) curl w], 
its integral vanishes, so that 7- curl cannot be a constant other than zero.) 


We thus find the eigenvalues 


(2.9) 1 ee le a es 


nm—1’ 
with the corresponding eigenvectors 
(2.10) u =? X grad P™, 
Here P™ is any spherical harmonic of degree . That is, it is a homogeneous 
polynomial in the x, of degree m and is harmonic. The function (2.10) is easily 
seen to satisfy 7- curl #=—n(n+1) P” as well as (2.5) and (2.6). 
7 
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Note that the case m=1 in (2.9) is eliminated by condition (1.4). 


If we have any function of the form 


(2.11) u =? xgrady 
where 

(2.42) An=0, 

we may expand 7 in spherical harmonics: 

(2.13) ” ~Da, Pe. 
Normalizing the P™ so that ea 

(2.14) DF x grad PO”) = 4 


we find that if D(7 x grad 7) exists, 


D(? Xgrady) = > |a,|? 
(2.15) vee 


S (# x grad 7) = 4(5 = *) | ho 


If % is subjected to the condition (1.4), we hav 


(2.16) Havel cenanewenee— [f| Haar I ==Oy 


which implies that a;=0. Consequently, we find that 


(2-47) D(¥ x gradn) S$ 4S (7 x gradn). 


Finally we consider the eigenfunctions of (2.5), (2.6) for which 7- curl #=0. 
This means that # can be expressed as the sum of a gradient and a multiple 
of 7. Taking the divergence and curl of (2.5) shows that if K 1, 2 both diva 
and curl # are harmonic. It follows that # can be expressed in terms of three 
harmonic functions g, @ and t by 


(2.18) “@ = grad[y + 720] + (r2?—1) gradt, 
(2.19) Ag=Ao=Ar-=0. 


Inserting this result into (2.5), we find that (aside from an unessential additive 
constant) 


a) 0 , 
(2.20)  a(1—K) (ar +3)e+|(4—3K)r 2 42 K|r=0. 
The normal component of (2.6) gives, if K 1, 
a) 0 o 4) OT 
(2.21) rolr a thot (rs +2)(r= +1)e+2r3% =0 


on B. Since the left-hand side is harmonic, (2.21) holds in R. Similarly we find 
from the tangential component of (2.6) that 


(2.22) (1—K)(r = 1} L (1 K)(r ae 1)o ye eae 


iby, JRC. 
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Since the equations (2.20), (2.21), (2.22) involve only the differential operator 
- the harmonic functions g, 9 and t are multiples of a spherical harmonic P™. 


We find the eigenvalues 


__ 2(n+1)(2n+1) = 
(2.23) us a areas gees yy aa aS 


with the corresponding eigenvectors 


(2.24) “= grad [n{2n+3 —972(2n+1)} P] + (72 —1) grad [(3 n?+5n+1) P™]. 


Any linear combination @ of these eigenvectors can be written in terms of 
a harmonic function 9 as 


wm = grad 


(2.25) (2+ (17°) (2 i r i)}r 7 I 


+ (7? — 1) grad |(3 (r i ie Oy. 2 1) 9]. 


Conversely, if w is of the form (2.25) and D(#) exists, we may expand 9 in 
spherical harmonics 


(2.26) o~ Ss pad ee 
Then ie 
S (B) = Yi a2 SH) 
(2.27) ee a 
D(é) = >) a, DU) = > a, 2h 1) (a), 


if 
where #) denotes the eigenvector (2.24). Both series in (2.27) converge. Noting 
that the largest of the numbers (2.23) is given when n=3, we find that for w 
of the form (2.25) 


(2.28) Dw) S28 Siw). 

Now if # can be expressed in the form 
(2.29) “#=v+gradp+7xgrady +H 
where 


divv=0O inR, 
(2.30) p= 0. 20n-B, 


Aye in R, 
and w is given by (2.25) with 
(2.31) A o==0- in KR; 
it is easily seen by integration by parts that 
S (a) = S(#) + S(grad yp) + S(7 x grad) + S(®), 
222) D(%) =D(s) + D(grady) + D(x grady) + D(#). 


That is, the four components are orthogonal in the sense of the norms [S (#) ]# 
and [D (#)]?. 
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It is easily verified that v, grady, and w satisfy the condition (1.4) auto- 
matically. Hence (1.4) implies (2.16). 
Thus, it follows from (2.3), (2.4), (2.17), and (2.28) that 


(2.33) D (uv) = 73 S (a) 


for all % of the form (2.29). We shall show that any sufficiently differentiable 
vector field % is expressible in the form (2.29). 

Applying the operator 7-curl to (2.29) on the boundary B, we find the 
boundary condition 


DET a ee ee 
(2.34) hae (r al 1) 7-curlu. 


Since the left-hand side is harmonic, it (and hence 7) is determined by this 
boundary condition. 

Since 7-curl[“—7xgrady] vanishes on B, the tangential components of 
u—?xgrady coincide with the tangential derivatives of a function f on B. 
That is, there exist functions f and g such that 


(2.35) ué —7Xgradyn=gradf+gr7r on B. 
Using the form (2.25) of w and the fact that v vanishes on B, we find that 
veptt—mfrd +S 
030), Bee y 
ate 2+ (1 P) (27 a) aes eee tg onB 


Taking the Laplacian of the divergence of (2.29) we find that 
nD o oe) a 3 
(2.37) AAlp+ [2+ 9) (27 +1) 7 22h — A div ii. 


The differential equation (2.37) together with the boundary conditions (2.36) 
serves to determine the function in braces. Taking the divergence of (2.29), 
we find 


(2:38) 2(3(r5) +572 + 1)r 22 <div —Aly+ [2+ (1—r9) (ar 2 +4)] 7 22. 


. Since the right-hand side is already determined and is harmonic by (2.37), 
this equation determines the harmonic function @. This then determines y. The 
he field v is now obtained by transposing in (2.29). It automatically satisfies 

DO), 

Thus we have shown that (2.33) holds for sufficiently differentiable vector 
fields 7. Hence it holds for limits of such fields in the sense of the norm [D (uv) }°. 
These are precisely the vector fields with Dirichlet integrable components. 


The inequality (2.33) is sharp in the sense that equality holds for # of the 
form (2.24) with n=3. 


Remark. In the same manner we can obtain a sharp inequality 


(2.39) Du) < Ky El) 


| 
i 


a 
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with E(w) defined by (1.1) for any o> —1. We find that 


4, (pa 
an 
ENT 2s eee eae te 
(2:40) Ky = 20g a ise ee Een lus ], 235 =0>54 
| in + [37-200 + {529 4000+ 160007}! Sse ees. 
14+100 ! , 4 29335 
1 F 
Pepe i ot Oe |, 1A<oS— 


The different functions of o correspond to different types of SEN TCC OR for which 
the maximum eigenvalue is attained. In particular for ¢= ;’; equality is attained 
with 4 =7 x grad P®, 


3. Smooth simply connected regions 
et 
(3.4) y= 9" (1, Xe, X) 
map the given region R onto the unit sphere y'y'<1 when the x; and y, are 
considered as Cartesian coordinates. Let the functions y’ have bounded first 


and second derivatives, and let the Jacobian determinant be bounded away 
from zero. Then we can consider the y’ as curvilinear coordinates on R. 


If the vector field % has the components u; in the y-coordinates, 


S(#)= aSSJ Veg" e” (U1; + Uys) (Mays + Mx) 22 dy? d 8 


(322) a a ee 
D(ué) = JJ Vege! us, Up d yd y? dy? 
where 
Ou; 
(3.3) Uj; = = Sah us 


oy) 
denotes the covariant derivative, and the I;# are the Christoffel symbols [7]. 
The conjugate metric tensor g*? is defined by 


e ay? oy) 
V7 = 7 
B.4) ax? axP’ 


and g is the reciprocal of its determinant. 


We let 
Ou. Ou; \ [ OU; Ou. 
ag Fe 4 i 4 8 \ayidy2 dy, 
=e UU) eget : calle ae Es ” 


(3.5) Ou; OU; 
4 yld 2d 3 
oo oyI Aree ba y . 


the integrals being over the region y' y'<1. We have shown that if w; satisfies 
the conditions 


ra) Cae 
(3.6) Res (= id = )aytdyrdy? = i,j =1, 2,3, 


yi yi<l 


then 
(3-7) Dy (%) S73 So(%). 
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Since the mapping (3.1) is assumed to be non-singular and to have bounded 
second derivatives, there are positive constants «, f, and y such that 


(3.8) ab F< gigtEE SBE, E; 
and eae. 
(3.9) Ve ee” Tf The, 6, S77 86: 


for all vectors &, and all points of R. It then follows with the aid of the triangle 
inequality that ‘ 
[SPS alS@l—y| [fuu ay ayray'), 
yysl 
10) i \ 
: [D (VSP Doli) +y| Sffurmadydytays). 
a 


It is well known [2] that the quadratic form 
(3-11) SSf pay dy? dy? 
yt yts1 


is completely continuous with respect to the Dirichlet integral 


é re) 
(3.42) ff a gyi dy ay? dy’. 


ye ytS1 
This means that a finite number of linear functionals L,(y), v=1,..., N can 
be found such that the integral of gy? is less than any positive constant «? times 
the integral (3.12) whenever L,(y)=0, v=1,..., N. If we choose 


(3.13) e<alse]? 
and impose these constraints on all the components w;, and if in addition the u; 
satisfy (3.6), we have from (3.7) and (3.10) 


56(B+ ey)? - 

(3.14) DN iasyeeae 

Thus D(z) is bounded with respect to S(#) on the orthogonal complement 
of a finite-dimensional subspace A of the Hilbert space of vector fields satisfying 
(4.4) with norm [S(z#)]*. But D(#) is also bounded on the finite-dimensional 
space A. Otherwise there would be a vector field # with Dirichlet integrable 
components satisfying (1.4) for which S(#)=0, D(w)=1. But S(#%)=0 implies 
wu = constant and hence D(w)=0. Thus, there is a constant K, such that 


(3.15) D (®) < K, S(w) 


for @ in A. Denoting the constant appearing on the right of (3.14) by K,, we 
decompose any % into a vector @ in the orthogonal complement of A and @ in A. 
Then we find, using ScHwarz’s inequality, 


D(i) =D@+B) (14+ 34) DG) +(t+32)D@) 
(3.16) <(K,+K,)[S@)+S@)] 
= (K, + K,) S(#). 


Thus, a Korn inequality holds for any domain R that can be mapped onto 
a sphere by a non-singular mapping with bounded second derivatives. 
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4. More general regions 
Let R be the union of regions R,, R,...,R, on which Korn’s inequality 
holds. If % is subjected to the conditions (1.4) for each R,, we have, in the 
obvious notation 
Hence 


(4-1) Dé) SDD, (%) S DK, S,(u) SDK, S (di). 


Thus, D(w) is bounded with respect to S(%) when % is subjected to a finite set 
of linear homogeneous constraints. By the same argument as that used in the 
preceding section we then show that Korn’s inequality holds on R. 

Thus, we have shown that Korn’s inequality holds in the union of a finite 
number of domains, each of which can be mapped on a sphere in a twice differen- 
tiable non-singular manner. 

A sufficient condition for this is that the boundary B is covered by open 
subsets, on each of which one can introduce a coordinate system with a diffe- 
rentiable metric tensor, and a differentiable vector field whose angle with any 
inward-pointing vector is bounded away from zero. We map a neighborhood of 
each open set onto a slab by projecting along the vector field, using the coordi- 
nates on B and the distance from B as coordinates. We approximate these slabs 
by subdomains with smooth boundaries in such a way that the union of these 
subdomains covers a neighborhood in R of the whole boundary B. We can 
clearly cover the remainder of R with smooth simply connected regions. Thus, 
R is the union of a finite set of regions on which Korn’s inequality holds. 

Our condition on the vector field permits R to have edges and corners, but 
not cusps or cusp-like edges. 


5. The circle 


We consider now the analogous two-dimensional problem of finding the 
optimum constant K, in (1.8) where D(«#) and S (z#) are integrals over the interior 
of a circle, and 7 andj take on the values 1, 2. The vector # is again subjected 
to condition (1.4). 

In a manner similar to that used in deriving the Korn constant for the sphere, 
it can be shown that any vector # which satisfies (1.4) can be decomposed in 
the following manner: 


(5.1) “=v+grady+w 

where 

(5.2) div¥ =0 

in R and 

(5-3) v=0 

on C. Also @ is expressible as 

(5.4) @ = 2|r2 grad p— 277 - | (v2 — 4) grad (e+r52), 


where gy is a harmonic function. It is easily verified that v, grad yp and w are 
mutually orthogonal in both [D(#)]* and [S(#)]* norms. 
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As in §3 we note that for functions ¥ satisfying (5.2) and (5.3) 


(5.5) D(@) =2S(v). 
Also as before 

(5.0) D (grad yp) = S (grad y). 
It is easily verified that 

(5.7) D () = 45 (). 


We thus establish the curious fact that the spectrum of D(#) with respect to 
S(w%) under condition (1.4) consists of only three eigenvalues (A=1, 2, 4), each 
occurring with infinite multiplicity. We find the sharp inequality 


(5-8) D(é) S4S(u), 


so that the best K, for the circle is 4. BERNSTEIN & TouPIN [J] showed that 
Kee 
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Theorems in Linear Elastostatics for Exterior Domains 


M. E. GURTIN & ELI STERNBERG 


1. Introduction 


The conventional proofs of fundamental theorems in three-dimensional clas- 
sical elastostatics, such as the traditional proofs of the uniqueness and reciprocal 
theorems or of the minimum energy principles!, are confined in their validity 
to bounded domains, although this limitation is not always made explicit in 
the treatise literature. The present paper aims at generalizations of the fore- 
going theorems to unbounded domains exterior to a finite number of closed 
surfaces. 


The generalizations to which we have alluded are elementary in the presence 
of sufficiently stringent restrictions upon the behavior at infinity of the relevant 
elastostatic field quantities. Thus KircHHorr’s [4] uniqueness proof at once 
carries over to exterior domains provided the cartesian components of displace- 
ment and stress, #; and t;;, obey the conditions 


u,;(x)=c;+O0(r), 45 (%) =¢;; + O(r*) as 4—>0o, (is 


where 7 is the distance from the origin, c; and c;,; (c;;=c;;) are prescribed real 
constants, and the notion of “order of magnitude’ is used in its standard 
mathematical connotation‘. 


While requirements (1.1), together with appropriate smoothness hypotheses, 
are sufficient to assure the uniqueness of the solution to the fundamental 
boundary-value problems of elastostatics for exterior domains, such a priori 
assumptions as to the characteristic orders of magnitude at infinity are quite 
artificial: the vate at which the displacements and stresses approach their limiting 
values at infinity is an item of information which one would legitimately expect 
to infer from the solution to the problem, rather than a condition to be imposed 
on the solution in advance. By the same token, a uniqueness theorem resting 
on (1.1) leaves in doubt the status of solutions to exterior boundary-value problems 


1 See, for example, Love [1], SoKoLnikorF [2], TIMOSHENKO & GOODIER [3]. 
Numbers in brackets refer to the list of publications at the end of this paper. 

2 Infinite regions whose boundary extends to infinity are excluded from the 
following considerations. 

3 The single argument ¥ denotes the triplet of cartesian coordinates ¥;. Through- 
out this paper subscripts have the range of the integers (1, 2, 3); summation over 
repeated subscripts is implied and subscripts preceded by a comma indicate differ- 
entiation with respect to the corresponding cartesian coordinate. 

4 See Section 2. 
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since it does not preclude the existence of alternative solutions to the “same | 
problems”, whose displacements and stresses approach their values at infinity 
less rapidly; nor could such solutions be rejected on convincing physical grounds. 


As far as the first boundary-value problem (given surface displacements) is 
concerned, the preceding uniqueness issue was resolved in a recent paper by 
Durrin & Nor [5]. In this instance, as shown in [5]*, conditions (1.1) may, 
without loss in uniqueness, be replaced by the mere requirement that the 
displacements uniformly tend to given values as 7— oo. 

The method of attack employed by Durrin & NOLL is well suited to their 
specific purpose but is not readily adaptable to a similar extension of the un- 
iqueness theorem if the stresses are prescribed at infinity and in case the boundary 
conditions involve given surface tractions or are of the mixed type. Also, analog- 
ous generalizations to exterior domains of other basic theorems in elastostatics 
would appear to call for a different approach. 

In Section 5 of the present paper we prove two theorems concerning the 
behavior of elastostatic fields at infinity, which enable us to meet the broader 
objectives just mentioned. A specialization of the first of these theorems yields 
the conclusion that if the displacements tend to zero uniformly at infinity, the 
displacements and stresses necessarily obey (1.1) with c;=c;;—0, provided the 
body forces are suitably restricted. Subject to a similar limitation of the body 
forces, one infers from the second theorem that the uniform vanishing of the 
stvesses at infinity implies that the stresses and displacements, except for an 
additive rigid displacement field, are again governed by (1.1) with c;=c;;=0. 

The material preceding Section 5 is preliminary to the two theorems cited 
and is partly expository in character. In Section 2 we collect certain facts 
pertaining to solid and surface harmonics, and recall a known theorem on the 
behavior of functions harmonic in a neighborhood of infinity. A corresponding 
theorem for biharmonic functions is deduced in Section 3. In Section 4 we recall 
certain properties of elastostatic fields and deal with theorems concerning the 
representation of such fields in a neighborhood of infinity. 


Finally, in Section 6, we extend various basic theorems of classical elasto- 
statics to exterior domains by utilizing the results established in Section 5. 
Specifically, we give appropriate generalizations of the reciprocal theorem, the 
principle of work and strain energy, the uniqueness theorem, and the two 
minimum energy principles. The present extension of the uniqueness theorem 
assures the uniqueness of the stresses belonging to a sufficiently smooth solution 
of the general mixed boundary-value problem for an exterior domain if either 
(a) the displacements or (b) the stresses tend uniformly to pre-assigned constant 
values at infinity (regardless of the rate at which these limiting values are 
approached), provided in case (b) the resultant force and moment of the surface 
tractions on the boundary of the medium are also prescribed. This supplementary 
condition may, of course, be omitted if the surface tractions are given on the 
entire boundary. An analogous uniqueness theorem has long been known to 
hold in two-dimensional elastostatics. 


* See also [6], which contains a slightly more compact version of the proof 
presented in [5]. 
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2. Preliminary definitions. Behavior of harmonic functions near infinity 
Let (7,8, ~) be spherical coordinates, related to the rectangular cartesian 


coordinates x; through the mapping 


) 
| 


| 


Xx, =7sin? cos @, %_=rsindsin y, A= 1 COSY, 


2.1) 
OEE Has ost Oa =r; OSG 2m: ( 


A deleted neighborhood of infinity, henceforth denoted by @, shall be a region 
characterized by 77<7< oo, and thus bounded internally by a sphere of radius 7. 
We note that a function of position defined in &need not possess a limit as 7—>oo. 

If v(7,%,q@) is a scalar field defined in @, we write v(r, 3, y)=O(r") or 
u(r, 3, p)=o0(r") according as |r~-"v(r, #, y)| remains bounded (uniformly in 9 


and qg) or 7~"v(r, 0, y) tends to zero (uniformly in @ and g), when r->co. A 
function v(r, 3, y) defined and continuous together with all its partialderivatives 


of order m (n=O) in some region of space, will be referred to as being of class 


C) in that region. If u(r, 9, y) is of class C and satisfies Laplace’s equation 


V2v=0 in an open region R, it will be called harmonic in R; if v(r, 0, ~) is of 
class C and meets V4v=0 in R, it will be called biharmonic in R. 

We now summarize certain well known results in the theory of harmonic 
functions®, which will be needed later on. In this connection we recall first 
that the general solid spherical harmonic of degree k admits the representation 


(7, 3, gy) = S” (8, gp) #* (R= 0,-E4, +2, .:.), (2.2) 


where S) (9, pm) is the general surface spherical harmonic of the same degree, 
and is given by 


|r| 
S (9, v) => [a cosn gy + b™ sin n gy] Pi” (cos ¥), OO =O. (2.3) 
n=0 


Here Pl” designates the associated Legendre function of the first kind, of 
degree k and order n, while a, a”, b” (n=1, 2, ..., | R|), for fixed k, are 2k+1 
arbitrary constants. If k and are non-negative integers, 
5 a’ P(E 4 hig &2__4 k 
PME =1-ee SE, A= Se (2.4) 


2k! dé 


in which P, is the Legendre polynomial of degree k. Consequently 
PO(E=P 6), BPE)=0 (n>k2O0). (2.5) 
In addition we have the recursion relations 
Pec) bye (e) (2.6) 


which are valid for unrestricted k and n. As is evident from (2.3), (2.5), (2.6), 
every surface harmonic of degree —1—k is a member of the family of surface 
harmonics of degree k. Next we recall that any two surface harmonics of distinct 
degree are orthogonal over the unit sphere y=1. With the notation 


2a 7% 


Sy Sr) i) neat ™ (9, py) sinddddy (2-7) 


6 See, for example, Poincar [7], KELLoac [8], and Hopson [9]. 
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we thus have 


DS nS =O (Rm, k= —m— 1). (2.8) 


The solid harmonics H") (7, 0, p) are harmonic functions in &, which are homo- 
geneous of degree k with respect to the cartesian coordinates (%,, Xg, X3) 5 they 
are homogeneous polynomials of degree k in the x; if & is a non-negative integer. 
Further, if H” (rv, #, p) is a solid harmonic of degree k, 


HY, 8,9) =H°7,0,9) (k=0, £1, £2,...), (2.9) 


where H (7, 9, y) is a solid harmonic of degree k—1. Also any specific solid 
harmonic of degree k satisfies the identity 


x, H™ (7,0, 9) =kH™ (7, 0, 9) (k=O, +1, +2,...). (2.10) 


Turning from solid spherical harmonics to general harmonic functions, we cite 
the following theorem. 


Theorem 2.1. Let H(r, 3, y) be harmonic in &. Then 
(a) H(r, 3, p) admits the representation 


loc) 


H(r, 9, y) = 2 H (r, 8, 9), (2.11) 


k=—00 
where the H™ (r, 8, p) are uniquely determined solid spherical harmonics of degree 
k and the infinite series 1s uniformly convergent in every closed subregion of 6; 
(b) H(r, 8, yp) mm & has partial derivatives of all orders, series representations 
of which may be obtained by performing the corresponding termwise differentiations 
im (2.11), the resulting expansions being also uniformly convergent in every closed 
subregion of &; 
(c) of n ts a fixed integer, the three statements 
(2) H(7, 8,9) =O(r""), 
(B) H(r, 8, 9) =o0(7"), 
(vy) H® (7,8, p)=0 in (2.41) for R=n, 
ave equivalent and imply 


(0) 1 .(7,8,9) =O"). 


3. Behavior of biharmonic functions near infinity 
In this section we state and prove a theorem which is the counterpart for 
biharmonic functions of Theorem 2.1. 
Theorem 3.1. Let F(r, 0, py) be btharmonic in & Then 
(a) F(r, 0, yp) admits the representation 


F(r,9,9) = > 10,9,9) + YH (7,8, 9), 6.1) 


k=—oo k=—oo 


where oe (r, a, ?) and Hy, 9, p) are solid spherical harmonics of degree k and 
both infinite series are uniformly convergent in every closed subregion of €; 
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(b) Flr, 9, p) 1 & has partial derivatives of all orders, series representations 


of which may be obtained by performing the corresponding termwise differentiations 


m (3.1), the resulting expansions being also uniformly convergent in every closed 


subregion of 6; 


(c) af nts a fixed integer, the three statements 
(x) F(r, 3,9) =O(r"™), 
(6) F(r, 8, ¢) =o(r"), 
(y) h® (7,3, p) =H") (7,9, ~) =0 in (3.1) fork=n, 


ave equivalent and imply 


(6) F;(,8,9) =O"). 


Proceeding to the proof of the foregoing theorem, we show first that any 
function F(7, 3, y) which is biharmonic in &, admits ALMANSI’s representation 


Ely, 3, 9) =h(r, 9, p) +P Hr, 9, @), (3.2) 


where h(7, #, p) and A(z, %, m) are functions harmonic in &. A completeness 
proof for the representation (3.2), applicable to a region that is star-shaped 
with respect to the origin and contains the origin in its interior, is indicated in 
[10|". To establish the completeness of (3.2) in the present circumstances, it 
is evidently sufficient to exhibit a function H(7, 9, y) which is of class C®) in 
& and there meets 

Var Oss ee Veit), (3.3) 
or, equivalently, 


Vapor 4r E+ OH = VF. (3.4) 


Since, by hypothesis, V?F is harmonic in &, we know from Theorem 2.1 and (2.2) 
that it admits the expansion 
Pe ¥ S® (0,9) (3.5) 
k=—0oo 
where the S®) (9, p) are surface harmonics of degree k and the infinite series 
has the convergence properties asserted in parts (a), (b) of Theorem 2.1. Now, 
consider H(r, #, y) defined by 
SWB, 9) | 
xo) 6 
H(r, 3, @) => 5% SCO aaa (3.6) 
Clearly, the infinite series in (3.6) has the same convergence properties as that 
in (3.5). Hence the function H(r, 0, g) given by (3.6) is harmonic in @; further, 
this function also satisfies the second of (3.4), as is confirmed by direct sub- 
stitution from (3.6) and by recourse to (3.5). 

On applying (a) and (b) of Theorem 2.1 to # and H in (3.2), (a) and (b) of 
the present theorem follow at once. With a view toward establishing (c), we 
observe that (f) is immediate from («) by virtue of the meaning of the symbols 
“QC” and ‘“‘o’”’. We show next that (f) implies (y). Let s(#, y) and S® (8, ¢) 
be the respective surface harmonics corresponding to the solid harmonics 


7 See Frank & v. Mises [10], p. 848 e¢ seq. 
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hn (r, 0, ) and H" (r, , y) entering the infinite series in (3.1). We now multiply 
the right-hand member of (3.1) by 77s (9, g) (m=0, +1, £2, .--) and inte- 
grate termwise over OS 057, OS S22. Bearing in mind (2.2), the convergence 
properties asserted in (a), as well as the orthogonality relations (2.8), and using 
the notation (2.7), we find in this manner that the integral 


2a 7 


Inn) =r" S SFr, 0, 9) 8 (8, p) sn ddd dg, (3.7) 
0 0 


for all sufficiently large values of 7, is given by 


7 (7) —yn—n [s™), si] ae y nmi Ee si] + 


= Nt 2 Wigs s™] abe yo n—mrd emer s(™) d 


mn 


(3.8) 


It is apparent from (3.7), because of (8) and since s™(#, m) is bounded on 
0<0<a, 0S ySX2z, that I,,,,(7) tends to zero as y->oo for fixed m and n; 
consequently the same is true of the right-hand member of (3.8). Noting that 
the coefficients of the inner products appearing in (3.8) are four distinct powers 
of r, we conclude that [s™), s™]—=0 if m—n=0. Recalling that s™ (9, y) is 
real-valued and taking account of (2.2), (2.7), we thus confirm that s™ (3, y) = 
h” (r, 9, p) =0 for m=n. On multiplying the right-hand member of (3.1) by 
v-" S™) (9, p) and proceeding as before, we verify also that H™ (7, &, y) =0 for 
m=n—2. Hence (f) is a sufficient condition for (y). 

We have to show further that (y) implies («). This is readily seen to be true 
by observing that each of the series in (3.1), according to (a) and (b), represents 
a function harmonic in @, and by invoking part (c) of Theorem 2.1. Finally, 
(vy) implies (6) because of (b) of the present theorem and (c) of Theorem 2.1. 
Theorem 3.1 has thus been proved in its entirety. 


4. Representation of elastostatic fields in a neighborhood of infinity 


The complete system of fundamental field equations in the linear equilibrium 
theory of homogeneous isotropic elastic solids, with reference to rectangular 
cartesian coordinates x; and in indicial notation’, takes the form 


Oo 
Tip +h, =0. (4.3) 


Here u;, e;;, and t;; are the cartesian components of displacement, strain, and 
stress, f; denotes the components of the body-force density, «4 and o stand for 
the shear modulus and Poisson’s ratio, whereas 6;; designates the Kronecker 
delta. Equations (4.1), (4.2), and (4.3) represent respectively the displacement- 
strain relations, the stress-strain law, and the stress equations of equilibrium. 
Necessary and sufficient for the positive definiteness of the strain-energy density 
are 


Ore ag Pe Gees (4.4) 


8 See Footnote No. 3. 
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_ and, provided (4.4) hold, the relations (4.2) may be inverted to yield 


peal 
2 


oO 
C3; Vij 1-to 07) Tp» os (4.5) 
Also, if «;(«) is of class C®) in the region at hand, elimination of the stresses 
and strains among (4:1), (4.2), and (4.3) leads to the displacement equations of 
equilibrium 


1 1 

i NIRary cues CMa aye mee (4.6) 

For future convenience we introduce the following definition of a “simple 
body-force field’. 


Definition 4.1. A body force field f;(x) is simple in an open region R if f;(x) 
ts of class C® in R and 
hii; Exgelp, 7 =O mM Ke (4.7) 


In (4.7) €;;, designates the components of the usual alternator, so that /;(x) 
here is required to be both solenoidal and irrotational in R. From (4.7) and 
the smoothness of /;(x) follows 

27 ==) as (4.8) 


for a body-force field which is simple in R; thus /;(x), in the present instance, 
is harmonic in RF and accordingly has continuous partial derivatives of all orders 
in that region. We are now in a position to state concisely the subsequent known 
theorem of elastostatics. 


Theorem 4.1. Let u;(x) be of class C® in an open region of space R. Let 
U;(x), €;;(%), and t;;(x) 1m RF satisfy (4.1), (4.2), (4.3), subject to (4.4), the body- 
force field f;(x) being simple in R. Then 

(a) the functions u;(x), e;; (x), T;;(%) have continuous partial derivatwwes of all 
orders in R; 

(b) V2u;,;=0, ej pV2u, 5 =0; (4.9) 


(c) Viu;=Vie,,=Vit;;=0 mR. (4.10) 


Accordingly, any solution of the fundamental field equations of classical 
elastostatics which is of class C®) and corresponds to simple body-forces, is 
biharmonic, while the associated dilatation and the rotation vector are harmonic 
in the region under consideration. Part (a) of the theorem just cited was proved 
by DurFin [11] for the special case in which the body-force field vanishes identi- 
cally; the argument indicated in [//] is, however, readily generalized to accom- 
modate any simple body-force field. Once (a) has been established, (b) and (c) 
follow® trivially with the aid of (4.6), (4.7), (4.8), (4.1), and (4.2). 

We turn next to two representation theorems for solutions of the elasto- 
static field equations in a neighborhood of infinity. The first of these theorems 
concerns the general solution of (4.6) in the absence of body forces; the second 
aims at a particular solution of (4.6) corresponding to arbitrary simple body 
forces. 


® See, for example, [2], p. 77. 
Arch. Rationai Mech. Anal., Vol. 8 8 
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Theorem 4.2. Let u;(x) be of class C‘?) in & and there satisfy (4.6) with f(x) =0, 
the elastic constants being subject to (4.4). Then u;(r, 9, p) admits the represen- 


tation 
foe) 


u, (7,9, 9) = WY, 8,9) +? X AY (7,4, 9) (4.11) 
k 


k=--co =—oo 


where h® (r, 8, p) and H\ (r, 9, y) are solid spherical harmonics of degree k and 


( = : (k+ 2) ==; eeele . 
HY = EES EMC CSOT ChE he (k=0, +1, +2,...) (4.12) 
Both infinite series in (4.11), as well as the series resulting from any finite number 
of termwise differentiations of (4.11), are uniformly convergent in every closed sub- 
region of &. 

A somewhat less explicit version of this theorem was established by KELVIN 
[12] in connection with his work on elastostatic boundary-value problems for 
a spherical shell. The following proof, however, would appear to be more direct. 

By hypothesis and Theorem 4.1, u;(7, 8, y) is biharmonic in &. Theorem (3.1) 
therefore assures the existence of an expansion in the form (4.11) which has 
the asserted convergence properties. 

This leaves only (4.12) yet to be confirmed. From (4.11), (2.10), and the 
harmonicity of h™, H™ follows 

Vu, = 4, 7=2 2 (2k + 3) HY. (4.13) 


k=—oo 
Substitution from (4.13) into (4.9), in turn, leads to 


foe) 


>», (2k +3) AY) =0, : > (2k +3) £;;p 19), = in &. (4.14) 
k=—0o =—0o 


But H!"), according to (2.9), for fixed 7 andj, is a solid harmonic of degree k—1. 
Hence (4.14), by virtue of the uniqueness of the expansion (2.11) in Theorem 2.1, 
imply 

TO) ha (R=O +1, +2,...). (4.15) 


’ J J,% 


With the aid of (4.11), (2.10), and (4.15), we obtain 
mente: me (Ai); + 2(k +1) A), (4.16) 


and inserting (4.13), (4.16) in (4.6) with /;—0, there results 


wy AiO, + 2[k + 1+ (4 — 20) (2k + 3)] HP} =0. (4.17) 


Since /\"), is a solid harmonic of degree k — 2, (4.12) is a consequence ™ of 
(4.17) and Theorem 2.1. The proof of Theorem 4.2 is now complete. Expansions 


0 If g(%)=8(4%,, %,, ¥3) is a function of the cartesian coordinates %;, we write 


consistently g(v, #, y) in place of g[x(r, 0, p)]. This ambiguity in the use of the 
functional symbol g ought not to cause confusion. 


“ Note that the coefficient of H(*) in (4.17) cannot vanish because of (4.4). 
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for the components of stress and strain associated with the displacement field 
u;(x) in Theorem 4.2 may be deduced from (4.11) by means of (4.1) and (4.2). 
Theorem 4.3. Let f;(x) be a body-force field simple in & Then 
(a) f;(7, 0, vy) admits the representation 


co 


17,9, ¢) = 2,9, 9), (4.18) 


=—0o 


where the {\ (v, 3, y) are uniquely determined solid spherical harmonics of degree k, 


=e; =f) (ie = Opiate Te sel= 2 50m:)5 (4.19) 


and the infinite series (4.18), as well as the series resulting from any finite number 
of termwise differentiations of (4.18), are uniformly convergent in every closed sub- 
region of @; 


(b) the same convergence properties apply to the series 


a s i= 2G 
2, k+1+(1—20) (2k+ 3) 


u, (7, 8, 9) = i (7,3, 9), (4.20) 
provided the constants yu, o obey (4.4), and further the functions u; so defined satisfy 
(4.6) in 6; 


(c) af n ts a fixed integer, 
17,0, ~) =0(r") implies u,(r, 8, p) =O (r"*}). 


Part (a) of this theorem is readily inferred from the simplicity of /;(x) in @, 
Definition 4.1, (4.8), and Theorem 2.1. 

The uniform convergence of the infinite series (4.20), and of the corresponding 
series resulting from successive differentiations of (4.20), is apparent from the 
analogous convergence properties possessed by the series (4.18). Consequently 
u;(r, 0, y) has partial derivatives of all orders in & Substituting from (4.20) 
into (4.6), carrying out the required differentiations termwise, and using (2.10), 
(4.18), (4.19), one readily verifies that ,(v, 0, m) is indeed a particular solution 
of (4.6). This establishes part (b) of the theorem. 

To justify (c) we note from (a) of the present theorem and (c) of Theorem 2.1 
that /;(7, 8, y) =o0(r") requires f\*) (7, 8, p) =0 for k2=n. Bearing in mind (4.20) 
and appealing once more to (c) of Theorem 2.1, we see that (c) of the present 
theorem holds true, and the proof is complete. 


5. Behavior of elastostatic fields at infinity 

At this stage we may turn to our main objective and examine the behavior 
at infinity of elastostatic fields which conform to the fundamental field equations 
in a neighborhood of infinity. 

Theorem 5.1. Let u;(x) be of class C® in &. Suppose u;(x), e;;(*), 74; (x) 
in & satisfy (4.1), (4.2), (4.3), subject to (4.4), the body forces being simple in &. 
Then, if n ts a fixed integer, 

u;(x) = 0(7") (5-4) 
g8* 
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Se u(x) =O"), ej(2) =O), 


(5.2) 


By hypothesis and Theorem 4.1, «;(x) is biharmonic in 6. On applying 
Theorem 3.1 to u;(x), we draw from («), (8) in part (c) of this theorem that 
(5.1) implies the first of (5.2). The remaining equations (5.2) then follow im- 
mediately from (4.1), (4.2), (4.3), and (6) in part (c) of Theorem 3.1. This con- 
cludes the proof of the present theorem. 

Taking »=0 in Theorem 5.1, we find that if the displacements vanish unt- 
formly at infinity, they must vanish to the order O (r+), while the associated strains 
and stresses must decay to the order O(r-®), and the body forces are necessarily of 
the order O(r-*) at infinity. It follows, in particular, that the dilatation and the 
rotation are both O(r-) in the present circumstances — a conclusion which 
DurFin & Nott [5] reached by entirely different means. 


Theorem 5.2. Let u,(x) be of class C) in &. Suppose u,(x), e;;(%), T;;(%) 
in & satisfy (4.1), (4.2), (4.3), subject to (4.4), the body forces being simple in &, with 


jilx) =0(r). (5.3) 

wi r,;(2) =0(1) (5.4) 

as u,(2) =8(x) 4009), e;;(2) =00%), a 
rj(%) =O), fix) =O), 


provided u,(x) represents an arbitrary rigid displacement field and hence has the 
form 


,;(x) = A; +} € 57. @; X, (A;,@;... constant) . (5.6) 


In view of the assumed simplicity of f;(x) and because of (5.3), we infer 
from Theorem 4.3 the existence of a particular solution w;(x) to (4.6) which is 
of class C® in & and conforms to u;(x)=O(r4). Accordingly, by virtue of 
Theorem 5.1, the strains and stresses associated with this particular w;(x) in 
the sense of (4.1), (4.2), obey the second and third of (5.5), while /;(x) must 
conform to the last of (5.5). We are thus able to exhibit a particular solution 
0 (4.4), (4.2), (4.3) which has the required smoothness and possesses the orders 
of magnitude (5.5) at infinity. Therefore and since the system (4.1), (4.2), (4.3) 
is linear, it remains to be demonstrated merely that Theorem 5.2 is true when 
/;(*)=0 in @. 

Suppose now /;(x)=0 in &. Then, by hypothesis and Theorem 4.2, u, (x) 
admits the representation 


u; (x) =u,(x) + 4,(x), (5.7) 


fe} a . 
where u; and w, are defined by 


u; (7,9, p) = Ds (1,0) Blt 7? Ds EEO a), 
ei aa (5.8) 
u,(7,0, 9) = dh (7, 8, y) +> H (7, 8,9), 
k=—oo 
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and the solid harmonics h\"), H satisfy (4.12). Further, each of the infinite 
series in (5.8) is uniformly convergent in every closed subregion of & and may 
be differentiated termwise arbitrarily often without loss of convergence. 


Next, define functions ¢;, and é@,;, through 


(5.9) 


From (5.7), (5.9), and (4.1) we have 

€;j(*) = e;;(x) + @;;(%)- (5.10) 
The two infinite series entering the second of (5.8) each represent a function 
harmonic in &. The second of (5.9) and part (c) of Theorem 2.1 therefore entitle 


us to write 
“(x)=O(r>), é;,(%) =O(r-*). (5.41) 


Because of (4.2), the order of magnitude at infinity of the stresses 1; ;(x) must 
be the same as that of the strains e;;(«). Also, if é;;(%) vanishes identically, 
it follows that “;(x) must have the form (5.6) appropriate to an arbitrary rigid 
displacement field. Consequently, and by virtue of (5.11), (5.10), (5.7), (4.2), 
all we have left to show is that ¢;;(”) =0 in the present circumstances. 


To this end we first conclude from (4.12) that 
Ip ee (eNO aesly ae Py on) (5:12) 
and hence reach with the aid of (5.8), (5.9), 


6:7 DME 2 De. (5.13) 


Eliminating H") from (5.13) by recourse to (4.12), and at the same time observing 
on the basis of (2.9), (2.10), (2.3), (2.2) that 


4%, Wet®) — (k+ 1 eo eres iy =0, (5.14) 
we arrive at 
= (1—2o) ( a (R) 
PS, PE sera oe (5.15) 
Now, by (5.4) and (4.5), 
é;;(*) =0(1), (5.16) 
so that from (5.10) and the second of (5.11) 
é;;(%) =o(4). (5.17) 


On applying Theorem 2.1 to the harmonic function é;;(«) given by (5.15), and 
bearing in mind (2.9), we draw from (5.17) that 


=e ior er 16 (5.18) 
which, together with (4.12) and the second of (5.14), yields 
H®=0 for k=—2. (5.19) 
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Inserting (5.19) in the first of (5.8), and using the first of (5.9) and the second 
of (5.14), we obtain 


22,,= > [A + hi). (5.20) 


Recalling (2.9) again, we gather from (5.20) that é;;(X) is harmonic in &. Also, 
in view of part (c) of Theorem 2.1, (5.20) is consistent with (5.17) only if 


é;;(*) =0 in @, (5.24) 


and the proof of the present theorem is complete. 


6. Extension of basic theorems in elastostatics to exterior domains 


We now apply the results established in Section 5 to the generalization for 
infinite regions of fundamental theorems in the linear equilibrium theory of 
homogeneous and isotropic elastic solids. In this connection we assume that 
the region occupied by the solid, whether finite or infinite, is a regular region 
of space, 1.e. a closed region D+ B whose boundary B consists of a finite number 
of non-intersecting closed “‘regular surfaces’, the latter term being used in the 
sense of KELLOGG [8]. Accordingly, the boundary 6 of a regular region of space 
D+ B may exhibit corners and edges. Further, D need not be simply connected. 
If D is infinite, B remains bounded, so that D in this instance is necessarily an 
exterior domain}. We turn first to an extension of Betti’s reciprocal theorem. 

Theorem 6.1, Let D+ B be a regular region of space. Let u;(x), e;;(x), T;;(*), 
f,(x) and u;(x), e;;(x), Ti; (x), fi(x) be two sets of displacement, strain, stress, and 
body-force fields having the following properties: 

(a) u;(x) and u;(x) are of class C® in D+B; 

(b) each of the two sets of fields in D satisfies (4.1), (4.2), (4.3) for the same 
elastic constants yu and o, the latter obeying (4.4); 

(c) if D is infinite, f;(x) and f;(x) are simple in D, in the sense of Definition 4A, 
and either 


u;(x)=o(1), wi(x) =o(A) (6.1) 
or 
t,;(%) =0(1), Ti; (x) =o(1), 
Epa ‘ (6.2) 
fi(x)=o(7%), f(%) =o(r), 
and 
ftdA+ffdV=0, fidA+ffdvV=o, (6.3) 
B D B D 
i tiie x;tyadA +S ein x;f,dV =0, og 
, t 4 
J Exip X tp aA et Eiip x;f,aV =0, 
where 
b, Teg; t= ij 0,, (6.5) 


while n; are the components of the outward unit normal to B. 


12 Cf. Footnote No. 2. 
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Then 
{umd + Jima — fimdd + [hu,aV = J 4; Gj AV — St e,,dV. (6.6) 


47°44 
For the case in which D is finite, the proof of this theorem is well known ®. 
Suppose, therefore, D is infinite. Let R(o) be the region whose boundary is 
B+Q(o), where Q(g) is a sphere of radius 9, centered at the origin, and con- 
taining B wholly in its interior. On applying (6.6) to the finite regular region 
R(o) + B+ Q(o0), we have 


A ea A a apt dV, 


6.7 
epee u,dA cape Pe eale terete ne) 
B 2(Q) R(e) R(Q) 
Sf ti;e,,4V =f 74, 6;4V. (6.8) 
R(Q) R(o) 
Assuming first (6.1) to hold, we draw from Theorem 5.1 (with 7=0) that 
u;(%)=O(r), u(x) =O(r), (6.9) 
es (%) =O(r), — t,;(*) =O"), F(x) =O(r), 6.10) 
exj(*)=O(r), (x) =O), fi(x) =O). 
According to (6.9), (6.10), and (6.5), 
ft,udA>o0, ftuj,dA>o0 as 9>~, (6.11) 
2(Q) 2(Q) 


while the volume integrals over R(o) in (6.7), (6.8) are convergent in this limit 
and tend to the corresponding (improper) integrals over D. Hence (6.6) follow 
by passing to the limit as goo in (6.7), (6.8). 

Next, assume that (6.2), (6.3), (6.4) hold instead of (6.1). In connection with 
this alternative hypothesis, it should be noted that (6.3), (6.4) are a consequence 
of the stress equations of equilibrium (4.3) and of the symmetry of the stress 
tensor if D is bounded, but constitute an independent assumption in the present 
circumstances; conditions (6.3), (6.4) could, however, be replaced with the require- 
ment that, as ooo, 


ft,dA 0, ft,dA 0, 
2(Q) 2) (6.12) 
Mc Secrmety aby tit x;,t,dA 0. 


In any event, Theorem 5.2 implies that 
u,;(x) = %;(x) +4;(*), ui (%) =m; + H;(%), (6.13) 
in which “;(x) and “;(x) are rigid displacement fields and 
&;(x)=O(r), Hi; (x) =O(r); (6.14) 
further, (6.10) continue to hold true. 


13 See, for example, [2], art. 109. 
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Since #,(x), ¢,;(x), Tj (x), fil”) and (x), e:;(%), Ti;(x), f(x) have properties 
(a), (b), and conform to (6.1), it follows from what has already been shown that 


[aa $f f,@av flags +f faaV fo ejaV =f vijejdV. 6.5) 
B D B D D D 


On the other hand, bearing in mind that 


or 


M1; (x) i oe hp,  Wi(%) =A; + b5 7p ®; Xp, (6.16) 
where /,, ®; and 7, &) are constants, we conclude from (6.3), (6.4) that 
fidA [fav =o; f[ta,dA+f fa,av =o. (6.17) 
B D B D ' 


Equations (6.15), (6.17), because of (6.13), imply (6.6). This completes the proof 
of the theorem. Theorem 6.1 at once yields, as a special case, the following 
analogous extension of the principle of work and strain energy ™. 


Theorem 6.2. Let D+B, u,(x), €;;(x), 1;;(%), f;(*), as well as w and o, 
satisfy the same hypotheses as in Theorem 6.1. Then 


B D D 
where W(x) ts the strain-energy density given by 
1 
Wiser silane! ee C5 5 O55 C55 O55 ; (6.19) 


and t, are the components of the surface traction on B, given by the first of (6.5). 
The truth of this theorem is immediately confirmed by setting w;(x) =u; (x) 

in Theorem 6.1 and by invoking (4.2). Theorem 6.2 serves as a basis for the 

subsequent generalized version of the uniqueness theorem of elastostatics. 


Theorem 6.3. Let D+ B be a regular region of space. Let u;(x), e;;(x), ti; (%) 
and u; (x), e;;(x), ti;(x) be two sets of displacement, strain, and stress fields having 
the following properties: 

(a) w;(x) and u;'(x) are of class C® in D+ B; 

(b) each of the two sets of fields in D satisfies (4.1), (4.2), (4.3), for the same 
body-force field f;(x) and the same elastic constants u and o, the latter obeying (4.4); 
(c) u;(%) =u; (x) on B,, 

, , ” ” (6.20) 

t; (x) = ti ;n,; = 4,0; =t; (x) on Ba, 
where B, and B, are complementary subsets of B, n, denotes the components of the 
outward unit normal of B, while t; and t;' are the components of the respective 
surface tractions; 


(d) af D ts infinite, either 


u;(%)=c;+o0(1), 4; (x) =e; + 0(1) (c;... constant) , (6.21) 
or 


Tse) — Cp 0(1) ay (a C,;+0(1) (c;;... constant) , (6.22) 


1s See, torexample, [2], art. 27, 
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with c;;=C;;, and 


a td A rail t; dA, (6.23) 
af Eijp %;tpdA =i &:;>%,t; aA, (6.24) 


1.e., the two systems of-surface tractions on B are statically equivalent. 
Then 


U,(%) = 4; (x) +4; (x), ej (*%) =e;;(*),  ti;(x)=Ti,(x) m D+B, (6.25) 


provided u;(x) designates a rigid displacement field. 


We emphasize that Theorem 6.3, in contrast to the preceding two theorems, 
involves no explicit restrictions of the body-force field. A proof of this theorem 
is readily carried out along the lines of the traditional uniqueness argument 
appropriate to finite domains. Thus let 

U;(%) =u; (x) — 5 (x), — e54 (x) = e5j(*) — e4(2), (6.26) 

4 5(%) = ti; (%) — 14; (%). . 

Clearly, w;(x), e;;(*), t;;(*%) meet the hypotheses of Theorem 6.2 with f;(x)=0. 
Therefore, and by virtue of (6.20), (6.26), 


{Way —0), (6.27) 
D 


where W(x) is the strain-energy density associated with e;,;(x), and the integral 
in (6.27) is improper but convergent if D is infinite. The second of (6.25) is a 
consequence of (6.26), (6.27), and the fact that W, in accordance with (6.19) 
and (4.4), is a positive definite quadratic form in the e;;. The remaining equations 
in (6.25) then follow from (4.1), (4.2), and this completes the proof. Evidently, 
u;(x) in (6.25) vanishes identically if B, contains at least one subregion of B 
with non-vanishing area or — in the event that D is infinite — if the regularity 
conditions at infinity take the form of (6.21). 

When D is infinite, B, is empty, and (6.21) are required to hold, Theorem 6.3 
yields in particular the uniqueness theorem for the first boundary-value problem 
which DuFFIn & Nott [5]'* established by different means, except that their 
argument remains valid even if ~ and o merely meet 


Mae, Wa corso <a5, Lis _20, (6.28) 


rather than the more restrictive assumptions!” (4.4). A uniqueness theorem 
analogous to Theorem 6.3 is known to apply to two-dimensional elastostatics '8. 
We also draw attention to the fact that the extension of the present theorem 
to mixed-mixed boundary-value problems is entirely elementary. 

If D is infinite and B, is empty, 7.e., if the surface tractions are prescribed 
over the entire boundary (second boundary-value problem) of a medium occupying 


18 See, for example, [2], art. 27. 

16 See also [6]. ‘ 

17 Theorem 6.3, of course, continues to hold if ~>0 in (4.4) is replaced by w<0, 
but the inequalities on o in (4.4) are essential to the argument we have used. 

18 See MUSKHELISVILI [13], art. 40. 


4114 M. E. Gurtin & ELI STERNBERG: 


an infinite regular region of space, conditions (6.23), (6.24) in Theorem 6.3 may 
be omitted because they are implied by (6.20). On the other hand, if B, is not 
empty (first or mixed boundary-value problem), (6.23), (6.24) assume the status 
of independent additional hypotheses. Further, in this instance the two hypo- 
theses under consideration, which require the prescription of the resultant force 
and moment (about the origin) of the surface tractions on B, appear at first 
sight artificial since the tractions over at least a portion of the boundary are 
here not known beforehand. One is therefore led to ask whether (6.23), (6.24) 
are necessary for the truth of the uniqueness Theorem 6.3 when D is infinite 
and the regularity conditions at infinity are taken in the form (6.22). That this 
is indeed the case is clear from the next theorem. 

Theorem 6.4. Theorem 6.3 is false if either the hypothesis (6.23) or the hypo- 
thesis (6.24) ts omitted. 

For the purpose at hand, let D+ B be the region exterior to the unit sphere 
about the origin, assume that B= B, (B, empty) in Theorem 6.3, take c;;=0, 
and suppose that w;'(x), ¢/(x), t;;(*), and f;(x) vanish identically in D+ B. 
Conditions (6.23), (6.24) then become 


{aA —0, (6.29) 
B 
J e:3p %jt4pdA =0. (6.30) 
B 


To prove Theorem 6.4, it clearly suffices to show that if either (6.29) or (6.30) 

are violated, there exist fields w;(x), e;;(x), t;;(x), which conform to hypotheses 

(a) and (b) of Theorem 6.3 with f;(~)=0, and are such that the strains e;;(x) 

do not vanish identically in D+ B, while 

u,(x)=0 on B(x;x,=1), (6.34) 
Ti; (*) =0(1). (6.32) 


We now demonstrate with the aid of two examples! that such fields do in fact 
exist. 


To this end consider the displacement field 


ae My y2—4 (3) 4 (A; = 0), (6.33) 


Y 2(5—60) 
in which the A; are constants, and the displacement field 


ES FO BN, vy2—4 Ej pq Wp % 

hammer murray war bed eee ame CO 
in which the w; are constants. Both of the foregoing displacement fields, together 
with their associated fields of strain and stress, possess all of the required pro- 
perties, as is readily confirmed. The leading two terms in (6.33) represent the 
displacements induced in the medium by a rigid translation (with components 4,) 
applied to the spherical boundary B, if the body is constrained against displace- 
ments at infinity; the last term in (6.33) corresponds to a rigid translation (with 


1° These examples were suggested by Professor R. T. SHIELD. 
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_ components —4;) of the entire medium. On the other hand, the leading two 


terms in (6.34) represent the displacements generated in the medium by an 
infinitesimal rotation (with components w,) applied to B, if the body is held 


fixed at infinity, while the last term in (6.34) corresponds to an infinitesimal 
rotation (with components —q@,) of the entire medium. Furthermore, if L; and 


M;; respectively denoté the components of the resultant force and of the resultant 
moment about the origin of the surface tractions on the unit sphere B, we find 
from (6.33) that 


pak (07 Goa 
whereas the displacement field (6.34) gives rise to 


Le 0.) VE = 870. (6.36) 


_ Thus, the solution of the field equations characterized by (6.33) violates (6.29), 


whereas that characterized by (6.34) fails to meet (6.30). The proof of Theorem 6.4 
is now complete. 


Our next objective is a generalization to infinite domains of the principle 
of minimum potential energy. For this purpose it is helpful to introduce the 
subsequent definition. 

Definition 6.1. Let u;(x) be a displacement field with the following properties: 


(a) u,;(x) is of class C® in a regular region of space D+ B; 
(b) i; =uj(x) on By, (6.37) 


where B, is a subset of B and ux (x) is a triplet of functions prescribed on By; 
(c) tf D ts infinite, 
7 C2 oe OO i 2 C22 el Oa A (6.38) 


where @;;(x) are the strains associated with %;(x) im the sense of (4.1). 

The totality of all such %;(x) constitutes the class of kinematically admissible 
displacement fields [%] corresponding to D+ B, B,, and uj (x). 

Theorem 6.5. Let D+ B be a regular region of space. Let u;(x), e;;(%), T;;(%), 
};(x) be a set of displacement, strain, stress, and body-force fields having the following 
properties : 

(a) u;(x) is of class C’) in D+ B; 

(b) ;(%), €;;(%), Ti; (x), and f;(x) im D satisfy (4.1), (4.2), (4.3) for elastic 
constants obeying (4.4); 

(c) U;=uy(x) on By, 470; =tf (x) on Ba, (6.39) 


where B,, By, and n; have the same meaning as in Theorem 6.3, while u(x) and 
t*(x) are functions prescribed on B, and By, respectively; 


(d) if D is infinite, f;(*) is simple in D, in the sense of Definition 4.1, and 


u,; (x) =o(1). (6.40) 
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Let [a] be the class of kinematically admissible displacement fields corresponding 
to D-+-B, By, and u¥(x). Let O{ii} be the functional of u;(x) defined by 


Piii} = J WE) dV —ff,u,dV —Jifu,aA, (6.41) 
D D By 


where W(e) is the strain-energy density given by (6.19) and @;;(x) are the strains 
associated with ti;(x) in the sense of (4.1). 


Then 
Pfu} = min O{ii} , (6.42) 


and this absolute minimum is assumed by D{ii} only if @;;(x) =e;;(x) m D+B. 
In outlining a proof of this theorem, we may confine our attention to the case 
of an infinite2° D. Then @f{i} is given by 


Pi} = lim ®, fi, (6.43) 
o> 0 = 
with 
®, {ii} =f We) dV — J f,a,dV — Jt a,dA, (6.44) 
EQ, R(e) By 


provided R(g) is once again the region whose boundary consists of B and the 
sphere Q(0) introduced in connection with the proof of Theorem 6.1. In view 
of (6.38), (6.40), and Theorem 5.1, the existence of ®{i#} is assured. By the 
present hypotheses and Theorem 5.1, u;(x) is a member of [#]. 


Now let 


u,(%) — u,(%), 6; (x) = €;;(x) — €;;(%)- (6.45) 


s 
S 
| 


From (6.44), (6.45), with the aid of (4.1), (4.2), (4.3), (6.37), (6.39), and the 
divergence theorem, follows the identity 


®, {ii} — O, {u} =f We')dV +f r,,n, u:dA, (6.46) 
R(Q) QQ) 


in which n, designates the components of the outward unit normal of Q(o). 
But according to (6.45), (6.38), (6.40), and Theorem 5.1, 


u;(x) =O(r),  e;,(x) =O (r=), T;;(%) =O(r?). (6.47) 
Hence, passing to the limit as @->oo in (6.46), we have 
Diu — Diu} se W (e') aV. (6.48) 


Since W(e) is positive definite, we draw from (6.48), (6.45) that 
Di} = Olu} (6.49) 


and that the equality prevails in (6.49) if and only if é;;(*)=e,;,;(x) in D+B, 
which is the desired result. A related minimum principle — also valid for both 
bounded and exterior regions, but applicable only to the first boundary-value 
problem (B, empty) — was established in [6]. We turn next to an extension of 
the principle of minimum complementary energy. 

20 A proof valid for bounded regions may be found, for example, in Hitt [14] 
p. 60 et seg., except that in [1/4] the body forces are assumed to vanish. 
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Definition 6.2. Let ;;(x) and f;(x) be a field of stress and a body-force field 
with the following properties: 


(a) 7;;(x) ts of class C in a regular region of space D+ B; 
(b) 7; ; (x) and f;(x) satisfy (4.3) in D; 
(c) ; t,,1% =tF(x) on By, (6.50) 


where B, 1s a subset of B, n,; denotes the components of the outward unit normal 
of B, while t¥ (x) is a triplet of functions prescribed on By; 


(d) if D ts infinite, 


qv 


ai) 0%), (6. 
ftidA=L, fe x,t,dA=M, (6.52) 


wD Oi G) 
in P 


where i= T;,;0;, while L; and M, are prescribed constants. 

The totality of all such 7; ;(x) constitutes the class of statically admissible stress 
fields [t| corresponding to D+ B, f,(x), By, t*(x), and (if D ts infinite) to L; 
and M,. 

Theorem 6.6. Let D+ B, u,(x), €;;(%), t;;(x), f;(*), u(x), U(x), as well as 
the elastic constants, satisfy the same hypotheses as in Theorem 6.5, except replace 


NN aa te Sow SOP): (6.53) 
iA = Eee lier, x71,44 = M1 ,, (6.54) 


where t;=t;;n;, while L; and M; are prescribed constants. 


Let Ae be the class of statically admissible stress fields corresponding to D+ B, 
};(x), Ba, t*, and (if D is infinite) to L; and M,. Let P{i} be the functional of 
T;;(x) defined by : 

Pir} = f We) dV — f t,uk dA, (6.55) 
D B, 


where W(e) is the strain-energy density given by (6.19) and eé;;(x) are the strains 
associated with 7; ;(x) in the sense of (4.5). 
het Ga t= min FAT. (6.56) 
[v] 
and this absolute minimum is assumed by V{t} only if T;;(x)=1;;(*) m D+ B. 


It will again be sufficient to sketch a proof on the assumption that D is 
infinite?!, In this instance Y{7} is given by 


WT} = lim YT} (6.57) 
“iia = = W(é) dV — g t,u* dA, (6.58) 


provided R(o) has the same meaning as in ‘fe proof of Theorem 6.1. It is clear 
from (6.51), (6.19), and (4.5) that {7} exists. Also, by hypothesis and Theorem 5.2, 


u,(*) =i; (x) +%,;(x), (6.59) 


21 See Footnote No. 19. 
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where %,(x) is a rigid displacement field, and 
ti;(x) =O(r), e;;(*) =O"), taj (%) =O). (6.60) 


In view of the present hypotheses and the last of (6.60), it is apparent that 
t;;() belongs to [7]. 


Setting ; 
T45(%) =%;;(%) — T4;(), 64; (%) = €;;(%) — e;;(%), (6.61) 
we draw from (6.58), (6.19), and (4.5) that 
Wi —W(p=fwe)dV +frjej;dV —Stup dA. (6.62) 
; K(g) R(e) B, 


Now, from the divergence theorem, (6.59), (4.1), and the fact that t;;(x) evidently 
meets (4.3) with /;(x)=0, 


Suipej,dV =f titdA +f ta,dA, (6.63) 
R(e) B Qe) 


provided B+QQ(o) is the boundary of R(o). Further, because of (5.6), (6.52), 
(6.54), and (6.61), 


Stu;dA=0, (6.64) 
B 
so that (6.63), in view of (6.50), (6.39), the first of (6.61), and (6.59), implies 
Stj;e,;4V =f tuk dA +f ta,dA. (6.65) 
R(Q) By 2Q(Q) 
By (6.62), and (6.65), 
Wit} —iink=fWwe)dV +f t,a,dA (6.66) 
R(o) 2(Q) 


or, on passing to the limit as @—>oo, and using (6.61), (6.60), (6.51), 
Pit} — Pfr} = f We’) dV. (6.67) 
D 


The desired conclusion is immediate from (6.67) by virtue of (4.5), the second 
of (6.61), and the positive definiteness of W(e). 


Both Theorem 6.5 and Theorem 6.6 are easily generalized to accommodate 
also mixed-mixed boundary conditions. It is also worth mentioning that the 
uniqueness Theorem 6.3 is a corollary of the foregoing two minimum principles. 

The results which have been established in this paper yield at once analogous 
extensions to exterior domains of the integral-representation theorems of classical 
elastostatics** and of Saint-Venant’s principle in the formulation given in [16]. 
On the other hand, the generalization of the present results to anisotropic media 


or to infinite regions whose boundary extends to infinity, does not appear to 
be elementary. 


The results communicated in this paper were obtained in the course of an in- 
vestigation conducted under Contract Nonr-562(25) of Brown University with the 
Office of Naval Research in Washington, D. C. 


22 See Love [1], art. 169. More explicit statements and proofs of these theorems 
can be found in [15]. 
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On the Scattering of Water Waves by a Circular Disk 


R. C. MACCAMY 


Communicated by C. TRUESDELL 


I. Introduction 


The boundary-value problem describing the diffraction of two-dimensional 
water waves by a dock of finite width has not as yet been solved in explicit 
fashion. It has, however, received considerable attention. The general problem 
of diffraction by surface obstacles was investigated thoroughly by JouN [/] and 
later by Peters & STOKER [2]. An existence proof for the finite dock problem, 
using the calculus of variations, was given by RuBIN [3]. Finally it was shown 
by the author [4] that the integral equation approach of [/] can be greatly 
simplified for the dock problem. 

It is the purpose of this paper to point out that the results of [3] and [4} 
can be extended to the three-dimensional problem of diffraction by a circular 
dock. Indeed it proves possible to reduce the three-dimensional problem to 
essentially the same (two-dimensional) problems studied in [3] and [4]. The 
technique is an extension of that used by HeErns and the author [4] to study 
the diffraction of sound waves by a circular disk. Briefly the solution is expanded 
in a Fourier series in the angular variable, and then the two-dimensional boundary 
problems for the Fourier coefficients are investigated. These coefficients satisfy 
rather complicated equations. However it turns out they can be expressed very 
simply in terms of certain harmonic functions of two variables, the latter being 
closely related to the values of the coefficients on the axis of the dock. When 
the boundary conditions on the coefficients are translated into conditions for 
the harmonic functions, one obtains essentially the problems of [3] and [4]. 

In addition to the fact that our methods yield an existence theorem together 
with a reasonably practical method for numerical solution of the disk problem, 
these methods seem to have an additional interest. It is the feeling of the author 
that many boundary-value problems for elliptic equations including as here 
those involving circular regions in three dimensions differ only in technical 
details from some analogous problems for harmonic functions of two variables. 
The connection is established by some integral representation of solutions of 
the equations in question. We use such a representation in Section V. Since 
a great deal is known about harmonic functions of two variables, one can hope 
to exploit this connection to obtain information about other equations. This 
paper is intended as an example of how one establishes the connection, and 
we illustrate its exploitation in the last section by finding the nature of the 
singularity produced at the edge of the dock. 
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II. General Formulation 
An inviscid, incompressible fluid occupies the region y<0, y=O being a 


_ free surface. We consider two problems. First suppose the motion is initially 


that of a plane progressive wave moving in the positive x-direction. A rigid 


circular disk is placed in the free surface and held in position, thus producing 
_ in addition to the original wave a diffracted wave pattern. All motions are 


assumed time-periodic, irrotational and of small amplitude. Suppose in the 


, second case the fluid is initially at rest, and then motion is produced by moving 


the circular disk vertically with forced periodic vibrations. By combining the 
two problems one could determine the oscillations of a freely floating disk in 
an incoming wave. 


If we divide all length dimensions by the radius a of the disk, the mathe- 
matical model of the above problems is as follows (see [7] and [2] for details) : 


Problem I (P. 1). Find a function u(r, , y) continuous in yO such that 
(A) 4,7 U, +7? Ugo4+%,=0 in y<O, 
(B) u(r,3,y) ~ AO) r-texp[ky+ikr] as r—+>oo, y bounded, 
©) 4, \7,9,.0)—2(7,9)— for 0=r<1, 
(D) 4, (7,9,0) —Ru(r,8,0)=0 for r<1. 


— 


Here x and z are co-ordinates in the free surface and 7?= x?+ 22, # = arctan 2/x. 
The parameter # is o%a/y where o is the frequency of the motion and y the 
acceleration of gravity. gis to be an analytic function of x and z on x?+4 271. 
For the diffraction problem we have g=exp(i7cos#@), while for the forced motion 
problem g=1. 

We shall impose on the solution a regularity condition at the edge of the 
disk*. If c, denotes the half-tube, (7—1)?+ y?=e, y<0, we require 


(II.1) lim ff uu,ds=0, | m=normal to c,. 
e—0 Ce 


With this condition it is shown in [2] that the solution of (P. I) is unique. 
We present now an integral formulation of (P.I). In a standard way we 
introduce a Green’s function, G(r, 0, 3, y, y), by 


G(r, 0,9, ¥, ¥) 
(11.2) =R2—kf (2 +n°)-bexp[hly —n)] dn +iake? HY (bu) =RALG, 
i 
R?=72+ 0?—2recos(@—y)+y?, w=r?+o2— 2recos(i—y), 


H) (kw) denoting the Hankel function. It is easily verified that G is harmonic 
for (r, 8, vy) = (@, y, 0). Further, G satisfies (B) for @ bounded and 


(11.3) Gy—kG= maeak 


x The author believes this condition can be shown to follow from conditions (C) 
and (D) and the continuity of wv in ySo. 
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For any function f satisfying a Holder condition on 0X1 we introduce 


the operators 


2a 1 
wrayvN= [ [Rev ededy, 
00 
ant 
LOA yN=z. ff ev Gl.d,» ov) ede dy, 
0 0 


Ut, 9,95 J) a sei 


Then for any such function f, U is a harmonic function in y<0, is continuous 
in y<0 and satisfies (B). Moreover from (II.3) we have 


é eo 
(11.4) ey k)U a W in y<0. 
It is immediate that U satisfies (D). On the other hand, by a formula of potential 
theory, (11.4) yields 


U, (r, 9,0; f) —kUG, 9,0; =f(,9), OSr<t. 


Accordingly U(r, 9%, y; f) will be a solution of (P.I) if f is a solution ot the 
integral equation 


(11.5) tir, 0 +k[W7,8,0;N+Lr,0,0;N\l=e,9) OSr<i. 


It can be shown that U(r, #, y; f) satisfies (11.1) for any 7. If the homogeneous 
equation corresponding to (II.5) had a solution f)(7, #), then clearly u(r, 3, y) = 
U(r, 8, y; fy) would be a solution of (P. I) with the right side of (C) replaced 
by zero. Hence, by the uniqueness theorem quoted, u(7, #, y) =0, and then 
by (11.4) 

fo(r, 8) = lim (5 — k) U(r, 8; fo) =0. 
Thus the homogeneous equation can have no non-trivial solutions. From (II.2) 
we see that the kernel of equation (II.5) has a singularity of the form w+-c log 1; 
hence the Fredholm theory yields a (unique) solution f(v, #) of equation (II.5). 
We summarize the results of this section. 


Theorem 1. The integral equation (11.5) has a unique Hélder-continuous 
solution f for any function g analytic on OSr<1, and the unique solution of 
problem (P.1) is given by 

u(r, 3, y) = U(r, 0, y; fy. 


III. Complementary Problems 
In this section we use the methods of [4] to transform (P.1I) into a simpler 


problem. We start with formula (II.4). Applying the operator —s-(s+4) 
to both sides yields dy \ by 


(III.1) (-jrt¥)2=- pat (3 +k) W in y<0. 
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a . . . 2 
Now both U and W are harmonic functions in y<0; hence we may replace uc 


oy? 
C2 Uae (a) 1 é ar 5 
b A= x : 
y apt or a obtaining 


oy oy 


Suppose now that U(r, #, y; f) is a solution of (P. 1). By (C) we see that the 
harmonic function 


(III.2) (A + 2) —=A(Z - k) W in y<0. 


0U 
V(r, 0, y) a Oy (7, 0, y; /) 
satisfies 
V(r, 3, 0) = g(r, 9), OSSeaak 


Since g is an analytic function, it follows that V and hence U itself may be 
continued across the disk into y>0 as a harmonic function. From (III.2) then 
it follows that the harmonic function 


ACs +k) (7,8, 951) 


and hence W itself can be continued across the disk as a harmonic function. 
Thus U and W are analytic for yS0, r<1, and we may pass to the limit y=0 
in (III.2), obtaining 


(IIT.3) (A + F?) g(7,0) =A(W, + kW). 
Now let 41(h) denote a frxed solution of the two-dimensional Poisson equation 
Aw=h. 


Then (III.3) yields 
(11.4) (a +2) W = [1+ A ee a Aa: 

The operator W satisfies 

W, (7, 8,0; f) =0, yr>1, 

so that (III.4) suggests that we should consider the following complementary 
problem. 

Problem II (P.1I). Find a function v (7, 8, y) continuous in yO such that 
(A’) Ue T=, Det Te Up prt Uy =0 in yx O, 
(Be => Oyfasyar72= 7S co, 
(C’) v,(7,8,0) +kv(7, 8,0) =A(7, 0), 7<1, 
(2) (7,0; 0) =O. 1 1. 
If we further require the condition (II.1) for the function v, one readily finds, 


using Green’s theorem, that the solution of (P. II) is unique. Also we can 
prove the existence of a solution, as for (P. I), by finding a unique f(r, #) satis- 
fying 

(III.5) (7,9) +kW(r, 8, 0; f) =A(7,9), Pal 

and setting 


/ 


v(r, 3, y) = Wir, 8, y; f). 
o* 
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Note that 
(III.6) f(r, 9) =v,(r, 8, 0). 


The problem (P. II) is simpler than (P. I) because of the simplified kernel 
in the integral equation. It would be advantageous to be able to assert that 
the solution of (P. II) yields the solution of (P. I). However, so far we have 
shown only that the solution of (P. I) yields a solution of (P. II). Suppose now, 
as suggested by (III.4), we find a solution of (P.II) for h(r, #)=[1+hk?A4] g 
and set f(r, #)=v,(r, #, 0). Then by (III.2) we could deduce only that 


(A +k) U, (7, 9,0; f) =A + A) g(r, 8) = (4 + F) g(r, 9) 


from which (C) does not follow. 

The difficulty here is of course that the null space of the operator (4 + k?) 
is infinite-dimensional. To circumvent this problem we try to find approximate 
solutions in spaces in which the null space of 4-++A? has finite dimension, and 
this leads us to a Fourier series decomposition. 


Theorem 2. If the function g(r, 0) of - I) has the form, 
cre A= r= Sant )cosnd, 


there exists a uniquely determined function h(r, 0) such that if v ts the solution 
of (P. II), then the solution of (P. 1) ts given by 
u(r,9,y) =U, 8,951),  £(,9) =2, (7,8, 0). 
Let w, (7) be solutions of 

(III.7) wr tr w,—n?r2w, =a, (r) 
with w,,(0) bounded. We find 

r Q 

Wy (r) =r" fo?" do f a, (zt) dt. 
0 0 


Convergence of these integrals is assured by a lemma of [5] which shows that 
the analyticity of g(r, #) implies that the a, (t) must be O(rt") near r=0. Let 


h’ (ry, 9) denote the sums 
N 


h® (r, 8) = = 2 (an r) + kw, (r)) cosnd. 


Let bas y) be the solution of (P.II) for h(v,9)=AN(r, 8) and v er, Owe 
n=0,1,...,N be solutions of (P. II) for h(7, 9) =r" cosn #, respectively. Follow- 
ing our previous remarks, we have 
We.d.ysf, wv =W,8,93 7%), 
f+ Wir, 6,0; =I, 8), 
fh + Wr, 9,0; f) =r" cosn®. 


By straightforward calculation we may verify that the functions f and f” must 
have the form 


aly 
f=21G,(r)cosnd, jf" =F (r)cosn@. 
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We show that the (unique) solution of (P. I) with g(r, #) equal to G(r, 8) 
can be written in the form 


~ N 
u(r, 8, y) = UY,8, 9; f) + 24,U(7;9, 9; 
n=0 


for some appropriate constants A,. From (III.2) we have for such a u, 
3 eH 
(A + #4) u, (7, 8,0) =A(5- +k) (We, 9,0; +2 A, Wr, 0,0; Pt 
y n=0 
N 
= AL (7,8) + YA, cosn dh. 
n=0 
But A(r” cos n #) =0, and by (III.7) 


A tO) ) = A (a 7) + kw Shen tA aes wh cas 
Hence a 
(A + #4) u, (7,8, 0) =(4 + 8) G(r, 9), 


and accordingly 
(III.8) u, (7, 3, 0) = G(r, 8) + K(r, 9), 
where K(r, #) is a solution of 
(III.9) (A +k?) K= 
Again it is not difficult to verify that 


U(r, 8, 0; h(e) cosny) = B(r) cosnd. 


Fe AN 

It follows that the tunction U(r, 8,0; f)+)>) A, U(r, 8, 0; 7”) has only N non- 
n=0 

zero Fourier coefficients. Since this fact is also true of G(r, %), we see by 

(III.8) that it must be true of the function K(y, #). But K(r, 8) is also a solution 


of (III.9), and hence we must have 


N 
=>) a, J,(kr) cosn®, 


n=0 
the J,,’s being regular Bessel functions. 
The course is now clear. We choose the constants A, so as to make the 


coefficients «, in K(r,#%) all zero. Since j” O=s and j{”(0)=0 for n<k, 
the conditions 
o” 


K=0 at P=0. TOs son, INE 
or" 


will insure that K(r, #)=0. Then according to (III.8) we choose the constants 
A,, as solutions of the eee 


P00.) 


” Orr — (0, 0), =O; sly veo Ne 


(III.10) & U(0, 0, 0; j+d4 


This is a ie eee aoe of linear equations for the constants 4A,,. 
It will always have a solution unless the associated homogeneous system should 
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have a non-trivial solution. Suppose the latter were the case, and let A® denote 
that solution. Then the above argument shows that 


N 
u(r, 0, y) => ALUG RYT) 
n=0 


is a solution of (P. I) with the right side of (C) replaced by zero. Hence by the 
uniqueness theorem w°(r, 0, y)=0. Recalling that 


U, (7, 0,0; f) —k U(r, 3,0; ) =f(7, 9), Osr<i, 


we infer then that 


0=W (7, 3,0) —kw(r, 9,0) = > Ag? 


But by construction we have 
N N \ W 
¥ Ao php Wr, 8,0; YAP) =D ASP cosnd. 
n=0 n=0 n=0 
Combining the last two formulas, we conclude 


N 

» 49" cosn = 0 

n=0 
from which A? =0 ~=0,1,..., N. The contradiction shows that the system 
(III.10) always has a solution. The proof of Theorem 2 is now complete. 

We make some further remarks concerning the coefficients in (III.10) in 
the Appendix, but here we have at least in principle reduced the solution of 
(P. I) with right side of (C) equal to GY (r, #) to the solution of (P. II). Trans- 
lated into integral equations this means that we can express the solution of 
(II.5), with right side equal to GY (7, 9), in terms of solutions of (III.5). 

Let {* (7, ®) be the solution of (II.5) with right side equal to GY (r, 9). Since 
the functions G(r, #) converge uniformly to g(r, %), one may expect that the 
functions f% (7, 9) converge to a function f(r, #) which is a solution of (II.5) 
itself. We verify now that such is indeed the case. 


Theorem 3. Let GN(r, 8) be a sequence of functions converging uniformly to 
g(r, 0) in OSrS1, and suppose the functions f(r, 9) are solutions of the integral 
equations 

fX (7,8) + RW, 8,0; P) +.L0,8,0;P)]=@*0,9), o<rei. 
Then the functions f® (r, 8) converge to a solution f(r, 8) of the equation 
f(r, 3) + kW, 3, 0; f) + Lr, 0, 0; J] =e(r, 9), OS 7k 
We need first a lemma expressing the complete continuity of the operator U. 


Lemma. Let /"(r, 3) be a uniformly bounded sequence of functions continuous 
and satisfying a Holder condition on OXr<1. Then the functions U(r, 8, 0; f”) 
form an equi-continuous family. 

For ease of writing we denote the integration points (0, p,0) by Q, the 
variable points (7, 8, 0) by P, and the distance from P to Q by PQ. Now by 
(II.2) 


" 1 
C= + KP, 0) 
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where K is square integrable. Since L was defined on p. 122 as an integral 
_ operator with square integrable kernel, it follows easily by Schwarz’ inequality 
and the boundedness of /” that the difference 


| LBP) — LPs?) 
can be made uniformly small by choosing P, P, sufficiently small. The diffi- 
in G. If |f’| <M, we have 


=| [fre gon RQ) @% 


ak 


culties arise from the term 


1 
PQ 


2x |W(R; f") — WR; f)| 


SLIT.44.) 
dSo. 


alle 
EQ #2 
Now let d=, P,, and let s, denote the intersection of a disk of radius 26 and 
center P, with the disk r<1. We have 


wi vmilieell: 


Sy (y<1)— 


By a known result of ele theory ([7], p. 129), for any region s of the plane, 


of area A, we have 
[rs <2(nA)! 


independently of P€s. Thus we find 


(11.12) sy | 450 S 4 (402 02)! = 8x0. 


Suppose on the other hand that Q€[(r<1)—s,]. Then we have 


LEO Olin Ours bo 0 aah 0ien 0, 

P,Q—-BO dso 
dso= {fron [t= If pod» 
Let us choose new polar co-ordinates (9, ) with origin at R,AQ=e. Then 
our last inequality becomes 

22 R( oF a 2a1 ee 5 

Cae QaK {2 
ey lenges <6] | g—o = 2216|log 20 ) 
(r<1)— 0 26 


If we enter “the estimates (I11.12) and (III.13) in the inequality (III.11), we 
obtain 


and 


ROBO 


(r<1)— 


22 |W(P;; f") — W(P;; f")| <M0(8x + 2m log*=**), 


The right side becomes uniformly small with 6, so that the family W(P; 7”) 
is equi-continuous, and the lemma is proved. 
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We return now to the equation 
(111.14) f?(P) + U(P; f’) =G"(P), 
in obvious notation. We wish to prove that the sequence /” is uniformly bounded 


and equi-continuous. Suppose first that the f” are uniformly bounded. Given 
any e> 0 we can find by the lemma a 0, so small that PP,< 6, implies 


|U(P; f) — UB Pyl<ze 


uniformly as to ”. From the uniform convergence of the GN (P) to g(r, 8) it 
is clear that we can also find a 6, such that for PR.< od, 


|G(P) —G(P)| <4. 
Then if 6=min(6,, 6), we have 
| (P)—f(B)|<e for PA<6 uniformly inn. 


It follows that the 7, are equi-continuous; hence a subsequence converges 
uniformly to f(r, %), a continuous function. It is well known (and proved 
essentially as the lemma) that an integral operator of the form of U is continuous *. 
Hence we have 


g(r, 9) = lim GY (7,8) = lim (7% + U(r, 8, 0; f%)) =f+ U(r, 8,0; f). 


We must still consider the possibility that the /” are not uniformly bounded. 
Let o” be the maximum of /"(7, 8) on r<1 and g”(7, #)=f"/o”. Then 


(111.15) ge" (P) + U(P, g") =G"Jo". 


If the o” are not bounded, there exists a subsequence m; tending to infinity 
so that G”/o” tends to zero. We could then repeat the steps of the preceding 
case to infer that there is a further subsequence {g”*} converging to a continuous 
function g which by (III.15) must clearly be identically zero; that is, 

lim g”*(r, 3?) =0. 


k—>0oo 


But max g"*(7, #)=1 for all ~; hence we have a contradiction, showing the /” 
must be uniformly bounded. Theorem 3 is thus proved. 

The results of this section show that at least up to a convergent approxi- 
mation method we can reduce the solution of (P.I) to that of (P. II). In 
Section IV we point out that the latter can be formulated as a variational 


problem for a positive definite functional. In Section V we present an alternative 
method of attack for (P. II). 


IV. A Variational Problem 
If one attempts to formulate a variational problem equivalent to (P50; 
two difficulties are encountered. First the functional involved is not positive, 
and second the conditions at infinity are such that this functional is not bounded. 
These facts are to be expected since the problem, like that of diffraction of 
acoustic waves, involves infinite energy. The same difficulties were observed 
by Rupin [3] in two dimensions. He was able to overcome them by considering 


* Continuous with norm taken as the maximum of the function on 7X1. 
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a complementary problem related to our (P. II). We show in this section that 
the same technique applies in the three-dimensional case. 

Let F denote the class of functions v(7, 9, y) which satisfy the following 
conditions: 


1) v continuous in y<0, 


fl 
rc 


2) v of class C® in y<0 for (7, y)= (A, 0), 

3) v0 as 72+ y2>0, 

4) v satisfies condition (II.1), 

5) MEN a3 (grad v)2dt+ ff (ku—h)?do<~x. 
y< fal 


y=0 
The volume integral in (5) is the limit of an integral over a finite volume. 
Letting C, denote the half-tube (y—1)?+ y2=e, ySXO of condition (II.1), 
the hemisphere 77+ y?=R?, y<0, and D, z the region bounded by C,, J; and 
the plane y=0, we can define the volume integral by 


SSF ( (grad v)?dz = lim Sff (grad v)2 dr. 


y<0 0D 
pues pe 


We now pose a variational problem. 
Variational problem. Among all v€F find that one minimizing F[v]. 


Theorem 4. The variational problem and (P.II) possess the same uniquely 
determined solution v. 

Suppose that v is a solution of the variational problem. Let 7 be any 
function of F which vanishes identically outside the region D, p above. Then 
the vanishing of the first variation of F at v implies, by Green’s theorem, 


1—e 2x 


— Sff (4 V22) arr f an (v, st Bay Are ea eens ar 0. 


De,R 1+e0 

Since 7, ¢ and R are arbitrary, it follows that v satisfies conditions (A’), (C’) 
and (D’) of (P. II). The continuity of v in y<0 and conditions (B’) and (II.1) 
follow from the fact that v€ F. This shows that every solution of the variational 
problem is also a solution of (P. II). Since we know, however, that the solution 
of (P. II) is unique, it follows that there can be at most one solution of the 
variational problem. 

We have shown in Section III that (P. II) does possess a solution v which 
satisfies conditions (1)—(4) of the variational problem. Since it was found in 
Section III that the solution v could be written in the form 


v(r, 0, vy) = Wr, 8, y; f) 
for some suitable function f, condition (3) may be sharpened to 
(IV.1) v=O((7?+y%)-#), gradv=O((?+y%)"#) as 7+ y2> 00. 


Accordingly we can show that F[v] for this (unique) solution of (P. IT) is finite. 
From Green’s theorem and the fact that v is harmonic in y<0 we have 


(yah (gradv)tdr= [fv do, 


Deg, R oDe, R 
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where 0D, p denotes the boundary of D, x. Using conditions (C’) and (D’) of 
(P. II) thus yields 


1—e 22 


[iliseerie tae tile Bee pete] | 0 Fa ziaee 


De, Rk 
Since v is continuous in yO, the ie integral has a limit as e—0. The first 
and second integrals approach 0 as ¢-+0 and R->oo by conditions (11.1) and 
(IV.1) respectively. Hence F[v] is finite. Thus the solution v of (P. II) is an 
admissible function for the variational problem. We show that it actually yields 
the minimum for F[v]. Let 7 be any other admissible function; then (with 
the same limit interpretation of the volume integrals as before) 


Flv +] =F[v] +8 [fda +k II (grad y)2d7 + 


(IV.2) “i. Bh aes ha ee atame (kv —h)do. 


y<0 ee I 
y=0 


From Green’s theorem and conditions (II.1) and (IV.1) we have 


2k fff (gradv-grady) dt = 2k ff nuv,do. 


y<0 (AG 
y=0 
Thus by (C’) of (P. II) equation (IV.2) yields 
Flv +4) = Fle), 


equality holding only if 7=0. This shows that v is a solution of the variational 
problem and completes the proof of Theorem 4. 


V. Solution of the Complementary Problem 
If we study the solutions of (P. II) necessary to solve (P. I), we see that it 
suffices to consider (P. II) for 
(V.14) h(r, 0) =r" h,(r\ cosn®, 


the functions h,,(r) being analytic functions of 72. We show in this section that 
this problem may be further reduced to a boundary-value problem for harmonic 
functions in two variables. We state the results here. 


Problem III (P. III). Find a function o(y, e) continuous in y<O such that 
yi: Poo = 0. nian 0, 
(BY) w+>0 as 92+ v2 00, 


(C") Py (0, @) + # P(0, e) = m(o) 
—1i<e<+1, m/(o) analyticin —1S 95 +1, 


(D") (0,e)=0 |o|>4. 


Theorem 5. Problem (P. III) possesses a unique solution 


P(¥, Q) = Ply, @; m) 
for each m/(o). 
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Theorem 6. Let u(r, 3, y) be the solution of (P. 11) for 
h(r, 0) =r"h,(r)cosn8, —h,(r) an analytic function of 7?. 


Then there exists a unique function m(o), analytic in —1S0S+1, such that 


Y 


(V.2) v(r, 0, y) =2c,7-” Sf (72 — 0?)"-3 Dy, 0; m) de|cosn#, 
0 


the constants c,, being independent of m(o). 


| We prove Theorem 6 first. That a solution of (P. II) exists follows from 


‘Section III. Also as remarked in that section the solution, for the choice (V.1) 
of h, must have the form 
v(r, 8, vy) = Wir, y) cosn®. 


Since, however, v(r, #, y) is harmonic and hence analytic as a function of x, y, 2, 
it follows from a lemma in [5] quoted above that 


Wir, vy) =r wir, y). 
From (A’) and (C’) follows 


(V.3) D201) rw) -— w,,— 0 iney <0, 
(V.4) Wy (7,0) + kwi(r, 0) =h, (7), Or An 
Let us recall that v(z, 3, y) can be represented in the form 
v(7,0,y)=W,d, 951), (7,0) =v,(7,0,0). 
Accordingly we have 


i= ee ee 


out lod 
w(r, y) cos 8 = cosn 9 [fe Ze AALS 
oer —2rocos+y23 


From this formula follows the relation 


1 
é F(g) "t+ de 
V. n! w (0, y) = sal fey <0. 
A proof of a more general result including (V.5) is given in [4]. 
So far all quantities introduced have been real. However w(r, y), being 
y~"v, is analytic for real y<0. The formula (V.5) serves a definition of w(0, y) 


for complex y=é+77, <0. Furthermore we can use (V.5) to continue w (0, y) 
é é 


as a function of complex y. For y=in, |n|>1 we have B63) BaP? and 
from (V.5) 
1 
[72 (er ees ffs Cre 
[o2-+y2? ig ane 

0 
Hence, 
(V.6) Rew(0,i7)=0 on |y| >1. 


It follows that w(0, y) n| > into E>0 so as 
to become an analytic function on the complex y-plane slit from —7 to +7 
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along the imaginary axis. Since w (0, y) is real for y real, we have also 


(V.7) w (0, y) =w(0, 9). 
We write 


w(0,€ +19) = (En) ++ (Em). 
p(é, 7) is then harmonic in §<0, and from (V.5) and (V.6) we have 
(V.8) y(0,y) =0 on |n|>4, 
(V.9) g(§ —n=eEn), ve —n =—v¥lén)- 
Finally from (V.5) we observe that 
(V.10) plEn) +0 as E4720. 


We have shown now that (&,7) satisfies conditions (A’’), (B’’) and (D”) 
of (P. III). In order to obtain condition (C’’) and to work toward the specific 
formula (V.2) we make use of the following fact concerning equation (V.3). 
It is well known (see for example [7]) that solutions of (V.3) are uniquely de- 
termined wherever they exist by their values on any portion of the axis of sym- 
metry, y=0, provided those values be analytic. The answer can be written 
down explicitly in the “Laplace’’ integral 

+9 


w(r,y)=c,7->* f 7—o%)"4Fw(0,y+to)de, ¢,=I(n)/l(n+3)I). 


This formula is valid in our case for y< 0 since we have observed that w (0, € +7) 
is analytic for €<0; by (V.9), it may be rewritten as 


te 


(V.11) w(r,y) =c,7~ >” f (x? — 9?)"~4 p(y, 0) do. 


0 
Substituting in condition (V.4), we find 


Y 


"Vy (1) = 2c, J (7° — 9°)" * [oy (0,0) tke e)]de, OSr<t. 


Setting 
r—=a, oe =B8, [9,(0,6*) +k —(0, B)] = x(p), 


we express the last equation in the form 


a hy el) = 2, J (a —By'4x(B) dB, O<a<1, 


or, differentiating ” times with TESPECtLONoe 
Pee ee eee ko. = 
(V.12)  A(u) = 5 (a" h, (a3) == ft —p)'x(8)dB, 0Sa<1. 
0 
Equation (V.12) is an Abel equation which can be inverted to yield 


ee zz f @ B)-*A(B) dB, OoX<a<1, 


0 
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for, setting D= = and returning to the original variables, 


f 
T 
d Q 
| Py (0, e) +k p(0, 0) =m(e) = sl ag | 8 ee D*(P*h,(0)) at, 
(V.13) “0 
| Oe. 


Since by (V.9) (y,@) is an even function of o, (V.13) must also hold on 
/—1<o0%. It is proved in Section VI that m/(g) is an analytic function of 0. 
The last member of (V.13) yields the desired function m(o) of Theorem 6, and 
equation (V.11) yields formula (V.2). Thus Theorem 6 has been proved. 
Having shown that (P. II) can be reduced to (P. III), we turn to Theorem 5, 
which asserts the existence and uniqueness of the solution of (P. III). We omit 
the rather simple proof of uniqueness, which, like that of (P. IJ), follows directly 
from Green’s theorem once account is taken of the singularities at the edges, 
(y, 0) =(0, +1). The necessary estimates near these edges are contained in 
formula (VI.5) of the next section. 
| In order to construct a solution we proceed as with problem (P. II). We 
seek a solution in the form 


aul 
(v.14) P(e) =—s5 5. [ wl loge — 0? + y%] dt 


for some function u(t), Hdlder-continuous on —1S#<-+1. This form insures 
that conditions (A”), (B’’) and (D”) are satisfied no matter what the choice 
of uw. By a familiar formula we have also 


y—>90 


41 
(0,0) =tim | FF fw ogile— "+ 98124 = u(t), me Ae Fee A, 
=i 


Assume now that w'(¢) exists and is Hélder-continuous on —1<#<-+1 and 
that 


(V.15) (+1) =u(— 1) =0. 
Integrating by parts, we have 


+1 


Sil 5 
sor J vl log[(g@ — #)? 4 y]}d=— [ ul) Feloglle —1)* + y*] ae 


f (t) (@—?) 

ENO gy 
fears ; 
=i 


Thus if the form (V.14) is to yield a solution of (P. III), w(¢) must be chosen 
so that 


(V.16) wt feOat=m@, -1<e<+1, 
=il 


the integral to be interpreted as the principal value. 


* Since h,,(t) is an analytic function of t?, D” (t?" h,(t)) is an even function of T. 
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Equation (V.16) is a special case of the Prandtl integro-differential equation 
of airfoil theory and has been extensively studied. It is shown in [8] that 
(V.16) can be reduced to a regular Fredholm equation of second kind, and 
from this it is established that there exists a solution yu (¢), with Hélder-continuous 
derivative, which satisfies (V.15). It follows then that (V.14), with w chosen 
as the solution of (V.16), is indeed a solution of (P. III). 


VI. Singularities at the Edge 

It is intended that some numerical computations based on the preceding 
considerations will be presented at a later time. Here we discuss a different 
problem. One feature of two-dimensional problems is the existence of a fairly 
extensive theory of the singularities produced at the confluence of two analytic 
boundary conditions, the theory having been initiated by Lewy [9] and con- 
tinued by LEHMAN [10]. Since we have reduced our problem to a two-dimensional 
one, it is natural to apply this theory and then translate back into original 
problem. This we do now following initially the pattern of [9]. 

Let us note at the outset that the function m(o) of Section V is analytic. 
For as remarked the /,,(7) are analytic functions of v?; hence in (V.13) we have 


t Dl has\a= 5) as 
m=0 
Since 
T 


Q 
i pil (gt i78) hg ot tr ie oe) see eee ee Ry <i meply 
0 
we have accordingly 


(V1.1) m(p):=> one 
m=0 
a convergent series for all a. 


From (V.7) we have, writing {(y)=w(0, y) for complex y=&+7y, 
i) =a SiG 
Thus we have 


292(5,9) =f(y) +79) =f) +F(—9), 
27 (E.) =f(y) +7) =H(y) —#(—9), 
and condition (V.13) becomes 
(V1.2) f(y) +7") + ALF) — f(y] =2G(y) for argy—=a/2, |y|<+2. 


Here ye denotes the point obtained from y by a counterclockwise circuit about 
the point y=7. If we set 


y 


(VI.3) Fy) = f(y) — f Ble) do, 
we can rewrite (VI.2) in the form ; 


(V1.4) EN(Y) AEP et aE) Ae 


again for arg y=27/2. 
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While we have derived (VI.4) initially for argy—a/2, we may use it as a 
_ functional relation to continue F(y) onto the portion |i— y| <2 of a Riemann 
surface with y=7 as a logarithmic branch point*. Observe that if F(y) is given 
on the sector —a/2<arg(y—1)<+3.2/2 of this surface, equation (VI.4) can 
_ be considered as an ordinary differential equation for F(y*) and 3 2/2< arg (y*—1) 
_<7a/2. By continuing in both directions F(y) is continued analytically to all 
values of arg(y—7) on the Riemann surface, the equation (VI.4) holding when- 
| ever |y—i| <2. 


It is proved in [//] that a solution of (VI.4) must have the form 


(VIS) = F(y)=c+> dian, (y —1)"*4[(y — 4) log(y—1)]", ¢ a constant, 
0 


m—0 2— 
the series converging absolutely and uniformly for |y—7| <p, @9 any positive 


number less than 2. Consider now the representation (V.11). Recalling that 
w(0,¢)=7(C), we have from (VI.4) and (VI.1) that 


TB OLN SAN pd ote ee 
m=1 
But Y m rv 
Te emer n) oer) SAAD: Oregon 0 Edo 


4 {3 
m 


= Pee ag yan—2 ke 
k=0 


Hence we can write (V.11) in the form 


+r 


(VI.6) Pwr, y) = 20,7" f (P?— oe)" 4Fly +10) do + S(r,y) 


=f 


where S(z, y) is a power series in 7 and y. Now if we introduce the variables 
CV U0, 2 = Vat, £=V— 27, 


the integral on the right side of (VI.6) can be written as 
sere 


f— ey Fy +10) do=S (ef CDI FO at =I), 


roams 


the integration being along any path joining z to z without crossing the slit, 
Ref€=0, —1<Im€<+1. With the correct choice of path we have 


i(e,2) =2Re f [(e-0 CDP IFO dt =2Re lL, , 7) 


along a path joining 0 and z without crossing the slit. Setting §=¢ —7, y=z—1, 
we have 


(V1.7) I; =f ln — 8) (§—H +20) FE +7) dé. 


Suppose now that || <1 and that we choose for the path of integration in J, a 
straight line joining =— 7 to =y. Then for & on this path we have |§| <|7—22 


, 


* In this relation y* of course continues to have the same meaning as above; 
that is, it is an image of y on the Riemann surface. 
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or, since || <1, 


(VI.8) €—g+i iS Diane 
If we enter the expressions (VI.8) and (VI.5) into (VI.7), we find 
=> SBS — ret ogey as 
=} > B, (7) iret f ( i) ili Ai 2s Se es Se) MC 
r=0 s=0 es 


the coefficients B,, being analytic functions of 7. From this it follows with a 
little calculation that Pa 


I, —H Sy Sie 7] ™ (n log n)’ 


k=0 m=0 7=0 


if |7| <1. Thus, if we collect our results, we find that 


w(r,9)=Red YY aus le — i" (2 +4)” [ze — 8) log (2 — a], 
(VI.9) k=0 m=—0 7—0 


z=y+r, 


the series converging for |z—7|<1. Note that the function S(r, y) of equation 
(VI.6) has been incorporated into the series (V1.9). 

Let us now recall the role played in the solution of (P. I) by solutions of 
(P. Il). We first found a finite number NV of solutions of (P. II) having the 
form v"(r, y)cosn@. Then for a suitable choice of constants A; the function 


eS thes vy (7,0) cosnd 
n=0 
was found to be equal to 


U, (7, 8,0) —k U(r, 8, 0), 


where U is a solution of (P. I) for the right side of (C) equal to G(r, 8). Since 
each function v”(r, y) is of the form 7"w (7, y), where w is a function of the type 
we have studied in this section, we conclude from (VI.9) that 

%7,0=2 2D Fm(i—7* 1 —7) log (—7)]™. 
Thus we have ee eD 


(VI.10) U, (7,8, 0)—k U(r, 8, 0) = > ys rm(1—7)*[(14—7) log (1 —71)]” cosn®. 
n=0 k=0 m= rs 
Finally we remark that in this solution U,(r, 8, 0) equals GY (7, 9), and 
GY (r, 9 => (> cx(1 —7)") cos nd. 
n=0 me 0 
It follows that the function U(r, #, 0), which except for a numerical factor is 
the fluid pressure on the disk, has a development of the same form as (VI.10). 


Here the singularity at the edge has exactly the same form as that obtained in 
[9] for the two-dimensional problem. 
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Appendix 
The purpose of this Appendix is to study the coefficients in the linear 
equations (III.10). At first glance these equations would appear to present 
difficulties since they require differentiations of the singular operator U. We 
shall show, however, that these difficulties are only illusory. 
To begin with we remark that the kernel G(r, 3, 0, 0, y) of U(r, 8, 0; f) can 
be evaluated explicitly. By [72] (p. 196) we have 


[oe 


| sane =F (Sotho) — Yoltw)) 
0 


m+ uel? 


where Y, is the singular Bessel function of first kind and S, is the Struve function. 
Thus by (11.2) we can write 


G(r,9, 0, 0,9) = Blu) = 4 —*2 [Sob n) +¥o (kn) +4 I(bH)], 


(A.1) 
=r + 0% — 2ro.cos(?—y). 

Let us now study the form of the functions f(7,¥), 7,(7, 0) appearing in 
(IV.10). As remarked in Section V they must all have the form 7”F(r) cosn@. 
Thus in (111.10) we are concerned with functions of the form, 

221 
x(r, = s | [ 0% Flo) cosmy Bu) de dy 


0 0 
2m1 


(A.2) -4( [ferro cos nad) dda) cos 


g=rto?—2rocosa, A(t) = Bir), 


and the problem is to differentiate such functions with respect to 7. We obtain 
the desired result from the following relation, easily verified by induction: 
[n/2] : — 
= > CB ee tary ope 
j=0 


o” 
er" 


Then from (A.2) we have 
[ 2a 1 


es sa an gi F(o) cosny B(2r—29c08y)"-*!-21dedyhcosn d 
or 
[n/2 2a 1 
= x ae [fe en21+* F(g) B"1(Q?) do dy. 
a 
=! 0 0 
But 


cos my + a, cos(m —2)y+-::, 


cos “my a = 
so that the only term remaining after the y integration is that for 7=0, and 


ra} y” 


(A.3) CP UNOie ye =( [Fle B” (0?) do) cosnt. 
0 
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Now the Struve function S,(z) is an analytic function of z, while Y(z) equals 
227 J,(z) log z plus an analytic function of 2’. Hence by (A.1) we have 


or 


A(t) =44+(Diage)loge + Dit + ees 
j=0 j= i= 


A”) (9?) = Fe a » a; 77+ >’ B; 0” log @ +y 7 2i+t, 
i=0 


j=—n 7 ==" 


It follows that the integrals (A.3) are all convergent. 


This work was sponsored in part by the Air Force Office of Scientific Research 
and in part by the Society of Naval Architects and Marine Engineers. 
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Diffraction by a Quarter-Plane 


JAMES RADLOW 


Communicated by M. M. SCHIFFER 


1. Introduction 

We report on our study [6] of the diffraction by a quarter-plane. This 
problem differs from its direct ancestor the Poincaré-Sommerfeld half-plane 
problem in being essentially three-dimensional in character. It has been known 
for some time as an open problem of diffraction theory. 

Our analysis is in essence a generalization of the function-theoretical method 
of WIENER & HoprF [9] from one to two complex variables. Relevant notations 
of two-variable function theory and the two-dimensional Laplace transformation 
are given in Section 2. We state the physical problem to be considered, together 
with the corresponding boundary problem, in Section 3. Next, in Section 4, we show 
that the boundary problem is formally equivalent to a transform equation 
relating two known and four unknown two-variable functions. 

The solution of the transform equation requires a certain two-variable factor- 
ization. This part of the analysis (Section 5) directly generalizes the basic factor- 
ization process of the one-variable WIENER-HopF theory. In Section 6, the unknown 
functions are expressed in terms of the factors obtained in Section 5, and it is verified 
that the conditions of Section 4 are met, 7.e. the transform equation is satisfied, 
and the functions involved possess certain required analyticity and integrability 
properties. We remark that the relationship between the factorization of Section 5 
and the solution given in Section 6 is less straightforward than the corresponding 
relationship in one dimension. 

We summarize the mathematical argument, in Section 7, by proving a theorem 
which asserts that the solution of the transform equation given in Section 6 corre- 
sponds to an actual (7.e., not merely a formal) solution of our boundary problem 
on the quarter-plane. Physical properties of the diffracted field are deduced in 
Section 8: we show that this field is an outgoing wave at infinity, and discuss 
its behavior near the edges and the corner of the diffracting screen. 


2. Notations 


We shall work in a Euclidean space F,: the schlicht finite s-space, where 
s—o-+17 and o=(0;, 6,), t=(t,, tT.) are real points of E,. The L, norm of an 
analytic function o(s)=(s,, S2) is defined by 


Ip@k={ flo +i alan ary) (2.1) 


10 
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and is in general bounded only for restricted domains of o. A typical domain 
of boundedness is a tube (see [2]): a tube 7(A) in s-space is a point set 


0,CA, = GO, Te OO (252) 


(j=1, 2) where A, called the basis of the tube, is any set in the o, o,-plane. For 
brevity, we agree that in this paper the subscript 7 is always to run over 
the values 1, 2, as in (2-2). 

If a function y/(s) is analytic and of bounded L, norm in a tube with basis 


Dw piee pep (2.3) 


then by a theorem of BocHNER (see [2], p. 738) it can be written as a sum of 
four functions @, (s), analytic and bounded in the respective octant-shaped tubes, 
or quarter-spaces, 7(D, 1), with 


DAS Oy, O2>Y2; (2.4.1) 
DD 28 Bae Om 02> Ye; (2.4.2) 
DBO NGS ON Ga Os, (2.4.3) 
D Ay 0, =; ee Oo: (2.4.4) 


The functions q,,(s) are uniquely determined up to additive constants. We note 
that the subscript in this paper will (except when otherwise specified) always 
run over the values 1, 2, 3, 4. 

Consider a function, say V(s) W(s), analytic and of bounded L, norm in a 
tube T(D) as above. This function is by BoCHNER’s theorem uniquely decompos- 


able into a sum of four functions A,,(s), analytic and bounded in the respective 
tubes T(D, n); we shall denote 


[V(s) W(s)]n = An(s), (2.5) 
[V(s) W(s) ]or = Ax (8) + Aa(s), (2.6) 
[V(s) W(s) lox = 41 (8) + Aa(s)- (2.7) 


The indicated decompositions are carried out by means of CAucHy’s formula, 
in two variables or one, as follows. Let Z,, Z, denote vertical contours from — 7 co 
to +7200; then we have, e.g. 


AG ee i: [ HAO taba dl (2.8) 


ca, 
for o;>Re(z;)>y, (see pe Also, 
V( (eam az. 
Ald Ad) may f Haste 29 
for o,>Re(z,)>y,, and 
Valse nV (See 
A, (s) + 44(s) = 24 vi Si Za) TG 24) des (2.10) 


for o.> Re (z.)>yp. 
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We denote by ¥ the two-dimensional Laplace transform (see [3]): given a 
function h(x,, x) of a suitable class, we define the Laplace transform H(s) by 


HS) =f [h(x = Jf h() x) exp(—s-x)dx,dx, 


with s - ¥=s,%, +S), s;=0,+71, as before. A tube, as defined above, is clearly 
a typical domain of convergence for the integral on the right, and therefore a 
typical analyticity domain for a Laplace transform H(s). The inverse trans- 
formation is defined by 


h(x) = 2H) =p ar} { #6 exp(s- x) ds, ds, 


SS 


where S,, S, denote vertical contours from —ico to +100. 


3. Formulation of the problem 
We consider the diffraction of an E-polarized harmonic plane wave by a 
perfectly conducting screen which occupies the quarter-plane x,=0, *,=0, %,=0. 
The incident plane wave is 


Ug (%1, Xp, X3) = EXP (— Ay % — Ay Xp — Az Xs), (34) 
where 
a, =tkh sind, cos o, Cp — 160s Us a3 =k COS Hy Sin Po 


with k=p—1q(p>0, q>0) and 0<%Xa/2, 0<MpSa/2. Our problem may 
alternately be taken to refer to the diffraction of plane sound waves, with velocity 
potential u,, by a perfectly reflecting quarter-infinite screen. 

Let the total field be w)+-w; then the scattered field u is a three-dimensional 
wave function which cancels “) on the screen. Next, denote the surface current 
lensity by { (x)= (a5 %): 

ea hes 


ou 
0% 


Ou 
%,=0-+ OX 


(3.2) 


t= 0— . 


It is known that /(x) must vanish off the screen, 7.e., for E[4,<0U%,<0], and 
we also ask that /(x) have the behavior of u)(x,, *,, 0) at infinity. 

Two more conditions on w are needed. Our assumption that #& is complex 
means (see the discussion in [J], pp. 154—55) that SOMMERFELD’s radiation 
condition may be replaced by the condition that w be outgoing at infinity, 2.e., 
that it have the behavior of the exponentially damped free-space GREEN’S function 
(e~'*®) (4 R), for large values of R=+ (x}-+ 23+ %3)%. Also, to ensure a finite 
flow of energy into the scattering region, we require that w(%,, x), 0) be of 
integrable square (as a function of two variables) at the origin. 

Finally, we shall in the interests of uniqueness ask that f(x) be of integrable 
square at the origin. This condition is based on our assumption that the current 
density, in a physically acceptable solution, has the behavior f~r~**° (with 
yi gtik 0O<d<4) near r=0. 
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We now have the following mixed boundary problem on the quarter-plane 
to be solved for the scattered field u(x, %2, %g): 


(As,r% + #) u=0; (3.3) 
W(%4,%q,0) =— e744 —%, for x20; (3.4) 
f(x) =0, for [m4 < 0v%,< 0]; (3.5) 
T(H) = Ole 26"); S iotilarge 470. (3.6) 
u—O(e4*/R), for large R, (3.7) 
provided x,;=-0, while if %3=0, 
U(X; %qs0) = Ole 2 x)) (3.7.4) 
for large values of |x| =-+ (27+ x3), with x= (%,, %,) any point in E[4%,<0U%,<0]; 
u (4%, %2, 0) ELe(|xj| S 4;) (3.8) 
for any finite d,, d,; Hx) €L,(0S x <0) 3.9) 


for any finite.c, C,. 

Conditions (3.3) to (3.9) completely describe our problem. The physically 
(and mathematically) reasonable assumption that the field scattered by the 
quarter-infinite screen x,=0, %,=0, %;=O0 is indistinguishable for large positive 
x, from that scattered by the semi-infinite screen 


yo 20 < maa}. | oy 2 ON = 0 (3.10) 
and indistinguishable for large. positive x, from that scattered by 
Lah = Oni (27 OO Aa a Oe eae el (3.11) 


is not a logically independent requirement. The solution derived in Sections 4, 5, 6, 
and proved in Section 7 to meet conditions (3.3) to (3.9), will be shown in Section 8 
to correspond to a field which reduces, for large positive x, (x), to that scattered 
by 2 (2%). 
4. The transform equation 
A formal solution of (3.3) which is outgoing at infinity, 7.e., satisfies (3.7), is 


W(%4, a) =e [ f Fl) K(s) exp|s+ «+ || ds,d5sq (4.1) 
ce 


2K(s) 


where s+ x*=S, %,+ 5s, %; S,, S, denote vertical contours from —ico to +100; 
F(s) =F (s,, Ss.) is an unknown function; and 

K(s) = (#2) (st + 3 + R®)-# (4.2) 
with the branch of the square root determined by the requirement that it reduce 
to + on the manifold s{+s3=0. From (4.1), we have 


U(X, %, 0) Serle Kis) exp (8-4) ds, ds, (4.3) 
Siecs 
f (4%, %_) = eae { [Fo exp(s- x) ds, dsp. (4.4) 


SiS. 
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Now (4.4) defines {(*) as the inverse Laplace transform of F(s). We shall, 
correspondingly, define F(s) as the Laplace transform of f(x). From (3.5), we 
then have 


F() =f J fl %) exp(—s- 2) dma, (4.5) 


and may show (using (3.6), (3.9), and the triangle inequality) that f(x) exp(—s- x) 
is of bounded L, norm, for o;>— Re(a;). It follows from PLANCHEREL’s theorem 
(see [3], Prop. 4.1) that F(s) is analytic in the quarter-space T(y, 1) with basis 
Mi AG) 6; == Ke(a,) (4.6) 
and of uniformly bounded L, norm for o€ (y, 1). 
Next, we see from (4.2) that K(s) is analytic and uniformly bounded in the 
tube 7(5) with basis 
B: oj| <b, (4.7) 
where 
b>0, (bf+03)?<g¢=|Imk|. (4.7.1) 


Since the intersection Bo(y,1) is nonempty, we conclude that F(s) K(s) is 
analytic in the tube with basis Bo (y, 1) and of uniformly bounded L, norm for 
o€[Bo(y,1)]. From BocHNnEer’s theorem (see Section 2), we have a unique additive 
decomposition 


F(s) K(s) = G(s) + Ge(s) + G(s) + Ga(s) (4.8) 


with G,,(s) analytic and bounded in a quarter-space 7(6, ”), where the basis 
(6, m) covers the m“ quadrant of the o, o,-plane and also contains the interior 
of Bo(y, 1). 

Notice now that (4.3) defines u(x, 0)=wu(x,, %2,0) as the inverse Laplace 
transform of (Ff K); alternatively, we have 


F (sys) = [a (4, Oy (4.9) | 


and we see, from (3.4), (3.7.1), (3.8), that the right side of (4.9) is a function 
of bounded L, norm for o€[Bo(y, 1)], which is in agreement with our result 
for (FAK). Let us write 


u(x, 0) = > u,,(%, 0) (4.10) 


where w,(x,0) vanishes outside the m quadrant of the x,%,-plane; i.e., if 
Y(x;) denotes Heaviside’s function (Y(x)=0, x;<0; Y(x)=1, %;>0), define 
u, (x, 0)= u(x, 0) V(x) Y(%2), v2 (x, 0) =u(x, 0) Y(— %) Y(%2), ug (%, 0) =u (4%, 0) 
Y(—%,) Y(— %2), “,(u, 0) =u(x, 0) Y(x,) Y(— x). Then, from (4.8), we have 

G, (s) = L[u, (*, 0)], (4.11) 
and, in particular, from (3.4) 


G,(s) = — (8, + 4) (82 + 4). (4.11.1) 


We have now applied all the conditions of our boundary problem and may 
summarize the results as follows. The right side of (4.1) represents a function 
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u (4%, X_, %3) which satisfies Eqs. (3.3) to (3.9), provided that F(s) is analytic 
in T(y, 1) and of uniformly bounded L, norm for a € (y, 1) (see (4.6)), and that 


F(s) K(s) = G(s) + Ge (s) + G(s) + G(s) (4.12) 


where the known functions K(s), G,(s) are given by (4.2), (4.11.1); and the 
unknown functions G,, G3, G, are according to (4.11) related to w(x, 0) for 
E[x,<0U%,<0]. Our problem therefore consists of solving the transform 
equation (4.12) for the unknown functions F, G,, G3, G4; it is understood these 
solutions must be analytic and of uniformly bounded L, norm in the respective 
domains o€ (y, 1) and o€ (6, n) for n=2, 3, 4. F is required to be analytic in 
T(y, 1) in order that condition (3.5) be met. 

We may instructively compare the present problem with the simpler one 
we considered in [7]: a vector problem of diffraction by a semi-infinite screen 
was there reduced, by the introduction of double Laplace transforms, to the 
problem of solving a transform equation (see [7], Eq. (2.5)) whose form was 
precisely that of (4.12). There were, however, two decisive simplifications, whose 
effect may be understood by our supposing, for the moment, that (4.12) is to 
be solved under the following conditions: (1) the known functions include not 
only K, G,, but also G,; (2) F is required to be analytic only for 6, > — %, | o2|<a, 
where «20, #)>0. This simplified problem may be solved in the manner of 
[7], namely, on the basis of the factorization 


| K(s) = (i/2) M, (s) M_(S) (4.13) 
aye M, (5) = [sy + i(s + 29) (4.14) 
M_(s) =[—i(8 + A). (4.15) 


Let k=p—v71q as before, and make the choice 
($+ #9 = +h; 


then the factors M,(s), M_(s) are analytic (and nonzero), respectively, for 
0,20, |o,|<q; 0,0, |o,.|<g. We solve (4.12), under the simplified conditions 
assumed, by arguments of the type given in [7], with the results: 


ye ak FG Saas) 
Fon eeae)  lad ay 
Gals) +Ga(s) =—M_(s) [SEE (4.17) 


where the bracket-subscript notation is that of (2.6); it is understood that (2.6) 
has the counterpart 


[V(s) W(s) Jor = Aa(s) + Aa(s). 
Notice that the right side of (4.16) is analytic, as required, for 0,20, |o.|< a. 
We compare the simplified problem solved by (4.16), (4.17) with the problem 


we have actually to consider by first observing that we have more information 
in the simplified problem, 7.e., we are given 


[2*(s) K(s) lox = G(s) + Gy (s) (4.18) 
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with both G,, G, known, while in our actual problem we have only 
[F(s) K(s)], = G(s). (4.19) 


Also, the simplified problem asks, in effect, that the complete additive decom- 
| position of f(s) have the form 


F(s) =K(s) + Fy(s), (4.20) 
| while our actual problem makes the more stringent demand that F(s) have the 


complete decomposition 
B(s) = Fils). (4.21) 


5. Factorization 
Our solution is based on a two-variable factorization of K(s) (see (4.2)) : 


K(s) = (/2) [T M,,(s) (5.4) 


with the M,(s) analytic and nonvanishing in respective octant-shaped tubes 
(quarter-spaces) 7(B,) whose bases (B,m) have B as intersection (see (4.7) 
for B). The (B, ) are given by 


(BAS 6,> —b;, Op > — de, (5:24) 
(Bis aps by) C3 > 05; (38232) 
(B,3):~ 0 < by, dg < by, (5.2.3) 
(B,4)s 6, >= 0,, Oo < dy. (5.2.4) 


The problem of deriving the factorization (5.1) generalizes the classical one- 
variable problem of WIENER & Hopr [9]. As in the one-variable case, we may 
take our problem to be that of additively decomposing log K(s). 

The additive decomposition problem is solved as follows. Let us allow the 
parameter k to vary over a disk 7, with center k=(p, —q) and radius 7. From 
the analyticity of K(s), 7 may be so chosen that A(s) is an analytic function of 
s, k throughout 7,. It is then permissible to differentiate A(s) with respect to 
k (k€n,), and to form the logarithmic derivative (the prime denotes differentiation 


with respect to k): 
K’(s) eas 


“K(s) ~ s?-+s3+h? 6-3) 


whose analyticity domain 7(B) is the same as that of A(s), and which we shall 
take as the function to be decomposed additively. The decomposition is to 
have the form 

goa Mn (s) (5.4) 


with the u‘* term on the right analytic in T(B, ). Given (5.4), we would expect 
to be able to justify an integration with respect to k, with the result 
4 


log K(s) = Slog M,, (3), (5.5) 


n=l 


146 James RaDLow: 


which is equivalent (after constants of integration are properly accounted for) | 
to (5.1). 
We begin the calculation by observing that (5.3) may be rewritten as 


Ris) th (2) 6 
#1) _ ih gpg (bal) (6.6 


with |x| =-+ (x3+28)*; both sides of the equation are analytic in T(B). Te 
make use of (5.6) in deriving the required result of the form (5.4), we express 
F as the sum of its restrictions to the four quadrants of the x, x.-plane. We have 


G=>e, (5.7) 
where ie 
GL, = ff dxdt, V (ens %) Y (Ena 2) exP(— 8-2) (5.8) 
and ~ 
&qeA (a =A 4) (5.91) 
Se tlle ee 
Eno = +1 oo (5.9.2) 


| (n = 3,4). 


VYi=O0e £20; VO 0. 
It then follows from PaRsEVAL’s theorem, 7.e. the result (see [3]) 
L(g f)=[g] ar Z[f), (5.9.3) 


where += denotes convolution, that the operations Y, map H/)(k|x|) into 
functions analytic in the octant-shaped tubes 7(B, 7). The required additive 
decomposition is therefore given by 


a ik < 
Rey a aL la) (5.10) 


The first term in the summation on the right side of (5.10) is the known 
transform (see [8], Eq. (3.5)) 


a Y, [HY (k|x|)] = Als, 52) = — 4 K2(s) [—7+5,,+ 5, ] (5.14) 
P= t* {1+ (22/20) log [(1/R) (s; + &)]} (5.11.1) 

and 
t; = (s} + A)8, (5.44.2) 


Next, each of the remaining three transforms on the right side of (5.10) may 
be reduced to an integral over the first quadrant of the x, X-plane by changing 


the sign of the s-component corresponding to the negative x-component. Our 
result, in the notations of (5.9.1), (5.9.2), is 


L,, (Hy (A|%|)| = Ales $1) E,2S0) (5.42) 
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with A(s) defined by (5.11). The decomposition (5.10), of the desired form (5.4), 
is now explicit. 
Integrating (5.10) with respect to k, we have 


log K(s) = ES fea 2) (k|x|)] dk + o9(s) (5.13) 


Where «% (s) is a constant of integration (i.e., independent of k). From (4.2), 
log K(s) = (¢ /2) — log 2 — (§) log (st + 3 + #?) (5.14) 


in T(B). Comparing (5.13), (5.14), we infer that «,(s) is everywhere equal to 
(¢z/2) —log2, and the factorization 


= (2) LT exp {(— 4) J & 2, [HS (| x1)] a} (5.15) 


of the desired form (5.1) follows from (5.13). The functions M,,(s) of (5.1) are 
then given by 


k 
M,,(s) = exp {(— 1/4) J k &, [Hg 2) (| x|)] aah (5.16) 


(see (5.41), (5.12)) and, in accordance with our argument, are based on PARSEVAL’S 
theorem, analytic in the octant-shaped tubes 7(B,n); they are, as required, 
nonvanishing. Observe also that we have (with M, (s), M_(s) defined by (4.14), 
(4.15)) 


M,(s) M,(s) =™,(s), (5.17) 
iM, (s) M,(s) = M_(s). (5.18) 
If we denote 
M? (s) =[se-+ 4 (s+ )J-4, (5.19) 
M-(s) = [sy — é(s} + 2)4}-4, (5.20) 
we have, further, 
M,(s) M,(s) =M*(s), (5.21) 
M,(s) M,(s) =M~(s). (5.22) 


Growth estimates for the M,,(s) may be derived by carrying out the inte- 
gration on the right side of (5.16). The results obtained are 


M,,(s) ~ (si + s3)~* (5.23) 
(B, n). 


for large 


6. Solution of the transform equation; verification 
The transform equation (4.12) poses the problem of determining the four 
unknown functions F, G,, G,, G, from the information that 


[F(s) K(s)h = Gi(s), (6.1) 


e., from the transform equivalent of the physical condition (3.4). We also 
require that F(s) be analytic in T(y, 1), and of bounded L, norm for o€(y, 1), 
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so that its additive decomposition is given in full by 


F(s) =A(s), (6.2) 


which is the equivalent of the physical condition (3.5). | 
We assert that a function F(s) which satisfies (6.1), (6.2), and which corres 
ponds, by means of (4.1), to a solution of the boundary problem stated by 


Eqs. (3.3) to (3.9), is given by 


ee BOR 6: 
FAS) pred ene Om 
a ae sh (s) Mg (si, —4s) 
Et Ty oa ay 6.4 
a0) M,(—4, S2) M3(— 2, TSE M,(s) |1 (6.4) 


where the notation is that of (2.5), (2.6). This means, according to (2.8), (2.9), 
that the expressions on the right sides of (6.3), (6.4) are explicit integrals. The 
assertion we make in (6.3), (6.4) is then that 
F(s) = 2 , F, (2, $2) d% (6.3.1) 
$;—4) 1 


274 


Zz, 


for o,> Re(z,) > — Re(a,), with Z, as in (2.9), and that 
= G, (2) M2 (z) Ms(4, —42) d% a2, ‘A 
Sees Gua Jit S14 —%) (Sg—42) Mz(—4, 2 M 3(—@,, —42) M3 (z) M,(z) ie 


for o,>Re(z,) > —a,, with Z,, Z, as in (2.8). G,(s) is given by (4.11.1), and the 
functions M,,(s) by (5.16). 


Next, we assert that the functions G,(s), G,(s) are given by 
ee M,(s) Mg(s,, —a) 
@, (s) i. Gi (s) ea Sa) netic —@y) 1 ; (6.5) 
= M3(—4,, $2) M,(s) 
o (s) = a (s) M3(—a, —4,) M, (Sy, —@q) 1 (Gm 
d that . es 
eae G,(s) = F(s) K(s) — G(s) — G,(s) — Ga(s) (6.7) 


with F, G,, G, given by (6.3), (6.5), (6.6) respectively. 

To verify our assertions (6.3), (6.5), (6.6), (6.7), we must show that one of 
the four (independently of the other three) implies (6.1), (6.2); and we must 
also show that each of the four assertions implies the other three. It is however 
clear that (6.7) is correct if (6.3), (6.5), (6.6) are, so that it will be enough to 
verify the results for F, G,, G,. 

Our verification argument will be organized as follows: Section 6.1 gives the 
verification for G,(s), omitting a certain proof in order to maintain the con- 
tinuity of the essential argument; Section 6.2 contains the omitted proof; the 
results for F(s), G,(s) are verified in Sections 6.3, 6.4. The argument for F(s) will 
not only show that (6.3) is equivalent to (6.5) but also that (6.3), independently 

f (6.5), implies both (6.1) and (6.2). This is more than is logically required, 
and we therefore observe that the additional material is included for the insight 
it gives into the structure of our solution. 


Diffraction by a quarter-plane 149 


6.1. Now the result of substituting G,, as given by (6.5), into the transform 
equation (4.12) is 
SpA Ae GAS NLS) (sta) 
F(s) K(s ; —— e 
ead) M,(—4,, 53) M3(—a,, —4a) 


feta lS) 1 sGrais) - (6.8) 


From (6.8) we find, successively, 
G(s) M,(s) M;(s,, — 4) 


Dee Oo ARE Sea Tce meee (6.9) 

[FA] = [(F K)o3]o¢ = G(s). (6.10) 

The same result is obtained by taking components in the other order: we have 
[F K],+ = Gy (s) + G,(s), (6.11) 

[PK], =F )otlot = Gi (3), (6.12) 


so that (6.5) is seen to imply (6.1). 
Next, we show that (6.5) implies (6.2). Dividing both sides of (6.8) by M, My, 
we have from (5.1) (with the M,,(s) given by (5.16)) 
GA M, (Gi, as) Scme | EGGARG, 6.13) 
M,(—4, 52) M3(—a,, —a,) M,M, ' M,M, ° 


(S)F M, M,= 


We then have additive decompositions of the following forms (A,, Ay, C3, Cy 
are unspecified analytic functions, and the B, are known analytic functions; 
subscripts denote analyticity domains): 


(i/2) FM,M,=A,+ 4,, (6.14) 


G, M, M3(s,, —4p) 
M,(—4,, S2) M3(—a,, —a,) M, M, 


(Gs + Gy) (M3M,)1=C3+ Cy. (6.16) 


The decompositions (6.14), (6.15), (6.16) are unique (this is proved in Section 6.2). 
From the uniqueness, we have 


Ap Ag = BL Bi, (6.17) 
CPC eB, (6.18) 


= B,+-B,+ B,+ By, (6.15) 


where, since the B,, are known, we may solve for F from (6.14), (6.17) and 
may solve for G3, G, from (6.16), (6.18). We have 


(1/2) F¥M,M,=8,+ By, (6.19) 


which gives (using uniqueness again, which in the present case amounts to 
carrying out a one-variable WIENER-Hopr argument on the variable sj) 


= 
LS) Sa) 


B,(s) 
mata 28) 


with B,(s) given, according to (6.15), by 
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Our statement that (6.3) implies (6.2) now follows from (6.20) by inspection, 
since M,(s) is analytic and nonvanishing in T(B, 1), while B[,/M,],+ is analytic 
in T[(y, 1)>(B, 1)]. Comparing (6.21) with (6.4), we see that B, =A, 1.e., our 
result (6.5) for G, implies the result (6.3) for F(s). Notice, again, that the result 
for F(s) is explicit: see (6.3.1), (6.4.1). 

Finally, we use symmetry considerations to show that we must have 


G4 (S11 503 G1, 2o) = Ge (Sq; $15 Ges 2); (6.22) 
and, since it follows from the definitions of the M,,(s) (see (5.16)) that 
MM, (55:8) = 0 ls 3a) (6.23) 
M;(s,, — 4) = Mg(— 4, 52), (6.24) 
M, (— a373,)=M,6,,—@), (6.25) 


we conclude that (6.5) implies (6.6). 


6.2. We now complete the argument of Section 6.1 by proving that the de-. 
compositions (6.14), (6.15), (6.16) are unique. A proof that F(s), as given by 
(6.3), is of bounded L, norm for o € (y, 1) will be included. 

Consider the left side of (6.15): G,(s) is given by (4.11.1), and the growths 
of the M,,(s) are given by (5.23). Denoting the left side of (6.15) by P(s), we 


then have P(s) =0(|s|~4) (6.26) 


oe Is] =|[s2+ 333 (6.26.1) 


for large |s| in the analyticity domain of P(s), 7.e. in the tube T[Bo(y, 1)]. 
P(s) is therefore of bounded L, norm for o€[Bo(y, 1)] and has according to 
BOocCHNER’S theorem (see Section 2) a unique decomposition. 

As to (6.14), (6.16), it is sufficient to show that one of these decompositions 
is unique. Consider the function on the left side of (6.14). We had actually to 
assume the boundedness of this function’s L, norm in order to assert the uni- 
queness of the decomposition (6.14). Since this process led to (6.20), we are 
logically obliged to show that 

|F M, My |, < const. (6.27) 


for o€ [(y, 1)>(B, 1)9(B, 2)], with F(s) given by (6.20) (or, equivalently, by 
(6.3)). 

To prove that (6.27) holds, we remark first that (6.3) may be shown (and is 
shown, in Section 6.3; see Eqs. (6.35), (6.36)) to imply 


2255) G(s) M,(s) Mg (s,, —a,) 
JIVE a= 2 1 2 3\S1 2 
i ; M,(—4y, 8.) Mg(—a,, —a@,) M3 (s) M,(s) pe ; 028) 
From (6.28), we find 4 

FM,M,=0(|s|~#) (6.29) 
for large |s| in T[(y, 1)>(B, 1)0(B, 2)], and we see that (6.27) is satisfied. Also 
from (6.29) and (5.23) F ay 

(s) =O(|s|~4), (6.30) 


so that the function F(s) given by (6.3) is, as required, of bounded L, norm for 
aE (y, 1). 
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6.3. We have shown in Sections 6.1, 6.2 that F(s), as given by (6.3), is analytic 
in T(y, 1) and of bounded L, norm for o€(y, 1). This means that (6.3) implies 
(6.1). In this section, we prove that (6.3) implies (6.2) and also that (6.3) 
implies (6.5). The proof that (6.3) implies (6.6) is given in Section 6.4. 


Let us begin by writing z accordance with the notation of (2.6), (2.9)) 


[FM, Moa: nes ai | F(4,s 1 (21, Sg) My (2, 52) d% (6.31) 


$1 — 24 


for o,> Re(z,) > — Re(a,); Z, denotes a vertical contour from —ioo to +/o0. 
Making use of (6.3), (6.4), we rewrite (6.31) as 


W,, So) dw, dz. 
[F M, M,)].+ = ay (21)? abate Ge ai reas (w,, 5) (6.52) 


Z, W, 


with o,> Re (z,) > Re (w,) > — Re(a,); W, denotes a vertical contour from — ico 
to +100. Carrying out the integration with respect to z,, we have, upon deform- 
ing the contour to the left, 

Gao ee (6.33) 


24 . SD, 
W, 


for o,> Re (w,) > — Re(a,). A second integration gives 

[F M, M,),x = — 21 P, (81, 89). (6.34) 
But (6.34) implies that we have 

FM,M,=— 21/4 (s) + 4(s)] (6.35) 


with F,(s) an analytic function in 7(B, 2), and we see from (6.4) that 


=a G, (s) Mz (s) M3 (s;, —42) / 
etiyeat M,(—4, 52) M;(—4a,, —4@2) M;(s) M,(s) le: e0) 
Again, (6.36) implies that we have 
G, (s) M,(s) M3 (s,, —42) ay S P, (s) (6.37) 


M,(—4,, S2) M3(—4,, —ay) M3 (s) My(s 


where P,, P, are as before, and P;, P, are functions analytic in T(B, 3), T(B, 4) 
respectively. 

Now multiply both sides of (6.35) by (t/2) Mg(s) M,(s): the result, upon’ 
making use of (5.1) on the left side and of (6.37) on the right side, is 


G, M, M;(s,, —42) a ( 


eK 
M,(—4a, $2) Mg(—a, —4e) 


P,+P,) M,M,. (6.38) 


To prove that (6.38) implies (6.2), we proceed exactly as we did in showing 
that (6.8) implies (6.2). Since (6.3) has been shown to imply (6.38), our proof 
that (6.3) implies (6.2) is complete. 
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Notice next that the right side of (6.38) is equal to G,+G,+G,+G,. It 


is then clear that we have 


G, M, M;(s,, —42) vt 2 
Ml—a,, Sa) Ms(—4, —4s) ids Gs (6.39) 
— (P, + 4) M,M,=G,+G, (6.40) 


with (6.39) equivalent to (6.5), so that we have proved that (6.3) implies (6.5). 


6.4. In this section we prove that (6.3) and (6.6) are equivalent. Since we 
have demonstrated the equivalence of (6.3) and (6.5) in Sections 6.1, 6.2, 6.3, and 
have also shown in those sections that the results (6.3), (6.5) separately imply 
(6.1), (6.2), the arguments of this section will complete our verification. 


We begin by noticing that the symmetry of our problem gives 
F(s1, $2; %, 42) =F (Sq, 813 4g, &). (6.41) 


Taking account of (6.41) and the symmetry properties of the M, (s) (see (6.23), 
(6.24), (6.25)), we find that an alternate (and completely equivalent) form of 
our result (6.3) is 


— S2tALe(s) 
US Ne are Ea li ces 
with 
= G,(s) M3(—4,, Sz) M,(s) 
# (s) ae M;(—4, —4@z) Mz(s) Ma(s, —@,) Ja (6-49) 


The notation in the last two equations is that of (2.7), (2.5), which means (accord- 
ing to (2.10), (2.8)) that (6.42), (6.43) are explicit integral representations. 

An argument of precisely the type we used to show that (6.3) implies (6.38) 
will now show that (6.42) implies 


G, (s) M3(—4a,, Sy) M,(s) 


IOI = 
M;(—a,, —a,) M,(s,, —4a) 


(2 + Qs) M. Ms (6.44) 


where Q., Qs are analytic functions in T(B, 2), T(B, 3). We observe that (6.44) 
implies (6.2). 
ae we have as consequences of (6.44) the equations (analogous to (6.39), 
(6.40) 
G,(s) M3(—a,, v2) 2a 
M3(—4,, —4@,) M,(s,, —ag) 


— (Q2+ Q;) M,M,=G,+G; : (6.46) 
Since (6.45) is equivalent to (6.6), we have shown that (6.42) implies (6.6); also, 


since (6.3) and (6.42) are equivalent, we conclude that (6.3) implies (6.6). 


To show that (6.6) implies (6.42) (and therefore that it implies (6. 3)); we 


make use of (6.6) in our ae equation (4.12). The result, after division 
of both sides by M, Mg, 


=G6,+G, (6.45) 


v = 
SWE Von ——— So) M Geese 
(3) Wes M, M, M. Ges ar Mls =e M, ™, . (6.47) 
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An argument which follows the same course as that which led from (6.13) to 

6.20) th i ‘ 4 

( ) then gives us Fls) = Te [2 (s) 
M,(s) |M,(s) 


(6.48) 


with D, (s) = R, (s) (see (6.43)). The right side of (6.48) is therefore precisely 


the right side of (6.42), and (6.42) in turn is equivalent to (6.3). Our proof that 
(6.3) and (6.6) imply one another is complete. 


7. Solution of the boundary problem 


The results of Section 6 yield a solution of the mixed boundary problem on 
the quarter-plane formulated in Section 3. We have the following 


Theorem. The expression (4.1), with K(s) given by (4.2), represents a solution 
to the boundary problem stated in Eqs. (3.3) through (3.9), provided that F(s) is 
given by (6.3) (or, equivalently, by (6.42)). 


Proof. We have shown in Section 6 that F(s) is analytic in T(y, 1) (see (4.6) 
for (y, 1)) and that it is (according to (6.30)) of bounded L, norm for o€ (y, 1); 
also, we observed in Section 4 that K(s) is analytic in T(B) (see (4.7) for B) and 
uniformly bounded for o€ B. Since [2K(s)]+ has negative real part (see the 
remark following (4.2)), we conclude that the integral on the right side of (4.1) 
converges uniformly in any region interior to 7[Bo(y,1)], as do the integrals 
representing the second derivatives of wu. It follows that (4.1) represents an 
actual solution of (3.3), if F(s) is given by (6.3) (or by (6.42)). 


Since (4.1) represents sea Ks, Ks) we-have 
U(%, X_, 0) =, (ni) ae | [FO s) K(s) exp(s+ x) ds, ds, (7.1) 


by dominated convergence. For rp 0, 


es ) K(s)],exp(s- x) ds, ds, (7.2) 


UW Kas Xs Ol Gna? 


oC {Fo ())nexp(s-%) 45,ds,=0 es 


for n=2, 3, 4 and “E 0; the notation [F K], is defined by (2.5). We proved 
in Sec. 6.3 that our result for F(s) (i.¢., either (6.3) or its equivalent (6.42)) 
implies (6.1), and we may eae: use (6.1) together with the explicit form 
(4.11.1) of G,(s) to rewrite (7.2) a 


exp(s- %) ds, ds. 
u(x, X,0 (221)? af ie eta Ss. 3+) 2 (73) 
eS 


Since the right side of (7.3) is equal to —exp(— a, x; — ay %y) for x,20, we have 
shown that (3.4) is satisfied. 

To verify (3.5), we remark as we did above that F(s) is, from (6.30), of bounded 
L, norm for o€(y,1). PLANCHEREL’s theorem (see [3], Prop. 4.1) therefore 
shows that F(s) corresponds, under Laplace inversion, to a function /(x) such that 


I(x) exp(—s- hb =[F()lb (7.4) 
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for o€(y, 1). The function f(x) is given (with equality in the mean-square sense) by 


(2) =20F)) = Gane [ [FO exp) dsdso, (755) 

ae? 
so that (4.4) is rigorously established (7.¢., it is not merely a formal represen- 
tation), provided that F(s) is given by (6.3) or by (6.42). Also, our function 
F(s) has been shown to be analytic in T(y, 1), so that the inverse Laplace trans- 
form on the right side of (7.5) vanishes for E[x,;<0 Ux,<0], and (3.5) holds. 


Since /(x) has been shown to exist, we may verify (3.6) by showing that 


lim (s, + 4) (S, + 42) F(s) = const. (7.6) 
or, equivalently, that 
lim (s,; + 4) (Sp + a) F(s) K(s) = const. (7.7) 


It is understood that the results (7.6), (7.7) must be independent of the order 
in which the limits are taken. 


Notice first that 


tim (5 + a) F(s) K(s) = lim (5; + a) [Gi(s) + G69) (7.8) 
since G,, Gg are analytic at s, = —a,, so that 
_lim, (s, + a) (Gels) + Ga(s)] = 0. (7.8.1) 


Making use of (6.6) and (4.11.1) in (7.8), we have 


i z fe\t a oe M3(—4, $2) My(—4, 52) 
Race ee ACS) (Se+4@,) M;(—a,, —a_) M,(—a,, —@g) ’ (7.9) 
and it follows that 
tim (5, +4) (83 + a2) F(s) K(s) = — 1. (7.10) 


If we take limits in the other order, we have 


Lim (s_+ a2) F(s) K(s) = lim (sy + ag) [Gy(s) + G, (s)] (7.11) 


Sg a, a= 8, 


and therefore find, using (6.5) and (4.11.1), that 


aS + 4) (S. + a) F(s) K(s) =—14. (7.12) 


The results (7.10), (7.12) are together equivalent to (7.7), as well as to (7.6), 
and we conclude that (3.6) is satisfied. We remark that (3.7), (3.7.4) may be 
deduced from (3.6); these conditions may also be verified directly by an argument 
similar to that we have just given. An alternate verification of (3.7), (3.7.49 
will be found in Section 8. 
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Finally, we make use of (6.30) to prove that f(x) is of integrable square at 
the origin (7.e., (3.9) holds), and it is true a fortiori that (3.8) holds. In Section 8, 
we give a more detailed account of the behavior of the functions f(x), w(x,, %>, 0) 
near the origin. 


8. Properties of the diffracted field 


In this section we discuss the behavior of the diffracted field at infinity, 
at the edges and at the corner of our quarter-infinite screen. 


The total field (incident field plus diffracted field) in our problem is 
Utot(%1, Xp» %3) = Ug (%y» Xe, Xy) + U(X, Xp, Xp) (8.1) 


with wu, u given by (3.1), (4.1) respectively. We proved in Section 7 that (4.1) 
represents a solution to the boundary problem stated in Eqs. (3.3) through (3.9), 
provided that F(s) is given by (6.3) or, equivalently, by (6.42). Notice now that 
(4.1) gives (by an application of PARSEVAL’s theorem: see (5.9.3)) 


(79, +8) exp(—ik|R r° ax? dx 
U(X; Xo; aes eee P( |R— a I) a 2 (8.2) 


where /(x) is given by (7.5), and 


| R— RO] = + [(m— 28)® + (om — 28)? + (or). 


Physically, the statement made by (8.2) is the evident one that the diffracted 
field results from a distribution of point sources, with density f(x), over the 
quarter-plane %,=0, %,20, *;=0. It follows from the representation (8.2) that 
the diffracted field is outgoing at infinity; this, as we remarked in the third 
paragraph of Section 3, is equivalent to Sommerfeld’s radiation condition. It is 
also clear from (8.2) that our solution meets conditions (3.7), (3.7.1). 

Next, we consider the diffracted field along the edges, but at a great distance 
from the corner, of the quarter-infinite screen. We shall show that our result 
(8.2) for the diffracted field reduces, for large positive x, (%,), to the field scattered 
by 2,(2,) (see (3.10), (3.11) for the semi-infinite screens 2, 23). 

The diffracted fields u;, uw; corresponding to scattering by 2, 2, are given 
by expressions analogous to (8.2): 


gay exp(—7tk |R— is ano ax 
Uy (x1, Xo, i) dy i Px = 2) ane ie U a (8.3) 
0 —oco 
x9, x8) exp(—thk| R—R®|) dx? dx8 
Uy (%1) Xo Xs) = efaal Pu ( 1, 3) arene 4 |) ec (8.4) 


—co 0 
where (x), Yu(*) are current densities on 2, 2}; it is understood that the 
screens X,, , are perfectly conducting and that the incident field is the plane 
wave uw, of (3.1). Applying the appropriate boundary condition (that the total 
field vanish on X,, 2,), we have, as equations from which gy, yy, may be deter- 
11% 
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mined, 00" 60 dx dx 
1 ‘i 0 49) exp(—tk|v—/r |) dx} dx 
exp(— 4) fl | pr (4a, 48) eT ectes Ceat} (8.5.1) 
6 =69 
{ ro gee —ik\|v—vr|) dx? dx 
exp(—a- a) = fa] pr (44, #2) ae 1) 67484 oo (Sa 
—c 0 


where @- %=@,%,+45%q, |r —1| =+ [(% — 24)? + (%2 — %2) 7). 

We shall now calculate the Laplace transforms ®;(s), ®j;(s) of the functions 
x(x), Pr(x) which satisfy (8.5), (8.6) and shall show that F(s), as given by (6.3), 
reduces to @,(s), ®;(s) in the appropriate limiting cases. To begin with it is 
clear from physical considerations that the solution g;(x) of (8.5.1) has the form 


Pr (%1, Xa) = EXP(— A %) Yr (2) (8.6.1) 
and also that the solution gy; (x) of (8.5.2) has the form 


Prt (%1, %2) = EXP(— 4g X) Pr (%)- (8.6.2) 


The function w;(%2) is then to be determined from the equation 
CXP (— ay %_) = (6/4) J (x9) Ish (Ay |X» aX x9|) d x9 (8.7.4) 
where k, = (k2+ aj)?; and wy, (x,) is to be determined from 


exp(— a %4) = (/4) J yu (xq) Hy” (Fo | x ins xy|) a xy (8.7.2) 
where ky= (k?+ a3)}. 
The Laplace transforms 


YF (s.) = fe (x5) exp (— Sy %>) dX, (8.8.1) 
Ps (s;) = Jon (x,) exp(— s, 4) dx, (8.8.2) 

of the solutions of (8.7.1), (8.7.2) are given by 
%, (Sz) = 24 (Sy + aa) 1 (— a, — 1 hy)4 (Sy +7 hy}, (3.9.1) 
Yrs (51) = 26 (8, a) Gay Se ee ea (8.9.2) 

so that (from (8.6.1), (8.6.2) 

Py (s) = (5, + a) * 4 (s3), (8.10.1) 
Dry (s) = (83 + ay) 2 My (5,). (8.10.2) 


The results (8.10.1), (8.10.2) are known to be the transforms of the unique 
physically acceptable current densities, in the respective cases of diffraction by 
the perfectly conducting semi-infinite screens 21, 2,. Our task is to show thas 
the function f(x) given by (7.5) is indistinguishable from Yr (Yu) for sufficiently 
large positive x, (x,). This means, according to (8.6.1), (8.6.2), that we must 
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show that 


_ Tim exp (@ 4) f(x) = ile), (8.11.1) 

i lim. exp (42%) f(%) = yu (24). (8.11.2) 
or, equivalently, that 

lim (5, + 4) Pls) = %(s), (8.12.1) 
pF lim  (s, + a2) F(s) = Px (s1), (3.12.2) 


| 
where F(s) is given by (6.3) (or (6.42)) and Y%, Wy, by (8.9.1), (8.9.2). 
| From (7.9), we have (in the notation of (5.21), (5.22)) 


lim (s; + a) F(s) BU Sastsda) (8.13) 


S.A a M*(— ay, Sg) M~(—a,, —a,) © 


Since we have (with k= (k?+ aj)}) 
M*(— ay, S$.) = (Sg +7 h,)—3, (8.14.1) 
M~(— a, — a.) = (— a, —tk,)-3, (8.14.2) 


we see that (8.13) is equivalent to (8.12.1), with ¥%(s,) given by (8.9.1). A result 
analogous to (7.9) follows from (6.5): 
me, M,(s1, —4) M5 (51, —42) 


Pun ase z 42) BIS) iS) = (si +) M,(—4,, a) M, (=, =a) : Ce >) 


and we see from this that 


lim (s+ 4) F(s) 2S a) 7 (8.16) 


Sp ay M, (51, —4) M_(—a,, —a,) 


with M,, M_ defined by (4.14), (4.15), and related to the M,, by (5.17), (5.18). 
We also have (with k,= (k?+ a3)*) 


M.,, (81, — a2) = (8, +4 e)~4, (8.17.1) 
M(— a4. 2— &) = (— aq —th,)=4, (8.17.2) 


so that (8.16) is equivalent to (8.12.2) with ¥%;(s,) given by (8.9.2). 

The argument we have given shows that the diffracted field w of (8.2) reduces, 
for large positive x, (x2), to the diffracted field w of (8.2) (wy of (8.3)) correspond- 
ing to the incidence of #) on 2 (24). Also, the diffracted fields u;, uy are the 
unique physically acceptable fields associated with the diffraction by 2,, 2 
respectively. This means that w;, “;; have the known energy-conserving behavior 
near the respective edges %,=0, x,;=O (i.¢., uy~7? near %,=0, uy~7 near 
x,—=0). The behavior of our diffracted field w (as given by (8.2), with f(x) given 
by (7.5)) near the edges of the quarter-infinite screen may accordingly be described 
as follows: w~r? near x,—0, for sufficiently large x.>0; and w~7} near x,=0, 
for sufficiently large x,>0. Further, the current density f(x) (as given by (7.5)) 
has the corresponding v~? behavior near either edge, at sufficiently large distances 
from the corner. 
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We conclude by discussing the behavior of the diffracted field near the corner 
(and in the plane) of our quarter-infinite screen. The quantities of interest will ; 
be u(x,, %,, 0), evaluated as r= (x?-+ x3)? approaches zero along a path exterior 
to the screen, and /(x,, x2), evaluated as y approaches zero along a path interior | 
to the screen. 

From (7.1), the Laplace transform of u(x,, %2,,0) is PK; we accordingly 
deduce the behavior of w near the origin from that of (s,s,/ A) near infinity. 
Results for two special modes of approach follow at once: from (6.5), the result 
as y->0 along the negative x,-axis is seen to be 


u (x1, X2, 0) ~74, (8.18) 


and we also see, from (6.6), that the same result is obtained as y->0 along the 
negative x-axis. The corresponding results for /(%,, %2) are clear: we have 


f (x1, %_) ~r—* (8.19) 


as r->0 along the positive x-axis or along the positive x,-axis. 


An approximation method developed by NoBLe [4] has led to results which 
differ from (8.18), (8.19): the exponents obtained in [4] are .30, —.70 rather 
than .25, —.75. We have, however, been informed [5] that the exponents obtained 
by the method of [4] vary with the mode of approach to the corner and, in 
particular, that the results mentioned are valid for an approach along x,= x 
(rather than along x,=0 or %,=0). A similar variation in field behavior may 
(or may not) be discernible in the course of an exact analysis. We are presently 
continuing our investigation in an effort to clarify this matter. 


This research was supported by Contract No. AF 19 (604) 5238 in the years 
1956—58, and by Grumman Aircraft Eng. Corp. in 1958—59. 
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Fehlerabschatzungs- und Hindeutigheitssatze 
fir Integro-Differentialgleichungen 


KARL NICKEL 


Vorgelegt von L. COLLATZ 


Ubersicht 
In der folgenden Untersuchung werden die Lésungen w(t) der allgemeinen 


_ impliziten Operatoren-Differentialgleichung 


LEYS; £4) = 0 “In 0 <tr (1) 
unter der Anfangsbedingung 
u(0)=7 (2) 


abgeschatzt und Einzigkeitskriterien fiir das Anfangswertproblem (1), (2) an- 
gegeben. Darin sei Fu ein Operator, der eine auf 0<¢tS<T stetige Funktion u 
in eine dort definierte Funktion v=Fw transformiert; die iiblichen Integral- 
operatoren besitzen im allgemeinen diese Eigenschaft. Uber die Lésbarkeit des 
Problems (14), (2) werden jedoch keine Aussagen gemacht werden. 

Das Anfangswertproblem (1), (2) enthalt u.a. die folgenden Sonderfalle: 


a) Die gewohnliche explizite Differentialgleichung erster Ordnung 
u'(t)=g(t,u) mit u“(0)=y, (3) 
ferner die aus (3) folgende ,,integrierte Form“ 
t 
u() =n-+fe(s,u(s)) ds, 
sowie mit z(t) =w'(t) die Gestalt 
t 
2() =e(tn +S 2(6) ds). 
0 
b) Die Fredholmsche (nichtlineare) Integralgleichung zweiter Art 
7 
u(t) =f K(s,t, u(s))ds+g(t), 
0 


in der fiir spezielle Kerne K auch die Volterrasche Integralgleichung enthalten ist. 
c) Die aus der Anfangswertaufgabe fiir die explizite Differentialgleichung 


zweiter Ordnung 
u'(t)=g(t,u,w'), ulO)=n, w(0)=C 
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folgenden , integrierten Formen ‘‘ 


wl) =C+Se(o.m(sw's)) ds, uO)=n 
und (mit z(¢)=w’(d)) | 
2'() =e (tin +f 209) 4s, 200), ZO) as Ci. 


Wie man aus diesen Beispielen sieht, kann die Anfangsbedingung (2) unter Um- 
standen in der Gleichung (1) mitenthalten sein. 

Die der nachfolgenden Untersuchung zugrundeliegende Abschatzungsmethode 
ist im ersten Abschnitt I enthalten. Vereinfacht dargestellt lautet sie: Ist Fu 
ein ,monoton wachsender“ Operator und geniigt die Funktion f bestimmten 
Monotoniebedingungen, dann ist die Aufgabe (1), (2) ,,von monotoner Art” (vgl. 
L. Cotratz [1]), d.h. es gilt das folgende Lemma: 


Aus 
v(0) < u(0) =7 < w(0) 


Ue f(t,v',v, Fv) <flt,w,u, Fu) =0< f(t, wv’, wv, Fw) 


ees v(t)<u(t)< w(t) auf ganz OStST. 

Lassen sich also zwei solche Funktionen v und w mit den genannten Eigen- 
schaften angeben, dann sind v und w untere und obere Schranken fir jede 
Lésung w des Anfangswertproblems (1), (2). Fiir gewisse Volterrasche Integral- 
und Integro-Differentialgleichungen sind die Voraussetzungen dieses Abschnitts 
erfiillt, die Abschatzung ihrer Lésungen wird also durch dieses Lemma unmittel- 
bar geleistet. 

Der nachste Abschnitt II befaBt sich mit der Fehlerabschatzung einer Nahe- 
rungslésung v(t) von (1), (2), fiir die allein der ,,Anfangsfehler“ e=v(0) —y und 
der ,,Defekt“ d(t)=/(t, v’, v, Fv) bekannt ist. Gesucht ist eine Abschatzung fiir 
den Fehler v(¢) —w(é) allein aus der Kenntnis von e und d(t) (,,Defektabschat- 
zung''). Im allgemeinen hat der Defekt d(¢) nicht durchweg gleiches Vorzeichen 
und auch Fw und f geniigen fiir gewohnlich nicht den Monotonievoraussetzungen 
des Abschnitts I, so daB sich das Lemma nicht direkt anwenden laBt. Gibt es 
nun zu (1), (2) ein ,, Vergleichsproblem‘ 


o(t,e',e,2e)>d(t), e(0)>e, 


dessen Operator 2 und dessen Funktion w den Voraussetzungen des Lemmas 


geniigen, und sind w und Q in einem gewissen Sinne Majoranten fiir f und F, 
dann gilt 
lv¢) —u@®|<o@ in OStST, 


d.h. eine Lésung @ (¢) dieses Vergleichsproblems ist eine Schranke fiir den absoluten 
Fehler der Naherungsfunktion »v (¢). 

Lat sich nun zeigen, daB die von e und d(t) abhangige Funktion @ (é) fiir 
e>0 und d->0 ebenfalls verschwindet, so heiBt das, daB der Abstand v—u 
zweler Losungen v und w von (1), (2) beliebig klein wird, d.h. daB das Anfangs- 


wertproblem (1), (2) hdchstens eine einzige Lésung besitzt. Der letzte Abschnitt III 
gibt solche Einzigkeitskriterien an. 
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Die erwahnten Ergebnisse lassen sich sofort auf Operatoren-Differentialglei- 
chungen hdherer Ordnung sowie auf Systeme davon iibertragen, der Einfachheit 
halber soll jedoch von dieser Verallgemeinerung abgesehen werden. Eine Uber- 
tragung auf Funktionen mehrerer Veranderlicher ist ebenfalls méglich; eineTheorie 
fiir hyperbolische Differentialgleichungen — die im wesentlichen auf denselben 
Grundideen basiert wie die folgenden Betrachtungen — wurde vor kurzem von 
W. WALTER in dieser Zeitschrift verdffentlicht [7]. Bei der Ausarbeitung der 
nachfolgenden Ergebnisse war mir der Rat von Herrn Doz. Dr. W. WALTER 
sehr wertvoll. 


I. Ober- und Unterfunktionen bei Operatoren-Differentialgleichungen 

mit monotonen Operatoren 

Bezeichnungen und Definitionen. Es sei T>0 eine Konstante. Das abge- 
schlossene Intervall 0<¢<T der reellen Veranderlichen ¢ sei mit J, die links 
offene Strecke 0<¢XT mit J, bezeichnet. Samtliche im folgenden eingefiihrten 
Funktionen seien als reellwertig und eindeutig erklart vorausgesetzt. 

Der Operator F sei auf dem Raum C(J/) der auf J stetigen Funktionen « (¢) 
erklart, genauer ist w€C(Z) und v=F yu, so sei v (¢) definiert fiir alle €J,. Die Ge- 
samtheit dieser Operatoren F werde mit 7 bezeichnet. Die Teilmenge 3,3 [3_CB| 
der (schwach) monoton wachsenden [fallenden] Operatoren sei erklart durch die 
Eigenschaft: Sind wu, v€C(Z) und ist fiir eine beliebige Stelle s CJ, stets 


u(t)<v(t) in 0<t<s, (4) 


dann soll sogar Fu<[=]Fv fir t=s (5) 


gelten. Ist FER, , so ist —FE_. Fir die nachfolgenden Untersuchungen geniigt 
es daher, allein den Fall F€%, zu betrachten. 

Beispiel 1. Die folgenden Operatoren gehéren zu %,, jedes Beispiel enthalt 
eine umfassendere Funktionenklasse. Die verwendeten Integrale sind im Sinne 


einer der gebraéuchlichen (auch uneigentlichen) Integraldefinitionen aufzufassen. 
t 


a) Fu= fu(s)ds, 


0 


15 
b) Fu =f K(s,t)u(s)ds, 
0 
wenn K(s, ¢) =0 ist fiir s, € Jy) und wenn K(s, #) fiir alle ¢€ J) integrierbar tiber s ist. 
t 
c) Fu=f K(s,t,u(s))ds, 
0 


wenn K(s,¢, z) — K(s,¢,2) =0 bleibt fiir s, £€ J) und alle z= z und wenn (s; t,u (s)) 
fiir ¢¢ 7, und alle Funktionen w€ C() integrierbar iiber s ist. 

Definition. Die Menge Z =Z (I) der auf J zulassigen Funktionen sei die Gesamt- 
heit aller auf J stetigen, in J, differenzierbaren Funktionen 1 (i). 

Problemstellung. Die Funktion f (¢, x, y, 2) sei definiert fir t€l), —w<%, y, 2 
<-+ oo. Man wahlt einen Operator F€% und eine reelle Zahl 7. Ist die Funktion 
ne Z (I) und geniigt w(t) der Gleichung 

fi,w,u,fu) =O ind, (1) 
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und der Anfangsbedingung u(0) =n, (2) 


dann hei&t u Lésung des Anfangswertproblems (1), (2). 
Definition. Ist die Funktion f(t, x, y, 2) definiert fiir t€I,, —w<%, Y, 2< + oe 
und gilt dort sogar 
f(t, x,y, 2) —f(4,%, y,2) 20 sabe ABE? (ESTA (6) 
d.h. ist f dort (schwach) monoton wachsend in x und (schwach) monoton fallend 
in z, so soll dafiir {¢ M geschrieben werden. 


Satz 1. Es sei fC M, fiir die Funktionen v, wEZ gelte 
v(0) < w(0) (7) 
und mit dem Operator FE%, set noch 


f(é,v',v, Fv <ftw',w,Fo) my. (8) 
Dann 1st sogar 
v(t)< w(t) auf ganz I. (9) 


Beweis durch Antithese. Angenommen, die Behauptung (9) sei falsch, dann 
gibt es Stellen ¢/€I mit v(¢’)=w(¢’). Man bildet t=inf?#’. Wegen der Stetigkeit 
von v und w gilt dann 


v(t) = w(é), (10) 


nach (7) ist also insbesondere ¢>0. Wegen v(t)<w(é) in O<t<t und wegen 
Fes, ist dann 


FusFw fir t¢=t. (11) 
Ebenfalls gilt damit 
eee cal mee fiir allefaus O<t<é: 


fiir t+ folgt daraus wegen der Differenzierbarkeit von v und w auch 


v'() 2 w'(). (12) 
Wegen f€M findet man damit fiir die Stelle t—=t aus (10) bis (12) die Ungleichung 
f(t, v',v, Fv) =f(t, w', w, Fw) im Gegensatz zu (8), und aus diesem Widerspruch 
folgt die Behauptung. 
Aus dem Beweis sieht man sofort den 


Zusaiz. Hangt die Funktion /(é, x, y,z) von x nicht ab, d.h. entartet die 
Operatoren-Differentialgleichung (1) zu einer reinen Operatorengleichung, so 
genligt es fiir die Giiltigkeit von Satz 1, wenn die Funktionen v, w€C(ZL) sind. 

Man bendtigt diesen Zusatz z.B. fiir das folgende 


Beispiel 2. (Lemma von GRonWALL, vgl. G. SANSONE [5], $.30) Es seien 
K=20, AZO, B und y>0 vier Konstanten. Die Funktion veEC(I) geniige auf 
I der Integralungleichung 


t 


v(t) — J (Kv) + 4) ds—Bso. (13) 
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Um Satz 1 anwenden zu kénnen, wird gesetzt: 
t 


f(t, *,y,2)=y—z2—B, n=B, Fu =f (Ku(s) +A) ds; 
0 
es ist f/€M und FeEs,. Fiir die auf J zulassige Funktion 
w(t) = (At+ B)expKt+yexp2Kt 
gilt offenbar 
y Piney Koes 0 


t, , ple Tae 
f(t, x, w, Fw) A (exp Kt —1—K2) 4 + ({itexp2K?) fir K>0, 


K 


also f(t, x, w, Fw)>0 fiir ¢€ J, und beliebige Werte x; weiter ist w(0)=B+y> 
B=u(0). Also ist nach Satz1 stets u(t)<v(t) auf J. Dies gilt fiir alle Werte 
y> 0, in der Grenze fiir y +0 ist also 


v(t) S (At+ ByexpKt aufl 


fiir jede Funktion v€C(Z), die der Integralungleichung (13) geniigt. 

Wird in Satz 1 in der Bedingung (7) das Gleichheitszeichen zugelassen, so 
ist die Folgerung (9) i.a. nicht mehr richtig, wie schon einfache Beispiele zeigen. 
Um die Forderung (7) trotzdem abmildern zu kénnen, fiihrt man eine neue 
Bezeichnung ein gemaB der folgenden 


Definition. Sind zwei Funktionen v(t) und w(t) in J, erklart und existiert 
eine Nullfolge s, (v=1, 2,...), fiir die s,EZ,, lims,=0O und v(s,)<w(s,) gilt fiir 
alle »=1, 2,..., so soll dafiir abkiirzend on 


v(0+) < w(0+) (bzw. w(0+) >v(0+)) 
geschrieben werden. 

Beispiel 3. Sind v, we C(Z) und ist v(0)<w (0), so gilt offenbar auch v (0+) < 
w(0+). Sind v, w in einer Umgebung von ¢=0 differenzierbar und ist v (0) =w (0), 
aber v’(0)<w’(0), so ist ebenfalls v(0+)<w(0+). Weitere hinreichende Bedin- 
gungen lassen sich leicht angeben. 

Mit dieser Definition laBt sich die Voraussetzung (7) noch abschwachen, man 
kann v (0) =w (0) zulassen und allein v(0+)<w(0-+) fordern. Da man den Satz 1 
auf jedes der Intervalle s,<¢<7T anwenden kann, gilt (9) noch in ganz Jy. Auch 
die Bedingung (8) kann noch gemildert werden. Wie man aus dem Beweis sieht, 
geniigt es, sie allein fiir die dort benutzte Stelle ¢ vorauszusetzen. Damit gilt der 


Satz 2. Es sei fCM, fiir die Funktionen v, wEZ* gelte 


»(0-+) <w(0+) (7a) 
und mit dem Operator FE, set noch 
f(t, v', v, Fv) < f(t, wv’, w, Fw) (8a) 


* Offenbar wiirde es mit der Bedingung (7a) geniigen, wenn v, w als in J, stetig 
und differenzierbar vorausgesetzt wiirden. Die Stetigkeit der zuldssigen Funktionen 
an der Stelle ‘=0 konnte also aufgegeben werden. Da diese Bemerkung fiir die 
praktischen Anwendungen nicht sehr wichtig ist, und bei einem Teil der folgenden 
Satze auf die Stetigkeit der zulassigen Funktionen an der Stelle =0 nicht verzichtet 
werden kann, wurde auf die Mitnahme dieser Verallgemeinerung verzichtet. 
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fiir alle Stellen i€I,, fir die v (f) =w (2) und v(t)<w(t) fiir 0<t<t gilt. Dann 
BSUS v(t) <w(t) in ganz Ip. (9a) 

Bemerkung. Kennt man die in den Satzen 1 und 2 auftretenden Funktionen 
v und w explizit, so ist die Aussage (9) baw. (9a) i.a. uninteressant, weil man 
mit der Kenntnis von v und w natiirlich auch weiB (oder wenigstens stets prinzipiell 
nachpriifen kann), ob v<w ist oder nicht. Ihre Bedeutung erhalten die beiden 
Satze erst dadurch, daB man sie auf die i.a. noch unbekannten Lésungen des 
Anfangswertproblems (1), (2)* anwenden kann. In diesem Falle ist — analog 
zu einer Definition von O. PERRON [3] bei gewohnlichen Differentialgleichungen — 
die folgende Setzung zweckmabig : 

Definition. Eine in I zulassige Funktion v (¢) werde als Oberfunktion (beziiglich 
des Anfangswertproblems (1), (2)) bezeichnet, wenn 


v(O+) > u(0+) (14) 
ist fiir jede Losung u(t) von (1), (2) und wenn 
[0s Lo) O tee (15) 


gilt. Sind dagegen (14) und (15) beide mit dem <-Zeichen erfiillt, so soll v eine 
Unterfunktion heifen. 


Nach Satz 2 gilt dann das 


Lemma 1. Fiir jede beliebige Losung u(t) von (A), (2), eene Unterfunktion v (t) 
und eine Oberfunktion w(t) zu (1), (2) gilt stets 


vi)<ul)<w(t) my, 
falls fEM und FES, ist. 

Bemerkung. Dieses Lemma bildet die Grundlage fiir die folgenden Abschat- 
zungs- und Eindeutigkeitssatze; in der von Herrn CoLLAtz eingefiihrten Sprech- 
weise ist dann das Anfangswertproblem (1), (2) ,,von monotoner Art’. Zur 
Abschatzung einer Losung geniigt es danach, passende Ober- und Unterfunktionen 
anzugeben. Lat sich sogar je eine Ober- und eine Unterfunktion mit beliebig 
kleinem Abstand finden, dann hat die Aufgabe (1), (2) héchstens eine einzige 
Lésung. 


Beispiel 4 (Plattengrenzschicht). Zur Untersuchung der Singularitat der 
Blasiusschen Differentialgleichung der Plattengrenzschicht betrachtet man (vgl. 
B. Punnis [4]) das Anfangswertproblem 


git sey = 0, 15 (ON 4" (O) a0 Orie (16) 


In 0St<oo. Setzt man u(t)=y’'(t), dann geniigt w offenbar der mit (16) gleich- 
wertigen nichtlinearen Volterraschen Integro-Differentialgleichung 
t 


f(t, uv’, u, Fu) =u' — exp ft = s)u(s\ds=0 (17) 


4 (0) =0. (18) 


* Wobei die Gleichung (1) auch durch eine Ungleichung ersetzt werden darf, 
vel. Beispiel 2. 


unter der Anfangsbedingung 
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t 
| Setzt man Fu =exp f (t—s) u(s) ds, so ist FC%, und f(t, x, y, 2) =x —z geniigt 
0 


——— ee oe 


der Bedingung (6), d.h. es ist f€M. Wie man leicht nachpriift, hat eine Lésung 
von (17), (18) bei £=0 die Entwicklung u(t) =¢-+ ¢4/24-+4 --.. Eine einfache Ober- 
funktion ist 

W(t) =3(a —t)®9—3a2% mit «=y/6 


in dem Intervall O<t<«, denn es ist w(t) =t+3#2/2a+---, also u (0+)<w(0+), 
und 
f(¢, w', w, Fw) = 6(« — t)3[1 — exp —at(204+2/4]>0 in 0<t<a. 
Eine Unterfunktion ist 
t fir OSts@ 
v (t) ae __ f\-2 1 oes 
B( Qo)" So oltin, wt 20 
wegen v(0+)< (0+) und 
Asta salsa 
‘ACEC ae tap oo 
6 (2% — ¢)-$ [4 — exp(1+ 4a (¢ — «)?)] 
AN = ile Or DOG 


Eine singulare Stelle jeder Lésung w(f) von (17), (18) liegt damit zwischen 
V6 und 276, bessere Schranken lassen sich leicht angeben. 


Bemerkung. Es ware offenbar nicht nétig gewesen, zur Abschatzung der 
Lésungen von (17) die vorstehende Theorie zu benutzen, da die Funktion 


line ds 


0 
zu (17), (18) ist, wenn man die bekannte Theorie der Unter- und Oberfunktionen 
fiir Anfangswertaufgaben bei gewoéhnlichen Differentialgleichungen heranzieht 
(vgl. etwa M. MULrER [2]). Danach gilt z.B. fiir jede zweimal differenzierbare 
Lésung y(¢) der Anfangswertaufgabe 


Dy=y''(t) —glt,*%,y)=0, y(0)=n,- y'(0) =e (19) 
a()<yt<2@) aut J, (20) 


t 
fos) ds| Unter- [Ober-] Funktion 
0 


stets 


wenn x(t) und z(é) auf J zweimal differenzierbar sind, wenn gilt: 
HAO) yes 2 (0), (0) So 240), 
DAAAO 2 De ind, 


und wenn g(t, 7, s) in der Veranderlichen s schwach monoton wachst. 

Besitzt g diese Monotonie-Eigenschaft nicht, so ist die Abschatzung (20) i.a. 
nicht richtig, auch wenn (21) erfiillt ist. Es kann nun sein, daB zwar git, 7, s) 
in s nicht monoton wachst, daB also eine Lésung y(¢) von (19) sich necht durch 
Unter- und Oberfunktionen abschatzen laBt, daB dagegen g(t, 7, s) ein monotones 
Wachstum beziiglich 7 zeigt. Dann gestattet aber die mit (19) identische Integro- 
Differentialgleichung (w(t) =y’(t)) 


w(t) =g(tn + fu (s)ds,u(),  u(0) =E 


166 Kar NICKEL: 


Abschitzungen durch Unter- und Oberfunktionen » (/) und w (¢) nach dem Lemma 1. 
Damit ist dann natiirlich auch 


heat | +) 79) ds 


in Schranken 


x(t)=n Wiad und z(t)=7 rae 


eingeschlossen, obwohl sich diese Abschatzung auf dem direkten Wege nicht 
zeigen lieB. Wie dieses Beispiel zeigt, kann es also oft fiir eine Abschatzung 
zweckmabig sein, von einer gewohnlichen Differentialgleichung zu einer Integro- 
Differentialgleichung iiberzugehen. 


Gegenbeispiel 1. Zu einem Anfangswertproblem (1), (2) braucht es il.a. weder 
Unterfunktionen noch Oberfunktionen zu geben, selbst wenn f/€M und Fes, 
ist. Ist z.B. {CM und F€§%,, hangt f(é, x, y, z) von x nicht ab und ist f zusatzlich 
noch schwach monoton fallend in y, dann gibt es sicher keine Ober- oder Unter- 
funktionen zu (1), (2). Denn ware w eine Oberfunktion, dann galte nach dem 
Lemma 1 stets w>wu in I,, also auch Fw=Fu in I, wegen FE%,. Wegen der 
vorausgesetzten Monotonie von f ware dann auch 0=/(t, wv’, u, Fu) =f (t, w’, w, Fw) 
in J, in Widerspruch zur Bedingung (15) der Oberfunktionen. Entsprechend 
schlie8t man fiir Unterfunktionen. Das Lemma 1 und die daraus zu ziehenden 
Folgerungen lassen sich also auf diese Klasse von Operatorengleichungen nicht 
anwenden. Ein Beispiel dafiir ist die Klasse der Volterraschen Integralgleichungen 
erster Art der Gestalt P 


LNs: t, u(s))ds =gi(d). 


Es ist fiir die Anwendungen oft unhandlich, da8 in den Ungleichungen (7) 
und (8) des Satzes 1 (bzw. (7a) und (8a) bei Satz 2) das Gleichheitszeichen nicht 
zugelassen ist. Man ist dadurch gezwungen, bei der Definition der Unter- und 
Oberfunktionen in den Bedingungen (14) und (15) das Gleichheitszeichen eben- 
falls auszuschlieBen. Da diese strengen Forderungen in der Tat notwendig sind, 
da8 also etwa Satz 1 i.a. nicht mehr richtig ist, wenn das <-Zeichen in (7) und 
(8) durch < ersetzt wird — selbst dann nicht, wenn auch in (9) das <-Zeichen 
zugelassen wird — zeigt das folgende 


Gegenbeispiel 2. Es sei T=2, 
ee t 
f(t, x, ¥,2)=%—4 Vy2—4(t—1)? z, Fu= fu(s)ds; 
0 
offenbar ist f€M und FC%,. Man betrachtet also die Integro-Differential- 
gleichung 


t 
ftw’, u, Fu) =u —4Vu2— 4(¢— 1)? f u(s)ds=0. 
0 

Fir die beiden zuldssigen Funktionen v(¢)=(f—1) und w (4) =0 gilt v(0)= 
—1<0=w (0) und f(t/e 0, Fo) == (¢—1)®<0=f(t, vw’, w, Fw) in 0St<2, 
dabei tritt das Gleichheitszeichen allein an der Stelle t=1 auf. Trotzdem ist 
v(2)=+1>0=w(2), d.h. die Folgerung (9) des Satzes 4 ist falsch, obwohl alle 
Voraussetzungen von Satz 4 erfiillt sind bis auf die Bedingung (8), die in J sogar 
nur an der einzigen Stelle t=1 verletzt ist. 
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Dagegen lassen sich die Bedingungen (7) und (8) von Satz 1 bzw. (7a) und 
(8a) von Satz 2 noch abschwachen, wenn man die iibrigen Voraussetzungen ver- 
scharft. Eine besonders einfache Méglichkeit ist die folgende, bei der die schwache 
Monotonie bei der Definition von §, und M durch eine starke Monotonie ersetzt 
wird*. Dazu dienen die folgenden 


Definitionen. Die Teilmenge eas der stark monoton wachsenden ** 


Operatoren F sei erklart durch die Eigenschaft: Sind v, w€C(Z) und ist fiir eine 


beliebige Stelle s€J, stet 
eee Pees PASS ho See ey Ey MRO PSS, a 


d ll 
ann soll sogar Fu<Fv fir t=s (5a) 


gelten. Fiir die drei Operatoren aus Beispiel 1 gilt offenbar: Zu a): FES, , zu b): 
FEe%,, wenn sogar K>0 ist, zuc): FES,, wenn K(s, ¢, z) —K(s, t,Z)>0 bleibt 
fiir s,¢€J, und z>z. 
Ist f(¢, x, y, z) definiert fiir ¢€[), —coo<x, y,z<-+ oo und gilt dort 
TVA PUG Te) SOR TE aA, <2 2} 
d.h. ist f dort stark monoton wachsend in x und stark monoton fallend in z, 


so soll f€M geschrieben werden. 


Damit gilt der folgende Satz 3, in dem gegeniiber Satz 2 in den Ungleichungen 
(8) und (9) noch das Gleichheitszeichen zugelassen wurde, wahrend dafiir sogar 
fE€M und FER, gefordert wird: 


Satz 3. Es set fe M, fiir die Funktionen v, w EZ gelte 


v(0+) < w(0+) (7a) 
und mit dem Operator F EF, set noch 
{(t,v',v, Fv) S f(t, wv’, w, Fw) (8b) 


fiir alle Stellen t€I,, fiir die v(t)=w(t) und v(t)<w(t) fir O<t<t gilt. Dann 
cee v(t) <w(t) auf ganz I. (9b) 
Der Bewets verlauft ganz entsprechend zu dem von Satz 2. 

Bei dem Gegenbeispiel 2 war die starke Monotonie von f(f, x, y, z) hinsicht- 
lich z gerade an der ,,kritischen“ Stelle i—1 verletzt, simtliche anderen Voraus- 
setzungen zu Satz 3 waren erfiillt. 

Nach Satz 3 liegt es nahe, auch den Begriff der Unterfunktion und der Ober- 
funktion zu erweitern zu der 

Definition. Eine in I zulassige Funktion v(¢) werde als Oberfunktion im 
weiteren Sinne (beziiglich des Anfangswertproblems (1), (2)) bezeichnet, wenn 


v(0+) < u(0+) (14) 


* Es soll jedoch ausdriicklich darauf hingewiesen werden, daB sich auch noch 
andere hinreichende Bedingungen fiir diesen Sachverhalt angeben lassen, die im Falle 
einer reinen Differential- oder einer reinen Integralgleichung sogar ahnlich einfach 


werden. 
xx Auf die analog erklarten monoton fallenden Operatoren kann man wie oben 


fiir das folgende verzichten. 
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ist fiir jede Losung w(t) von (1), (2) und wenn 
T(b,0.0,1 0) 2 0 mk, (15a) 

gilt. Sind dagegen (14) und (15a) beide mit dem <- bzw. <-Zeichen erfiillt, 
so soll v eine Unterfunktion im weiteren Sinne heiBen. 

Nach Satz 3 gilt dann das 

Lemma 2. Fiir jede beliebige Lésung u(t) von (1), (2), eime Unterfunktion im 
weiteren Sinne v(t) und eine Oberfunktion im weiteren Sinne w (t) zu (1), (2) gut stets 

v(t) <u(t)< w(t) auf ganz I, 

falls {EM und Fe x, ast. 


II. Fehlerabschatzungen 


Ist v(t) eine Naherungslésung zum Anfangswertproblem (1), (2), so ist der 
,Anfangsfehler“’ e=v(0)—y#O und ebenfalls gilt fiir den ,,Defekt" d()= 
}(t, v', v, Fv) ~0 in Jy. Im folgenden soll allein aus der Kenntnis von Schranken 
fiir Anfangsfehler e und Defekt d(¢) eine Abschatzung fiir den Fehler v (f) — uw (¢) 
der Naherungslésung v(t) angegeben werden. 


Satz 4. Die Funktion vCZ habe die Eigenschaften 
|2(0) —n| Se, (22) 
\fG,v',v,Fv)|=d6@ mdy. (23) 


Die Funktion 6 (t) sea definiert in I,, f(t, x, y, 2) set erklart fir t€Iy, —w< %, y, 4 
<+oo, es ser FC. Werterhin gelte mit einer Funktion om©M[M] und einem 
Operator QE®,[G,] die Abschatzung 


f(t, x, wv, Fw) — f(t, %,0, Fw) Sa (t, x —%,w —W, Q(w — @)) (24) 


fiir beliebige Zahlen x,x und alle Funktionen w,@CC(I) mit der Eigenschaft 
ww auf I. 


Gentigt dann die Funktion 0€Z der Bedingung 


o(t,e',0,20)>[2]d@) inI, mit o(0)>e, (25) 
dann gilt fiir jede Losung u(t) des Anfangswertproblems (1), (2) die Abschdtzung 
Jo) —u@|<[S]oW auf 1, (26) 


fi : ne ast eine Schranke fiir den absoluten Fehler der Ndaherungslésung v (t) zu 
i), (2) 

Bemerkung. Wahrend bei den Abschatzungen des vorigen Abschnitts durch 
Unter- und Oberfunktionen nur Funktionen f/€M[M] und Operatoren FE, (3. 
zugelassen werden konnten, brauchen bei der Anwendung des Satzes 4 keine 
solchen Monotonie-Voraussetzungen getroffen zu werden, es gentigt vielmehr, 
wenn allein die ,,Majoranten‘‘ w und Q den entsprechenden Bedingungen gentigen. 


Damit ist der Anwendungsbereich von Satz 4 gegeniiber den Satzen von Ab- 
schnitt I stark vergr6Bert worden. 
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Beweis. Es soll zunachst v—u<[S]o auf I nachgewiesen werden. Man 
wendet dazu Satz 1 [Satz 3] auf den Ausdruck m mit dem Operator Q (statt f 
mit F) und die beiden Funktionen v—w und 9 (statt v und w) an. Nach (22) 
und (25) ist v(0)—w(0)Se< (0), wahrend wegen (23), (24) und (25) stets 


w(t, v' — wv’, d — u, Ov — u)) Sfi(t,v',v, Fv) — f(t,w', uv, Fu) 
< 6(4)<a(t, 0’, 0, Qo) 


‘in J, ist. Nach Satz 1 [Satz 3] bleibt dann auch v(t) —u(t)<[S]o (é) auf ganz J, 
denn w und Q geniigen den dort verlangten Monotoniebedingungen. Ganz ent- 
sprechend beweist man auch u(t) —v(t)<[S<]o(é), womit (26) vollstandig be- 
wiesen ist. 


Beispiel 5 (Lipschitz-Bedingungen). Im folgenden sei die Funktion L(t) in J, 
stetig, das Integral 


t 
f L(s) ds existiere (eventuell uneigentlich) 
0 


fiir ¢E I. 
a) Man betrachtet das Anfangswertproblem der gewoéhnlichen Differential- 
gleichung erster Ordnung 


u' — g(t, u)=0 u(0)=y, (27) 


verwendet also etwa f(t, x, y, z)=x—f(t, vy), F beliebig. Die Funktion g(é, y) 
sei dabei definiert fiir #€ I), — oo< y<_+ ov; die reelle Zahl 7 sei beliebig. Geniigt 
g der Lipschitz-Bedingung 


g(ty)—sby) SLO (y—y) fir 7€h, yv2y, (28) 


dann ist (24) (fiir einen beliebigen Operator 2€ 3,) befriedigt mit der Funktion 
w(t, x, y, 2) =x—Li(i)y, dabei ist @€ M. Die Ungleichungen (25) gehen damit 
tiber in die Differentialungleichung 0’ —L(f)o>dé(é) mit @(0)>e. Fir y>0 
und 6€C(Z) ist 
e(t) =y(1+4 +o(¢) 
mit 


o(t) =€Exp i fheas e+ f 8(s) exp Seine dtds (29) 


offenbar eine Lésung, fiir yO folgt daraus die bekannte Fehlerabschatzung 
|v(t) —u(f)| So) mit a(t) aus (29) (30) 


fiir alle Lésungen w(é) von (27) und alle Naherungsfunktionen v(¢), fir die 
|v (0) —u(0)|Se und |v’ —g(é, v)| Sd(4) in Jy gilt. 

Es soll dabei besonders darauf hingewiesen werden, daB es zur berdsettigen 
Abschitzung der Differenz v — wu nach (29), (30) gentigt, eine ezmseztege Lipschitz- 
Abschatzung (28) vorauszusetzen. Damit sind auch negative Funktionen L(¢) 
zugelassen (vgl. etwa M. MULLER [2}), in diesem Falle ist das Fehlerverhalten 
besonders ,,gutartig‘‘. In der Literatur wird statt (28) oft 


lf») -f69|SLOly—F| fiir t€L, —w<y,~<t+o 
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gefordert, wodurch fiir L <0 die Abschatzung (30) unter Umstanden ganz erheb- 


lich verschlechtert wird. 
b) Man betrachtet die (nichtlineare) Volterrasche Integralgleichung zweiter 


Art 
> J K(s,t, w(s))ds=g(d), (31) 


die in (2) vorgeschriebene Anfangsbedingung w(0)=7 ist mit g(0)=7 in (31) 
implizit schon enthalten. Man verwendet also etwa die Funktion f(t, x, y, z)= 
y —z—g(t) und den Operator 
t 
Fu=f K(s,t,u(s)) ds. 
0 
Es ist F€%, wenn K integrierbar ist fiir s,t€ J und we C(I). Gentigt der Kern 
K(s, t, 2) fir s, ¢€I der Lipschitz-Bedingung 
K(s, t, 2) — K(s, t, 2) S L(s) N@ (z—2) firalle 2272, (32) 
mit einer in J stetigen Funktion N(é), dann ist die Ungleichung (24) befriedigt 
mit der Funktion w (t, x, y, z) =y — N(¢)z und dem Operator 


Qu= jepe u(s) ds. 


Ist L=0O und N&O auf J, dann gilt m¢€ M und 2€%,, also kann Satz 4 ange- 
wandt werden. Die Ungleichung (25) geht iiber in die Volterrasche Integral- 
ungleichung zweiter Art 


S16) s)ds>6d(f) mit o(0)>e, (33) 
wo ¢ und d(¢) durch 
|v(0) —u(0)| Se, 


HO PG UIGN = s|<0@ ind, Oe 
definiert sind. , 
Fir y>0, d€ C(Z) und N(0) +0 ist 
et) =y Ni) (exp J E69) N (s) ds + a(t) 
mit 
o(t) =6() +N) exp f L(9 N(s) ds {[Max (c, 5(0)) — 6(0)]/N(o) + 
(35) 
+ Jo6) () exp — { L() Nip) dvds} 
eine Lésung von (33). Mit y—>0 folgt daraus die Fehlerabschatzung 
|2) — "| <a) mit o( aus (35) (36) 


fiir alle Lésungen w(t) von (31) und alle Naherungslésungen v(é), die den Ab- 
schatzungen (34) gentigen. 
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Fiihrt man die Differentialgleichung (27) durch Integration iiber in die 
Volterrasche Integralgleichung zweiter Art 


u (t) =i ets; (s)) ds=7 


entsprechend zu (31), so ist hier nur der Fall L=0 zulassig, d.h. der oben erwahnte 
Vorteil bei den Fehlerabschatzungen geht bei dieser Transformation verloren! 


c) Man betrachtet das Anfangswertproblem mit der Fredholmschen Integro- 


Differentialgleichung 
vis 


u' (t) — g(t, u(Z)) a t,u(s))ds=0, u(0) =n, (37) 
wahlt also 


ip 
1%, 4,2) =x—el,yy—2 und Fu= fiK(sit, u(s))\ds. 
0 


Ist K integrierbar fiir s,t€J, wEC(Z), so gilt FE %. Es sei 
g(t, y) —gy)S2Ly—y) fir tel, undalle y2y 
und fiir ¢€T gelte 


J [K(s,#, w(s)) — K(s,t,@(s))]ds< nf [w(s) —w(s)] ds (38) 


0 


fiir alle Funktionen w, w€C(Z) mit w(t)=w(t) auf J. Dann wird (24) befriedigt 
t 


durch die Funktion o (é, x, y, z)=x—2Ly—Nz, wenn Qu= f u(s)ds gesetzt 
0 
wird. Es ist QE, und w€M, wenn N=O gilt. Ist d(f)<6, dann geht (25) 
tiber in die Ungleichungen 
t 
eo’ —2Le—Nfoe(s)ds>6, e(0)>e, 
0 
die fiir alle y>0O durch 
e (f) =yexp(A+ 1) é+ o (4) 


mit 
o- fir—L? 4+ N=0 39) 
~ [lea +6) exp At — (eu +6) exppé]/(A—m) fir L2+N+40 
und 
A=L+VE24+N, w=L—VDR+N 
befriedigt werden. Fir yO findet man dann die Fehlerabschatzung 
|v(t) —uw(é)| So(¢) auf I mit o aus (39) (40) 


fiir alle Lésungen w(t) von (37) und alle Naherungslésungen v (¢), deren Anfangs- 

fehler durch ¢ und deren Defekt durch 6 beschrankt wird. Wegen 420 nimmt 

die Schrankenfunktion o(¢) nur fiir L=N=6=0 oder fiir e=d=0 nicht zu. 
Fiir beliebige Funktionen / und Operatoren /€% sind die Voraussetzungen 

des Satzes 4 im allgemeinen nicht erfiillbar, insbesondere braucht es keine 
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Funktion w und keinen Operator 2 zu geben, so daB (24) befriedigt ist. Das 
zeigt das folgende 

Gegenbeispiel 3. Es sei f(t, x, y, z) =y*, dann geht die Ungleichung (24) iiber 
in 

w? — 2 = (w —@) (20 + w —B) Zol,...,w—@B,...) 
fiir alle Funktionen w, @C C(J) mit w=. Ware @ eine festgehaltene Funktion, 
so kénnte etwa w(t, x, y, 2) =2wy+y? gesetzt werden; fiir beliebige Funktionen 
w€C(I) ist diese Ungleichung jedoch offenbar nicht zu betfriedigen. 

Im folgenden Abschatzungssatz 5 sind daher die in Satz 4 noch bendtigten 
Voraussetzungen sehr stark reduziert. Dabei werden besonders die speziellen 
Eigenschaften der Naherungslésung v(t) noch weiter ausgentitzt; auBerdem zeigt 
es sich, daB man durch einseitige Abschatzungen oftmals zu giinstigeren Fehler- 
schranken gelangt. Es gilt der 

Satz 5. Die Funktion f(t, x, y, 2) set in t€Ig, —w0< x,y, 2<+00 erklant, 
der Operator F sei €§. Die Funktionen u,v€Z(L) seien eine Lésung und eine 
Naherungslésung zu dem Anfangswertproblem (1), (2) und zwar médgen mut geerg- 


neten Konstanten e, € und in I, erklarten Funktionen 6(t), d(t) die Abschatzungen 
—€é<v(0) —u(0) Se, (41) 
f(t,u',u,Fu)=0, —d()<f(t,v',v,Fv) S60 ml, (42) 

gelten. 
Die Funktionen w,0C M [M], die Oferatoren Q, QE 19.) und die Funk- 


tionen 0,0 Z(L) mégen den folgenden beiden Systemen aus je drer Unglerchungen 
gentigen, es gelte namlich 


i(é,v', v, Fv) — f(t,’ @',v—@, Flv —@)) 200, 0,20) (43) 

fiir alle t€I, mit v=u-+o, falls fir 0Xt<t stets v—u<o bleibt, 
o(t,e',e,2e)>[2)6 m Ih, e(0)>e; (44) 
f(,0 +09 +0,F(v+0) —f(,v',v, Fv) S000, 20) (45) 

fiir alle t€Iy mit v=u—6, falls fiir O0X<t<t stets u—v<o Dleibt, 
@(t,0’,0,20)>[=]d@) inl, O(0)>8. (46) 

Dann gilt die Abschiétzung 

—e@@[S}<v) —v¥<[S]e (47) 


auf ganz I. 


Der Bewets verlauft mit Satz 2 [Satz 3] ganz entsprechend zum Beweise des 

Satzes 4. Satz 4 ist ein Sonderfall des allgemeineren Satzes 5. 

Offenbar ist es durch die Verwendung der Naherungslésung v und sogar der 
(im allgemeinen unbekannten) Lésung w in den Ungleichungen (43) und (45) bei 
konkreten Beispielen haufig méglich, zu erheblich giinstigeren Funktionen @ und 
@ als in Satz 4 und damit auch zu giinstigeren Schrankenfunktionen go und @ 
zu gelangen. Dies soll deutlich gemacht werden durch das folgende 
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Beispiel 6. Es soll die Lésung* u(t) €Z(Z) der Volterraschen Integro-Diffe- 
_ rentialgleichung t 
) (1+?) uv’ +3u—sinu : [ cos u(s)ds =—>— (48) 


é(se2)y. 1+ 
0 


unter der Anfangsbedingung w(0)=0 im Intervall 0<t< oo abgeschatzt werden. 


t 
Man wahlt T beliebig groB, Fu= f[ cosu(s)ds und 
0 


f(t, x, y,2) =(1+4) x+3y—siny — z/f(1+ 4) —3/(1+4). 


| Offenbar ist f¢ M und F4¢%,, so % 
-daB eine Abschatzung durch Ober- 
und Unterfunktionen mit den 
Methoden von Teil I nicht un- 
mittelbar méglich ist. Zur Be- 
_stimmung einer Naherungsfunk- 
tion v(t) ersetzt man sin w durch 
wu und cos u durch 1 in (48) und 
findet damit die explizite Nahe- 
rungslosung 


v(t) = 44/(1-+8)2. 


02 


Es ist 
v (0) = 4(0) =0, 01 

fé,v',v, Fv) =d(t) 
mit 
a(t) =v—sinv+ ' 

t kT, Ze RS ES BE CPE) DECLEOR SIO, 
=i Ge (1 +?) abt (1 COS v) ds Abb. 1. Zu Beispiel6: Defekt d(t) und Defektschranke 6(t) 
0 


(vgl. Abb. 1); man kann daher in (41) und (42) etwa e=@=0, —6(t)=6(t) =d(?) 
wahlen. 
Die Ungleichung (43) wird zu 


n,o,...) —f(t,v —o’,...) =4+2) 0 +3@— 2sin 0/2 cos(v — o/2) + 
t 
+ 2#4(1+2)1/f sin o/2sin(v — 9/2)ds= (1+2) 0’ +20, 
0 
falls nur 0=0 (49) 
und 20 —g=0 (50) 


bleibt in 0<t< oo. Diese Eigenschaft lat sich — fiir die noch zu bestimmende 
Funktion 9(t) — spater nachweisen, also kann o (é, x, y, z)=(1-+?)*+2y mit 
w€M gewahlt werden. Damit geht (44) tiber in die Differentialungleichung 
(1+24)0’+20>d(t) mit @(0)>0; eine Lésung ist @(t)=y(1+4)+o(¢) mit 

t 


o(t)=(1+4)*f (1+) d(s) ds (54) 


0 


* Ein Existenzbeweis kann — etwa mit dem Iterationsverfahren — leicht gefiihrt 
werden, die Eindeutigkeit einer Lésung wird in Beispiel 7 nachgewiesen werden. 
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fiir alle Zahlen y>0 (vel. Abb. 2). Wegen d(t)20 (s. Abb. 1) ist o20 und damit 
auch 9=0, also ist (49) erfiillt. Aus d(/)<(1 +12)+ (s. Abb. 1) folgt o(f)S4v(é), 
also auch 9 (t)<2v(¢) auf J, wenn nur y hinreichend klein gewahlt wird; mithin — 
ist auch (50) befriedigt. Fiir y +0 findet man aus (47) damit schlieBlich v (¢) — u(é) 
Soli). 

a Abschatzung von v — wu nach unten betrachtet man die aus (45) folgende 
Ungleichung 


fev +0...) (GU. Ja Une 3eo> 2sin 0/2 cos(v + @/2) + 


mit 
a(t) =2¢7(1+2)+x 


O15 


t 
x f sin 9/2 sin (v — 9/2) ds, 
0 


G10 


GOS 


bleibt in OXt<oo. Auch diese 
Eigenschaften werden sich an- 
| schlieBend nachweisen lassen. 
ais 
0 [i lk oe Ge oe Ge ee INigeyiieian 
Abb. 2. Zu Beispiel 6: Fehlerschranken o(¢) und — G(f) as (t, x, 2) = (1 | #) x4 3y a (0) 


gewahlt werden, die Ungleichung (46) geht damit iiber in die Differentialunglei- 
chun = - eed = 

: (1440 +30—a()>—d(f) mit 9(0)>0; 
eine Lésung ist @(f)=(1+24)y+a(é) fiir alle y>0, wenn 


| 


= t 
—G(t)=(1+4 5 J (1+s)?[d(s) — a(s)] ds 
0 
gilt. Wie man leicht nachpriift, sind die Ungleichungen (52) und (53) fiir alle 
t= 0 erfiillt, also folgt aus (47) fiir y—>0 die Abschatzung —@ (¢)<v(t) — u(t) mit 


—o(é)=(1 Hf (1+ s)?d(s) —2(1+s) of inate sin (v — o/2)dt|ds (54) 
(vgl. Abb. 2). : 


' Die Schranken v(/)—o(¢) und v(¢)+@(é) begrenzen in Abb. 3 einen (schraf- 
fierten) Streifen, der die exakte Lésung u(¢) von (48) enthalt. Man entnimmt 
Abb. 3 das verbliiffende Ergebnis, daB dieser Streifen auferhalb der verwendeten 
Naherungslésung v(t) verlauft, so daB also die Fehlerabschatzung nach Satz 5 
hier sogar zu einer Korrektur dieser Naherung benutzt werden kann, indem 
man etwa v=v- (6 —o)/2 als neue und bessere Naherung fiir w verwendet. Die 
Berechnung des Defektes d(¢) zu 0(t) und neuerliche Fehlerabschatzung fiir 3 
aus (51) und (54) wiirde einen neuen, schmaleren Streifen ergeben, innerhalb 
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dessen die exakte Lésung w (#) gelegen ist ; dieses Vorgehen lieBe sich noch iterieren. 
Ebenfalls lieBe sich eine Verbesserung der Fehlerabschatzung erzielen, wenn man 
die noch ziemlich groben Abschatzungen, die oben zu den Funktionen w und @ 
fiihrten, durch Verwendung der jetzt schon recht genau bekannten Schranken o 
und o verfeinern wiirde. Da in dem vorliegenden Beispiel nur das Prinzip einer 
solchen Fehlerabschatzung gezeigt werden sollte, wurde auf die numerische Durch- 
fihrung dieser Vorschlage verzichtet. 


g7 
hie l 


| | l 
0 FRI D Bi} 4 5 6 7 8 Os SO) 


Abb. 3. Zu Beispiel 6: Naherungslésung v(t) und Schranken v+G bzw. v—o fiir die exakte Lésung 1% (t) 


SchlieBlich ist noch zu sehen, daB sich leicht auch geschlossene Ausdriicke 
aus elementaren Funktionen fiir o und o angeben lassen, z.B. folgt aus 


v—sinv<v3/6, 1—cosv XS v?/2 
d(i) <8 (146048) (14+) *= 6 (55) 
(vgl. Abb. 1) und damit nach (51) die Formel 
o() S$ (21+? (1+4 — (1442-40144 
Padre pre Ol ag Orta Aste 


die Schranke 


(vgl. Abb. 2). 
III. Eindeutigkeitssatze 

Aus den Fehlerabschatzungen des Abschnitts II folgen in vielen Fallen Ein- 
zigkeitssatze. So haben etwa die in den Beispielen 5a) bis 5c) betrachteten 
Anfangswertaufgaben jeweils héchstens eine einzige Loésung, denn unter den 
Lipschitz-Bedingungen (28), (32) und (38) strebt in den Abschatzungen (30), (36) 
und (40) die rechte Seite gegen Null, also strebt v(¢) >w(t), wenn e +0 und 60 
geht. Etwas allgemeiner ist der folgende Eindeutigkeitssatz: 


Satz 6. Die Funktion f(t, x, y, 2) sed in t€Iy, —w< 4, y, 2<+ 0 erklart, 
der Operator F sei €§. Zu y>0 und je zwer Losungen u,vEZ des Anfangswert- 
problems (1), (2) gebe es stets eine Funktion w € M[M] und einen Operator 2€5,[8+] 
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derart, daB is (Z, v' —u',v—u, Q(v — u)) <0 (57) 

ist fiir alle Stellen t€ I, mit v(t) —u(t) =y und v(t)—ult)\<y m O<t<t. Weiter 

existiere eine Funktion 0 €Z, fiir die OS Sy in I, und 

w(t, e’,0,29)>[2]0 ml, (58) 
(0+) > v(0+) —“(0+) (59) 


ist. Dann besitzt das Anfangswertproblem (1), (2) héchstens eine einzige Lésung u€ Z. 
Der Beweis verlauft ganz analog zu demjenigen von Satz 4. 
Bemerkung. Die Voraussetzung (57) ist z.B. erfiillt, wenn f(¢, x, y, 2) einer 
Abschatzung (24) geniigt, denn dann ist fiir zwei Lésungen w, v von (1), (2) stets 


O=f(t,v',v, Fv) —fit,w',u, Fu) Solt,v' — u,v — 4, Q(v—u)). 


und 


Beispiel 7. Als einfache Anwendung von Satz 6 soll zunachst die Einzigkeit 
jeder Lésung der in Beispiel6 betrachteten Integro-Differentialgleichung (48) 
nachgewiesen werden. Offenbar gilt fiir die dort verwendete Funktion 7 und den 
Operator F die Abschatzung 


i(t, x, w, Fw) —f (t,¥, @, F@) = (144) (x—%)+2(w—@) aatsn ft (s)]ds 


fiir beliebige Zahlen x, ¥ und alle Funktionen w, ~€C(J) mit w=w@. Also ist die 

Bedingung (57) erfiillt mit der Funktion @(é, x, y, z)=(11+2)*+2y—2/t(1+2 

und mit dem Operator t 

Qw =f w(s)ds 
0 


dabei ist w¢€ M und Q€E§,. Setzt man nun 9 () =y in J, dann ist w(t, 0’, e, Qo)= 
y(1+22)/(4+2)>0 in J, also befriedigt o(f) die Ungleichung (58). Weiter gilt 
0 (0) =y>0=v(0) — x«(0) fiir je zwei Lésungen wu, v von (48) mit einer beliebigen 
Anfangsbedingung u(0)=y, also ist auch (59) erfiillt. Mithin besitzt (48) unter 
beliebigen Anfangsbedingungen nach Satz 6 héchstens eine einzige Lésung u€Z. 


Beispiel 8 (verallgemeinerte Nagumo-Bedingungen). 
a) Gewéhnliche Differentialgleichung erster Ordnung. Man betrachtet 
u' — g(t, u) =0, u(0)=y, (27) 


also f(t, x, y, 2) =x*—g(t, y), F beliebig. Die reelle Zahl y sei beliebig gewahlt, 
die Funktion g (#, y) sei ee fiir ¢€I, —co< y<-+ oo und sei noch stetig an der 
Stelle =0, y=. Dann gilt fiir jede Lésung w€Z des Anfangswertproblems (27) 
stets u(t) =n-+tg(t, n)+0(t) fiir ¢-0, also ist dann 


v()—u(t)=o() fir t>0 (60) 


fiir je zwei Loésungen u,v von (27). Weiter geniige die Funktion g der verall- 
gemeinerten Nagumo-Bedingung*. 


st») — 89) <|$4+L0] (y 9) (64) 
* Diese Bedingung verdanke ich Herrn Doz. Dr. W. WaLTER. 
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fiir ¢¢ J) und alle y=. Dabei sei die Funktion L(é) in J, stetig und das Integral 


fig ds 


existiere (eventuell uneigentlich) fiir #¢ J. Dann geniigt die Funktion o (f, x, y,z) = 


| x— (+. +L (6) y der Ungleichung (57) fiir einen beliebigen Operator 2€%, und 


es istwCM. Ist y>O0 fest gewahlt, so liegt die Funktion 


t 
o (t) =Bt(1+2) exp f L(s) ds 
0 
fiir hinreichend kleine Zahlen f>0stetsin0XoXy. Weiter ist 9 €Z und geniigt 
wegen 


t 
| altos} op..) =ptexpy L(syds> 07 inG 
6 


der Bedingung (58). SchlieBlich ist noch o(f)=ft+o(é) fiir t+0, also ist 

o (0+) >v(0+) —u(0-+) nach (60), womit auch Ungleichung (59) befriedigt ist. 
Damit hat das Anfangswertproblem (27) unter der Bedingung (61) héchstens* 
eine einzige Lésung u€Z. 

In dieser Aussage ist der oben erwahnte, aus Beispiel 5a) folgende Eindeutig- 
keitssatz nicht enthalten, obwohl natiirlich die Bedingung (61) auch die Lipschitz- 
Bedingung (28) mit umfaft. In Beispiel5a) war namlich tiber die Funktion 
g(é, y) nur vorausgesetzt, daB g fiir ¢€I,, —coo<y<oo definiert sein sollte, 
wahrend beim vorliegenden Satz g noch bei ¢=0 definiert und an der Stelle 
t=0, y= sogar noch stetig sein soll. 

b) Volterrasche Integralgleichung zweiter Art. Man betrachtet die nichtlineare 


Integralgleichung 
t 


UA Dios PAGS, 22(6)) ds=g(t). (31) 
Es sel 


f(t, %,¥,2)=y—z-g(), fi tee ae 


gewahlt. Die Kernfunktion K(s, ¢, w) sei integrierbar fiir s,t€J, we C(I); die 
Funktion g(t) sei stetig auf J und die vordere Ableitung d, g(é)/dt existiere bei 
t=0. Weiter geniige K der Bedingung 

K(s,t, y) — K(s, 4,9) S M(s) N(@) (y —¥) (62) 


fiir s,¢€ I, und y=y, wobei die Funktion M(s)=0 und auf I, definiert, N(¢)>0 
und auf J stetig sein soll; und in J, gelte noch 


MN) S++L@ (63) 


mit einer in J, stetigen, auf J (eventuell uneigentlich) integrierbaren Funktion 
L(t). L=0 entspricht offenbar der iiblichen Nagumo-Bedingung. 


* Unter den genannten geringen Voraussetzungen an die Funktion g ist die 
Existenz einer Lésung von (27) natiirlich nicht allgemein gesichert. 
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Zur Anwendung des Satzes 6 setzt man dann 
t 


o(t,%,,2)=y—z mit w€M und Qu = N(t) f [-+L()]u(s)as. 
0 
Der Operator Q ist zwar monoton, wegen des Terms 1/s im Integranden jedoch 
nicht auf alle Funktionen w€ C(I) anwendbar. Immerhin existiert jedoch Qw 
sicher fiir alle Funktionen w€ C(I) mit w(t) =O (t) fir +0. Offenbar gentigt es 
nun schon zum Beweis von Satz 6, wenn  allein definiert ist fiir die Differenz 
» —u zweier Lésungen u, v von (31) und fiir die gesuchte Funktion @ (¢). Wegen 


u(t) =g(0) + ¢[g'.(0) + K (0, 0, g(0))] +0 (4) 
ist jedoch 
v(t) —u(t)=o() fir t+0 (64) 


fiir die Differenz zweier beliebiger Lésungen u, v€Z von (31), also ist Q(v—w) 
auf I definiert und stetig. Wegen (62) und (63) ist dann die Ungleichung (57) 
erfiillt fiir alle Lésungen w, v€Z von (31). 


Man setzt nun 
t 
o (t) =Bt N(é) exp f L(s) ds. 
0 


Ist B>0 hinreichend klein und y>0 gewahlt, so gilt OXo (4) <y fir alle ¢cJ. 
Wegen o(t)=ft+o(t) ist 0(0+)>v(0+)—w(0+) nach (64), also befriedigt @ 
die Bedingung (59). Weiter ist damit der Operator 2 anwendbar auf 9 und es gilt 


t Ss 
wl,«..;0,26)=BNO Ss exp PL@) drds > 0. da, 
0 0 


d.h. auch (58) ist erfiillt. Damit hat die Integralgleichung (31) unter der Be- 
dingung (62), (63) héchstens eine einzige Lésung uw €C(Z). 

Die Einzigkeitsbedingung (62), (63) findet sich fiir V =1 und L =0 bei T. Sato 
[6], S. 282, allgemein ist sie in dessen Théoréme 13 enthalten. Auch hier ist der 
aus Beispiel $b) flieSende Eindeutigkeitssatz nicht in dem soeben abgeleiteten 
Ergebnis enthalten, da in Beispiel 5b) die Voraussetzungen iiber K und } milder 
sind. 


c) Integro-Differentialgleichungen. Gegeben sei die Integro-Differentialglei- 
chung 


T 
w'(t) — g(t, w(t)) — f K(s, t, u(s)) ds =0 (37) 
0 
mit der Anfangsbedingung w(0)=y7. Es sei 
2 i 
f(t, %, y,2) =x—g(t,y)—z und Fu = | K(s,t, u(s)) ds 
i) 


gewahlt; es sel F€%, das ist z.B. der Fall, wenn die Kernfunktion K(s, t, u) 
integrierbar fiir 8, t€L , wEC(L) ist. Die Funktion g(é, y) sei definiert fir ¢€ J, 
— =< y< oe, stetig an der Stelle =0, yn und geniige der Lipschitz-Bedingung 


8.9) —g”) SLO (y—y) fir t€l,, yey (28) 
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mit einer in J, stetigen Funktion L(é). Weiter gelte 


- t 

J [K(s,t, w(s)) — K(s, t, @(s))] ds < N( f M(s) [w(s):— (s)]ds (65) 
0 0 

fiir alle ¢¢J und alle Funktionen w, w€C(Z), fiir die w= auf J gilt. Dabei 
sel N, MEC(,) mit N20, MZO auf J und fiir alle ¢€ J existiere das Integral 


t 
Js M(s)ds_ (eventuell uneigentlich). 
0 


SchlieBlich sei noch 
N() +¢L@) S1+2M(t) in Ig. (66) 


Die Bedingung (66) ist z. B. erfiillt, wenn L und N in J, je einer Abschatzung 
L(t) SA/t+1, N(t) Sv +tn geniigen, wobei die Konstanten 4,/,y, n die Ungleichungen 
A+vs1, l+nsStM(t) in I befriedigen. Ist y=J=0, so ist dies die Nagumo- 
Bedingung fiir L(é). 

Man setzt nun w(é, x, y, 2) =x—L()y—N(é)z in J, und 


Qw = f M(s) w(9) ds. 


Wegen (28) und (65) ist dann (57) erfiillt. Es ist a € M, der Operator Q2 ist mono- 
ton, jedoch ist 2 wie im Falle b) unter Umstanden nicht auf alle Funktionen 
w€C(I) anwendbar. Genau wie dort tiberlegt man sich, daB es ausreicht, wenn 
{2 allein definiert ist fiir die Differenz v—wu zweier Lésungen wu, v€Z von (37) 
und fiir die gesuchte Funktion 9. 


Es ist 
i he 


u(t) =n +t|g(0,) + S K(s,0,) ds} +0(), 
0 
also gilt 
v(t) —u(t)=o(t) fir ¢t-0, (67) 
mithin ist Q(v—) auf J definiert und stetig. 


Man setzt nun 
t 
o(t) =Bt(1+4) exp f s(1+s) M(s) ds 
0 
mit einer hinreichend kleinen Zahl B>0. Es ist @€Z und o(t)=ft+o/(t), also 


gilt 0(0+)>v(0+)—u(0+) nach (67), folglich ist die Bedingung (59) erfiillt. 
Mit (66) und wegen M=O, N2ZO ist weiter 


Bie 2 ene 1 — NY pee (te OE) MEN) BN 


also befriedigt @ auch noch (58). Damit besitzt das Anfangswertproblem (37), (2) 
unter den Bedingungen (28), (65), (66) héchstens eine einzige Loésung u€Z. 

Analog zu oben ist in diesen Einzigkeitskriterien das Ergebnis von Beispiel 5 c) 
nicht als Sonderfall enthalten. Dagegen erhalt man fiir V =0 den Eindeutigkeits- 
satz von Beispiel 8a). 
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[Z] 
[2] 


be % 
Kart NickeL: Integro-Differentialgleichungen 
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